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ABSTRACT 

-  ^  An  experimental  apparatus  to  measure  the  particle  density 
distribution  in  the  freeboard  of  an  atmospheric  fluidized  bed  uas 
designed  and  constructed.  The  density  versus  height  measured  by 
the  sampling  apparatus  gives  a  similar  exponential  decrease  as 
previous  investigations  have  found. 

A  particle  trajectory  model  is  developed  which  calculates  the 
height  and  particle  density  distributions  above  the  bed  surface  of 
an  atmospheric  fluidized  bed.  The  parameters  input  to  the  model 
are  the  superficial  velocity,  initial  particaC  velocity,  gas  jet 
velocity  and  duration,  and  the  particle  size  distribution  of  the 
bed  mass.  The  model  uas  evaluated  using  the  experimental  data  for 
jet  velocity,  deration,  and  particle  size.  The  predicted  slope  of 
the  particle  density  versus  height  in  the  freeboard  agrees  with 
the  experimentally  measured  slope  within  Z0X. 

A  sensitvity  anal ysi  seising  the  trajectory  model,  resulted  in 
a  determination  of  the  particle  distributions  in  the  freeboard  of 
a  fluidized  bed  as  affected  by  varying  the  input  parameters  to  the 
trajectory  model.  The  most  significant  effects  were  achieved  when 
the  jet  velocity  or  duration  uas  altered.  (  7/,^,  ..  ^ 
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CHAPTER  I 


INTRODUCTION 

Fluidized  beds  have  been  used  in  industry  for  many  years. 
They  have  been  used  to  nix  and  dry  particulate  materials  and  are 
the  principle  process  in  catalytic  cracking  plants.  In  the  past 
decade  or  so,  the  use  of  fluidized  bed  conbustors  for  power 
generation  has  become  a  source  of  major  interest.  Prototype  coal 
burning  beds  have  already  been  built  which  are  comparable  to 
existing  coal  plants.  Fluidized  bed  conbustors  have  the  added 
benifit  of  low  N Ox,  SOj  and  hydrocarbon  emissions  and  the 
flexibility  of  being  able  to  burn  a  wide  range  of  fuels  ranging 
from  refuse  and  high  sulfur  content  coal  to  high  grade  fuels. 

A  fluidized  bed  [Fig  1]  is  composed  of  a  distributor  through 
which  an  air  flow  is  Introduced  through  thousands  of  small 
orfices.  This  air  then  passes  through  the  dense  zone  of  the  bed 
which  is  comprised  of  a  mass  of  particles.  The  air  velocity 
through  the  dense  zone  is  maintained  above  the  minimum 
fluidization  velocity  (Unf)  during  normal  operations.  At 
velocities  equal  to  or  greater  than  Unf,  the  frictional  force 
(Drag)  of  the  air  flowing  past  a  particle  is  equal  to  the  weight 
of  the  particle.  Under  these  conditions,  the  particle  mass 
behaves  very  much  like  a  fluid.  It  will  maintain  a  horizontal 
surface  if  the  container  is  tilted,  flow  out  of  holes  in  the 


container,  and  has  a  pressure  drop  across  any  section  of  the  bed 
apprxanately  equal  to  the  weight  of  air  and  particles  in  the 
section  Ill.  At  velocities  above  Umf,  the  excess  air  uill 
coalesce  and  form  snail  voids  or  bubbles  as  it  procedes  towards 
the  surface  of  the  fluidized  bed. 

As  bubbles  break  at  the  surface  of  the  bed,  the  solid 
particles  are  throun  up  above  the  bed  surface  and  are  entrained  by 
the  upward  flowing  gas  strean.  This  zone  above  the  bed  surface  is 
the  freeboard  zone.  In  the  freeboard,  sone  particles  are  carried 
by  the  gas  flow  far  above  the  bed  surface  and  are  renoved  fron  the 
fluidized  bed  (elutriated)  ,  while  the  reminder  fall  back  to  the 
bed.  In  general  the  an aunt  of  bed  solids  suspended  in  the 
freeboard  ( entrainnent >  decreases  exponentially  along  the  freeboard 
height.  This  distribution  is  sinilar  to  that  of  the 
Maxwell -Boltzmann  distribution  for  the  case  of  discrete  energy 
states  as  it  applies  to  the  Law  of  Atmospheres  IZ1. 


N(z)  =  No  exp<-mgz/kT> 


(1) 


Research  in  the  area  of  entrainment  by  Lewis  et  al.  (31,  Zenz 
and  Ueil  14]  and  others  has  resulted  in  the  following  correlation 
for  entrainment  as  a  function  of  gas  velocity  and  freeboard  height 


for  small -particle  beds  I  IK 


At  Uo 


2 

B  exp<  *  I ( b/Uo)  +  a  H  ]> 


(Z) 


where: 

F  =  Entrainnent 
At  =  Area  of  bed 

Uo  =  Superficial  velocity  through  bed 
B  =  Particle  dependent  constant 
b  =  Particle  dependent  constant 
a  =  Particle  dependent  constant 
H  =  Height  of  freeboard 


At  sone  point  above  the  bed  surface,  the  quantity  of 
entrained  particles,  becones  .constant.  At  this  point,  the  free 
fall  velocity  of  the  renaining  particles  is  equal  to  or  less  than 
the  uniforn  superficial  operating  velocity.  The  height  at  which 
the  entrainnent  becones  constant  is  called  the  transport 
disengagenent  height  (TOH)  (11. 

The  particles  that  are  thrown  above  the  bed  are  affected  in 
the  freeboard  region  by  hydrodynanic  paraneters  such  as:  bubble 
size,  bubble  frequency,  fluidizing  velocity,  height  above  the  bed, 
particle  size,  particle  density,  column  diameter  (SI,  and 
baffles  (Gl.  The  internittent  high  velocity  bursts  of  gas  which 
occurs  when  a  bubble  bursts,  inposes  a  flu ousting  and  highly 
irregular  tine  dependent  velocity  profile  over  the  cross  section 
of  the  bed  surface.  At  successively  higher  levels  above  the  bed 
surface,  this  velocity  profile  becones  nore  and  more  uniforn  until 


at  the  TDH,  the  flow  is  at  the  uniforn  superficial  operating 


Until  recently,  little  attention  was  paid  to  the 


understanding  of  the  freeboard  reactions  for  large  particle  beds. 
However,  due  to  recent  work  in  fluidized  bed  combustors,  the 
extent  of  particle  and  fines  loading  in  the  freeboard  has  been 
shown  to  significantly  affect  the  SO2  absorbtion,  NOx  reduction,  CO 
emission.  In  general,  the  last  5  to  10X  of  the  combustibles  will 
burn  in  the  freeboard.  It  was  shown  that  the  fine  sorbent 
particles  entrained  into  the  freeboard  will  enhance  sulphur 
capture  and  that  entrained  char  particles  will  react  with  NOx  and 
reduce  its  emission  [71.  _  Fines  reinjection  has  been  shown  to 
significantly  Increases  the  fine  particle  concentration  in  the  bed 
and  in  the  freeboard  with  the  conscience  of  further  enhancing 
char  oxidation.  However,  this  can  result  in  overheating  in  the 
freeboard  region  and  excess  SO2  and  NOx  emission  [71.  The  end 
result  is  that  the  potential  for  fluidized  bed  combustion  power 
plants  is  enhanced  by  their  ability  to  burn  high  sulphur  content 
fuels  and  maintain  low  SO^  emissions.  Further  research  in  the 
area  directly  above  the  bed  surface  is  required  to  properly  model 
the  reactions  occur i ng  within  the  freeboard. 

Extensive  studies  on  entrainment  rate  and  elutriation  have 
been  made  with  numerous  correlations,  some  of  which  are  proposed 
in  Cl,  3,  4,  5,  S,  7,  81.  However,  most  of  the  reported  work  on 
entrainment  from  fluidized  beds  has  been  carried  out  with  either 
a  closely  sized  fraction  of  particles  or  a  mixture  of  two  such 
fractions.  Virtually  all  of  this  work  has  been  conducted  on  bench 
scale  or  catalytic  cracking  fluidized  beds.  The  results  of  these 
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studies  have  been  shoun  to  produce  much  lower  entrainment  rates 
than  full  sized  beds  or  tend  to  operate  in  the  slugging  condition. 
Therefore,  the  entrainment  rates  and  transport  disengagement 
heights  (TOH)  for  fluidized  beds  are  generally  estimated  from 
empirical  or  seniempirical  correlations  obtained  from  this  data 
and  most  of  then  show  extreme  discrepancies  between  different 
experimental  results.  Extrapolation  of  these  empirical 
correlations  usually  leads  to  strange  results  151  with 
discrepancies  which  can  vary  by  two  orders  of  magnitude. 

*  —  -  -The  lack,  of  good: correlations  stems  mainly  from  the 
difficulty  in  obtaining  accurate  entrainment  rate  data,  host  of 
the  data  is  based  on  pressure  measurements  at  incremental  heights 
in  fluidized  cracking  plant  type  beds  (111.  The  effect  of  wall 
loading  by  particles  and  the  actual  relation  between  pressure  and 
particle  concentration  is  considered  to  be  major  sources  of  error 
when  using  this  method  with  large  particles.  As  a  result,  none  of 
the  correlations  are  widely  accepted  as  giving  accurate 
predictions  (81. 

A  complete  model  of  the  entrainment  process  from  fluidized 
beds  must  take  into  account  all  the  mechanisims  involved  within 
the  process.  The  arrival  of  bubbles  at  the  bed  surface,  ejection 
of  particles  from  the  dense-phased  bed  into  the  freeboard  region 
as  the  bubbles  erupt,  particle-particle  interactions,  and  the 
trajectories  of  ejected  particles  are  all  important  181.  Much 
work  has  been  done  concerning  bubble  growth,  velocity,  volume, 


etc.  and  their  behavior  is  fairly  uell  understood. 

The  nechanisn  of  solids  ejection  at  the  bubbling  bed  surface 
is  still  not  uell  understood.  The  origin  of  ejected  particles  is 
reported  to  be  primarily  due  to  iuo  sources.  The  particles  uhich 
have,  been  lifted  by  the  bubble  uake  and  throun  upwards  following 
the  bubble,  burst  _at  the_sutf  ace  is^  the  _f  irst  source.  This  theory 
is  supported  by  work  done  by  6eorge  and  6race  (81  who  performed 
experiments  uhich  concluded  that  the  vast  majority  of  the  ejected 
particles  did  not  originate  from  the  surface  layers  but  from 
bubble  wake  pick  up.  Uork  done~by  Page  and  Harrison  (61  also 
appears  to  agree  with  this.  The  second  theory  suggests  that  the 
ejected  particles  originate  at  the  nose  of  the  bursting  bubbles 
and  are  throun  outward  uhen  the  bubble  breaks.  Research  by  Roue 
and  Partridge  [8]  and  Glicksman  et  al  tlZl  have  shown  this  second 
mechanisim  as  being  the  dominate  particle  ejection  source  and  thus 
supporting  this  second  theory.  Their  uork  has  also  shoun  that 
under  the  conditions  in  uhich  Z  bubbles  coalesce  just  below  the 
surface  of  the  bed,  the  jet  of  gas  produced  can  result  in  a 
significant  amount  of  particles  being  ejected  from  the  uake  of 
the  first  bubble. 

The  effect  of  multiparticle  interactions  have  been  for  the 
most  part  ignored  except  by  Peters  and  Prybylouski  Cl 31.  The 
motion  of  any  individual  particle  is  influenced  by  the  presence 
of  other  particles,  i.e.,  through  direct  particle-particle 
interactions  and  deviations  in  the  fluid  drag  force.  The  major 


drawback  of  their  work  is  that  the  paper  compares  their  theory 
with  only  a  single  set  of  experimental  results  (3). 

Studies  to  model  the  trajectories  of  particles  in  the 
freeboard  have  been  conducted  several  tines.  The  work  of  Ualsh 
etalt71 ,  Georgeand  Grace  [81,  and  Peters  and  Prybylouski  [131 
name:  just' a  feu  ofi  the' latest  efforts.  All  of  these  studies 
relied  upon  experimental  data  to  develope  their  theories. 

However,  to  check  the  accuracy  of  their  theories,  more 
experimental  data  is  required. 

As  of  yet,  none  of  the  entrainment  models  available  can  be 
incorporated  into  fluidized  bed  combustion  models  with  sufficient 
accuracy  to  warrant  their  use.  This  is  due  to  a  lack  of 
experimental  information  on  entrainment  rate  as  a  function  of  the 
complete  fluidization  parameters  of  the  bed  to  test  the  models 
with.  As  a  result,  the  purpose  of  this  study  has  been  to  obtain 
particle  density  distributions  above  a  cold  atmospheric  fluidized 
bed  containing  a  continuous  particle  size  distribution  [Appendix 


There  are  many  methods  available  for  determining  particle 


distributions  in  fluid  flous.  The  more  commonly  used  methods  are: 

1>  Catching  mechanisms 
-r sirchr-r-i -arm  Z)  Trapping  mechanisms 

3)  Radiation  atteraiation  measurements 

4)  Optical  measurements 

5)  Capacitance  and  Inductance  measurements 

Catching  mechanisms  are  passive  devices.  That  is,  particles 
are  captured  merely  by  the  presence  of  the  catching  mechanism  in 
the  fluid  flou  containg  the  particles  to  be  sampled.  The  data 
obtained  using  this  method  is  position  dependent  and  produces 
average  values  for  the  particle  flux  loadings.  These  catching 
mechanisims  are  also  limited  in  that  they  can  only  catch  particles 
uith  particle  fluxes  traveling  in  a  single  direction.  The  device 
used  by  Ualsh  et  al  [101  only  captured  falling  particles  while  the 
device  used  by  George  and  Grace  C 81  required  the  upward  moving 
particles  to  deflect  off  of  a  baffle  surface  and  fall  into  a 


collecting  trough. 


Trapping  nachanisins,  unlike  catching  nechanism,  are  active 
particle  samplers.  Their  operation  involves  the  trapping  and 
isolation  of  a  finite  volume  of  the  fluid  flou  at  a  specific 
period  in  time.  This  sampling  technique  produces  time  dependent 
as  uell  as  position  dependent  data.  This  uill  allow  correlations 
between. bubble^eruption^and  particle. density  to  be  made  using 
multiple  bubble  conditions  rather  than  single  bubble  capture.  As 
the  number  of  random  samples  taken  by  this  method  increases,  the 
average  value  of  this  data  uill  approach  that  of  the  catching 
nechanism.  Trapping  mechanisims  also  capture  particle  fluxes 
traveling  in'multiple"  direetionsr~This~  ability  reduces  the  error 
inherent  in  measuring  only  the  downward  or  only  the  upuard 
particle  flux.  The  apparatus  used  in  this  paper  is  a  trapping 
mechanism. 

Attenuation  of  nuclear  particles  from  a  radioactve  source  can 
be  used  to  give  average  particle  density  distributions  across  a 
suspension.  However,  this  method  is  not  adaptable  to  density 
determinations  at  a  point.  An  average  time  dependent  density 
determination  can  be  achieved  with  this  method.  Another  drawback 
of  this  method  is  the  radiation  hazards  involved  with  the  use  of 
nuclear  particles. 


Optical  density  determinations  consists  of  two  seperate 
methods.  The  first  method  uses  a  very  small  light  beam  which  is 
eclipsed  by  the  transition  of  a  particle  through  it.  A  related 


nethod  uses  the  absorbtion  and  scattering  of  a  somewhat  larger 
light  bean  to  correlate  the  change  in  light  intensity  with 
particle  density.  This  method  has  been  used  frequently  in  the 
study  of  aerosols  but  requires  complicated  and  intricate  equipment 
[141.  The  second  method  involves  high  speed  photographs  of  a 
small  volume  of  space.  This  method  cannot  be  used  uhen  the 
particle  density  is  so  large  that  multiple  particles  eclipse  each 
other  frequently  enough  to  produce  unacceptable  error.  This  is 
the  case  uhen  the  probe  height  above  the  bed  is  less  than  7-15  cm 
(3-6  in). 

....--.Measurements  at  a  point  can  also  be  made  by  inserting  either 
a  toroidial  inductor  or  a  parallel  plate  capacitor  in  the  flou. 

The  presence  of  the  particles  changes  the  permeability  and  thus 
the  inductance  of  the  inductor,  or  the  dielectric  strength  and 
thus  the  capacitance  of  the  capacitor.  The  drau-backs  of  these 
methods  involves  the  unknoun  effects  of  particle  velocity  and 
external  particles  on  the  inductor  and  charge  transfer  to 
particles  from  the  capacitor  [141. 

Apparatus  Requirements 

The  goal  of  this  study  uas  to  determine  the  density 
distribution  of  particles  above  an  atmospheric  fluidized  bed  with 
particle  velocities  of  up  to  10  meters  per  second.  The  particle 
size  distribution  of  material  ejected  from  bubbles  is  required  for 


particle  trajectory  calculations.  A  correlation  between  the 


average  density  and  the  density  present  immediately  after  a  bubble 
bursts  from  the  bed  surface  uas  also  of  interest.  These 
requirements  dictated  that  the  method  used  for  measuring  densities 
have  the  following  capabilities: 

1)  Measure  densities  with  good  spatial  resolution. 

Z)  Measure  densities  at  specific  moments  in  time. 

_^3)T  Obtain,  particle  size  information. 

4)  Operate  under  extremely  dirty  conditions. 

5)  Easy  sample  removal  from  bed. 

S)  Remote  operation  of  sampler. 

The  radiation  attenuation  and  inductance/capacitance  methods 
can  not  determine  particle  sizes.  Therefore,  these  methods  uere 
no  longer  considered  as  possible  measurement  alternatives. 

Because  the  optical  methods  are  either  not  reliable  at  small 
heights  above  the  bed  or  their  use  is  too  complex,  they  uere  not 
used.  Catching  devices,  although  simple  to  use,  do  not  have  the 
ability  to  measure  data  at  specific  points  in  time  and  determine 
particle  density  loading  in  space.  As  a  result,  the  determination 
to  use  a  trapping  mechanism  as  the  method  of  measuring  particle 
densities  uas  made. 
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General  Design  Criteria 


The  following  criteria  was  used  in  determining  the  design  of 

i 

the  trapping  device. 

1)  The  closure  tine  of  the  trap  was  chosen  to  be  equal 
to  the  tine  for  a  particle  with  a  velocity  of  10 
— _l._l.'~_n/s'  to.  transit-1/10- the  length  of  the  sample 

container.  This  velocity ■ is  considered  to  be  the 
upper  Unit  of  the  particle  velocity  distribution 
present  in  the  test  bed,  based  on  the  work  of 
...  George  and  6race  (81. 

Z)  The  apparatus  must  be  capable  of  frequent  sampling 
without  requiring  access  to  the  sampling  device 
itself. 

3)  The  samples  trapped,  must  be  easily  accessable 
from  outside  the  fluidized  bed  uithout  interupiing 
the  bed  conditions. 

4)  The  apparatus  must  be  able  to  operate  in  the  high 
particle  flux  environment  of  the  fluidized  bed. 


5)  The  actuation  of  the  trapping  device  rust  be  able 
to  be  accurately  determined  to  allou  correlation 
uith  other  tine  resolved  measurements. 

General  Design 

The  apparatus  is  shown  in  Figs.  Z  and  3.  A  description  and 
list  of  all  components  is  given  in  Appendix  H.  The  sample 
container  is  mounted  on  an  extension  arm  to  minimize  the 
disturbance  to  the  air  flow  around  the  sample  trap  caused  by  the 
rest  of  the  mechanism.  The  sample  container  is  closed  using  tuo 
(Z)  paddle  arms,  one  above  and  the  other  below.  These  paddle  arms 
are  attached  to  aluminum  interfaces  which  are  used  to  connect  them 
to  a  rotary  solenoid.  The  solenoid  is  used  to  suing  the  paddles 
over  the  sample  trap  and  shut  it.  Not  shown  in  these  figures  are 
the  power  supply  for  the  solenoid,  the  vacuum  system  used  to 
remove  the  particles  from  the  sample  trap  and  the  water-proof 
nylon  shell  used  to  keep  the  particles  from  interfering  uith  the 
operation  of  the  solenoid.  All  of  these  systems  are  described  in 
greater  detail  in  the  following  sections. 

Sample  Container 

To  ensure  that  the  sampling  device  had  minimal  effect  on  the 
fluid  flow,  the  cross  sectional  area  presented  to  the  flow  had  to 
be  minimized.  This  constraint  required  that  the  sample  container 


Extension  Bar 


be  separated  f roe  the  rest  of  the  apparatus.  This  uas  benificiai 
in  the  final  design  because  it  helped  to  reduce  the  apparatus 
closure  tine. 

To  ensure  a  good  seal  uas  achieved  uhen  the  trap  uas  shut, 

i 

felt  uas  used  as  a  gasket  betueen  the  sample  container  and  the 
closure  paddles^!  Fig.  4  shous  the  final  design  used  for  the 
sampling  nechanisn.  It  is  constructed  of  1/1S  inch  aluminum  uith 
1/4  inch  square  stock  used  for  the  frame  and  mounting  structure. 
Epoxy  is  used  to  seal  the  sides  of  the  container. 

Closure  Paddles 

After  several  iterations  on  paddle  designs,  it  uas  determined 
that  the  paddle  construction  uhich  offered  the  greatest  stiffness 
for  the  least  ueight  uas  a  composite  laminate.  The  paddle,  shoun 
in  Fig.  5,  is  made  using  a  0.4  in  thick  foam  core  uith  1/3Z  inch 
thick  8assuood  laminations  on  the  exterior.  Harduood  (Maple)  end 
pieces  uere  used  to  provide  a  noncompress ive  connection  betueen 
the  foam  paddles  and  the  aluminum  interfaces.  The  aluminum 
interfaces  couple  the  solenoid  shaft  to  the  paddles.  Epoxy  uas 
used  to  join  the  laminate  materials. 

Actuator 

A  rotating  mechanism  utilizing  a  rotary  solenoid  uas  chosen 
to  shut  the  sample  trap.  A  rotary  solenoid  uas  selected  because 


Design  of  Closure  Paddles.  Dimensions  in  cn. 


of  the  problem  of  Maintaining  lou  friction  surfaces  for  sliding 
mechanisms  in  the  presence  of  the  particulate  natter  being 
sampled.  This  precluded  the  use  of  any  sliding  nechanisn  to  shut 
the  sample  trap.  The  calculations  found  in  Appendix  A  and 
Appendix  C  determined  the  size  of  the  solenoid  required  to  achieve 
the  desired  closure  tine. 

The  45  degree  stroke  solenoid  uas  chosen  to  place  the  paddle 
arms  far  enough  auay  from  the  sample  trap,  such  that  uhen 
de-energized  it  prevents  interference  uith  the  particle  flou. 

This  stroke  also  minimized  the  area  which  must  be  clear  of 
obstructions  to- the  traveLrOf-the  paddles.  The  solenoid  operates 
at  a  1/10  duty  cycle  pouer  rating  uhen  initially  actuated, 
providing  the  torcpje  output  shown  in  Fig.  6.  After  the  solenoid 
has  shut  the  sample  trap*  the  solenoid  is  operated  at  a  lower 
power  rating,  providing  a  holding  torque  of  5.5  in-lbs.  This 
decreased  rating  is  necessary  to  prevent  overheating  of  the 
solenoid.  This  assembly  is  encased  within  a  nylon  shell. 

Pouer  Supply 

Fig.  7  is  a  schematic  of  the  electrical  system  used  to  power 
the  solenoid.  Appendix  H  contains  a  list  of  all  components  used 
in  the  power  supply.  The  power  supply  plugs  directly  into  a 
standard  115  volt  AC  line  source.  Switch  SI  is  used  to  apply 
pouer  to  the  solenoid  Ml.  The  full  wave  bridge  rectifier  assembly 
converts  the  AC  line  voltage  to  DC.  The  rectifier  assembly  is 


protected  by  the  arc  suppressor  01  fro*  the  large  voltage  spike 
induced  in  the  solenoid  coil  uhen  the  solenoid  is  actuated. 


Resistor  R3  is  used  to  reduce  the  current  through  the  solenoid 
after  the  solenoid  has  been  closed  to  prevent  it  fro*  overheating. 
To  initially  shut  the  sa*ple  trap,  full  current  is  applied  to  the 
solenoid  by  shutting  switch  SZ,  which  bypasses  resistor  R3. 
Resistors  R1  and  RZ  for*  a  voltage  divider  network  to  provide  a 
low  voltage  (5.5  volts  AC)  trigger  source  for  an  oscilloscope. 
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Vacuum  System 


To  remove  the  particles  in  the  sample  trap,  a  vacuum  system 
was  developed.  This  system  is  shoun  in  Fig.  8  uith  a  list  of  the 
components  given  in  Appendix  H.  A  vacuum  is  produced  by  allouing 
air  from  a  100  psi  air  source  to  flow  through  valve  VI  into  the 
venturi  eductor  PI.  The  vacuum  places  a  suction  on  the  sample 
container  CZ  via  a  fine  mesh  screen.  The  purpose  of  this  screen 
is  to  prevent  particles  from  escaping  the  sample  container.  The 
suction  is  applied  to  the  sample  trap  Cl  through  1/4  inch  polyflou 
tubing.  It  is  through  this  tube  that  the  particles  are  removed 
from  the  sample  trap  and  collected  in  the  sample  container.  An 
equalization  and  agitaion  line  is  connected  to  the  opposite  side 
of  the  sample  trap.  This  line  serves  tuo  purposes.  First,  it 
ensures  that  the  vacuum  system  does  not  pull  in  particles  From 
outside  of  the  sample  trap.  Second,  it  allous  a  flow  of  air  to  be 
introduced  uhich  stirs  up  the  particles  trapped  inside.  This 
helps  push  them  into  the  suction  line  and  reduce  the  remaining 
particles  to  a  minimum. 

Apparatus . T??ting 

Two  tests  uere  run  to  determine  the  effectiveness  of  the 
system.  The  first  test  determined  the  closing  time  of  the  sample 
trap.  The  second  evaluated  the  error  from  the  loss  of  particles 


uhich  were  left  in  the  sample  trap  by  the  vacuum  system. 


Schematic  Diagram  For  Vacuum  System. 

(J  Reference  number  far  component 
list  in  Appendix  H. 


The  procedure  and  method  used  to  determine  closure  time  for 
the  sample  trap  is  given  in  Appendix  C.  From  the  results  of  these 
tests  given  in  Appendix  C,  the  closure  time  uas  determined  to  have 
an  average  value  of  1.44  milliseconds.  This  is  equal  to  a  10  m/s 
particle  traveling  1.44  cm  or  approximately  IS  X  of  the  sample 
trap  length.  The  average  particle  velocity  uill  be  less  than  Z 
m/s  and  uill  therefore  introduce  an  average  error  of  less  than  3  X 

Using  the  sample  trap  closure  time,  a  dynamic  analysis 
determined  that  the  total  time  from  initially  applying  power  to 
the  solenoid  until  it  shut  the  sample  trap  is  4Z.G  milliseconds. 
These  calculations  are  given  in  Appendix  □.  This  actuation  time 
is  important  for  determining  whether  a  specific  bubble's  debris 
uas  uithin  the  vicinity  of  the  trap  at  the  time  of  closure. 

The  procedure  and  results  for  the  testing  of  the  sample 
removal  vacuum  system  are  given  in  Appendix  8.  The  average  amount 
of  particles  lost  by  the  vacuum  system  uas  determined  to  be  0.5Z  X 
of  the  initial  sample  placed  in  the  trap.  The  maximum  error  uas 
0.33  X.  The  error  from  the  vacuum  system  is  therefore  considered 
to  have  an  insignificant  effect  on  the  data  obtained. 


CHAPTER  III 


EXPERIMENTAL  PROCEDURE 

Fluidized  Bed  Configuration 

The  M.I.T.  ainospheric  fluidized  bed,  in  uhich  the  sampling 
device  uas  used,  is  a  model  of  the  20  MU  atmospheric  fluidized  bed 
combustor  prototype,  jointly  sponsored  by  the  Tennessee  Valley 
Authority  and  the  Electric  Pouer  Research  Institute.  The 
fluidized  bed  model  is  described  in  Lord  et  al  [15]  and  Jones  et 
al  C1E3.  Uhile  using  the  sampling  apparatus,  a  different  heat 
exchanger  tube  bundle  configuration  uas  used  than  is  described  in 
ones  et  al  CIS].  The  heat  exchanger  configuration  used  is  shoun 
in  Figs.  9  and  10. 

The  heat  exchanger  used  during  this  uork  is  made  of  1.Z6  cm 
(0.5  in)  0.0.  tubing  arranged  in  4  rows  of  ZZ  pipes  each.  The 
tubes  are  aligned  as  shoun  in  Fig.  10.  Each  pipe  is  spaced  uith 
a  vertical  center  to  center  distance  of  5.08  cm  (Z  in)  and  a 
horizontal  center  to  center  distance  of  3.91  cm  (1.5  in).  The 
distance  from  the  distributor  to  the  center  of  the  upper  most  tube 
is  Z7.SZ  cm  <10.875  in).  A  5.08  cm  (Z  in)  spacing  separates  the 
front  and  back  ualls  of  the  fluidized  bed  from  the  end  tubes  of 
the  bundel.  The  cross  sectional  area  of  the  bed  is  1.079  sq  m 
(11.51  Sq  Ft).  The  particulate  material  used  in  the  bed  is  a 


Mixture  of  steel  grit  abrasive  having  a  specific  gravity  of  8.1. 


Appendix  £  lists  the  size  distribution  of  the  steel  grit  used 
during  the  sampling  operation.  The  bed  uas  operated  without 
recycling  the  fines  captured  in  the  cyclones.  The  static  bed 
height  of  the  Material  uas  ZZ.S4  cm  (8.875  in)  throughout  the  data 
collection  period. 

Fnui nment  Set -u  a 

Fig.  11  shows  the  placenent  of  the  sanpling  device  inside 
the  fluidized  bed.  The  sanple  trap  uas  positioned  such  that  it 
uas  directly  over  a  spacing  betueen  tubes  1ZZ.85  cm  (3  in)  fron 
the  center  of  the  front-nost  tube!  and  41.31  cm  (16.5  in)  fron  the 
left  uall.  The  height  of  the  sanple  trap  above  the  distributor 
uas  varied  during  the  sanpling  process  as  discussed  in  the  section 
on  sanpling  procedure. 

Fig.  1Z  shows  the  placenent  of  the  bubble  probe  and  the 
anenoneter  probe  with  respect  to  the  sanple  trap.  The  bubble 
probe  uas  placed  directly  below  the  sanple  trap  and  Z6.67  cn  (10.5 
in)  above  the  distributor.  The  probe  extension  uas  placed  at  an 
angle  so  as  not  to  interfere  with  the  sanple  trap  operation.  To 
protect  the  anenoneter  uire  fron  particles  inpacting  it,  a  special 
shield  consisting  of  4  3Z0  nesh  screen  and  an  alunirun  frane  uas 
placed  around  it.  The  anenoneter  probe  uas  attached  to  the  bubble 
probe  extension  with  the  entrance  to  the  anenoneter  probe  Z3.Z1  cn 
(11.5  in)  above  the  distributor.  This  placed  the  entrance  to  the 
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(a)  Right  Side  of  Heat  Exchanger  Tubes 


<b)  Front  of  Heat  Exchanger  Tubes 


Fig.  11  Position  of  Sample  Trap,  Kibble  Probe, 

and  Anemometer  Probe  Above  Heat  Exchanger 
Tube*.  All  (Dimensions  in  cm. 


Sample 

Trap 


anemometer  probe  Z.54  cn  (1.0  In)  auay  from  the  center  of  the 
sanple  trap  and  0.64  cn  (0.25  in)  outside  the  area  directly  below 
the  sanple  trap.  As  a  result,  the  probe  has  a  Mininun  effect  on 
the  air  flow  fron  the  bed  to  the  sanple  trap,  but  will  only 
measure  the  gas  velocity  at  the  edge  of  the  sanple  trap  perimeter. 
The  positioning  of  the  two  probes  above  the  distributor  remained 
constant  throughout  the  sampling  evolution. 

Fig.  13  is  a  block  diagram  showing  the  eqiipnent  used  diring 
the  sampling  operations  and  their,  interconnections  .  Table  1  is  a 
listing  of  the  equipment  used.  The  oscilloscope  time  base  was  set 
for  MANUAL  TRIGGER,  SINGLE  SWEEP  node  and  a  sweep  tine  of  50 
ns/div.  The  channels  of  the  dual  trace  amplifier  were  set  at  5 
volts  per  division  for  the  bubble  probe  and  2  volts  per  division 
for  the  sample  trap  inputs.  On  the  differential  amplifier,  one 
channel  was  not  used  and  the  second  channel  was  set  at  1  volt  per 
division  for  the  anemometer  probe  input. 

I 


6 


Dual  l 
Trace 
Amplifier 


□iff erential 
Amplifier 


Solenoid 

Pouer 

Supply 


Optical 

Detector 

Pouer 

Supply 

Light 

Source 


Bubble 

Probe 


Anemometer 


A 

Anemometer  Probe 


Fig.  13 


Block  Diagram  of  the  Equipment  Used 
During  the  Sampling  Operations. 


TABLE  1 


Oscilloscope  TEKTRONIX  5111  StoraQe  Oscilloscope 

TEKTRONIX  5A18N  Dual  Trace  Anpl 
TEKTRONIX  5AZ1N  Differential  Anpl 
TEKTRONIX  5B18N  Tine  Base 

Anenoneter  Thermal  Systens  Inc. 

10S1-Z  Monitor  and  Power  Supply 
10S4-A  Linearized  Anenoneter 
Anenoneter  uire  u/  #3Z0  screen  guard 

Optical  Signal  Oetector 
Optical  Signal  Detector  Power  Supply 
Optical  Source  and  Power  Supply 
Oscilloscope  Camera 


Listing  of  equipment  used  during  particle  sampling 
operations. 


To  determine  the  fluidization  conditions  uithin  the  bed,  a 


set  of  nanometers  uere  used.  These  nanometers  measured  pressures 
uithin  the  bed,  at  heights  from  4.13  cm  (1.6Z  in)  to  37.15  cm 
(14.SZ  in)  above  the  distributor,  in  Z.54  cm  (1  in)  increments. 

The  pressure  data  corresponding  to  each  trap  position  and  bed 
velocity  is  listed  in  Appendix  I. 

To  determine  the  gas  flow  conditions  uithin  the  bed,  an 
orfice  flou  meter  uith  10  -  1/ZD  taps  uas  located  upstream  of  the 
distributor.  The  computer  program  listed  in  Appendix  G  uas  used 
to  convert  the  pressure  tap  data  to  mean  air  velocities  uithin  the 
bed. 

SgpoUnq  Prpggdupq 

Data  uas  collected  for  four  (4)  mean  bed  velocities  at  six 
(6)  different  sample  heights.  The  sample  trap  uas  placed  at  a 
given  sampling  height  (measured  from  the  distributor  to  the  bottom 
of  the  trap),  and  ten  (10)  samples  uere  collected  at  each  desired 
velocity.  The  trap  position  uas  then  changed  to  a  neu  height. 

During  certain  sampling  conditions,  those  which  involved  low 
sampling  heights  uith  the  higher  air  velocities,  the  paddle  arms 
uould  occasionally  impact  the  sides  of  the  sample  trap  and  not 
close  the  sample  trap  completely.  It  is  assumed  that  this  occured 


when  a  large  bubble  erupted  directly  under  the  paddle  arn  and 


deflected  the  paddle  arn  into  the  side  of  the  trap.  Whenever  this 
occured,  the  trap  uas  de-energized  and  the  closure  cycle  repeated. 

Each  tine  the  height  or  velocity  uas  changed,  a  conplete 
sanpling  cycle  uas  conducted  and  this  sanple  discarded.  This  uas  to 
prevent  any  accumulation  of  particles  (in  the  entrance  to  the 
vacuun  or  purge  lines  on  the  sanple  trap)  fron  being  added  to  the 
first  sanple  at  the  neu  height  or  velocity. 

For  each  set  of  data  at  a  given  height  and  velocity  the 
follouing  information  uas  recorded: 

1)  Fluidized  bed  height  above  the  distributor 

determined  visually  and  by  pressure  measurements. 

Z)  Pressure  upstream  of  orfice  plate  (PI) 

3)  Pressure  difference  across  orfice  plate  (delta  P) 

4)  Air  temperature  in  bed 

5)  Pressure  distribution  in  bed 

6)  Height  of  sanple  trap  above  distributor 

The  follouing  procedure  uas  used  during  sampling: 

For  each  sanple  to  be  taken  uithin  a  data  set. 

A.  Initial  conditions 

1)  Vacuun  air  supply  OFF 

Z)  Sueep  trigger  on  oscilloscope  RESET 


50 


3)  O-scope  memory 


ON-CLEAR 


4)  Solenoid  pouer  supply  suitch  SI  OFF 

5)  Solenoid  trigger  suitch  SZ  ON 

B.  Sampling  Procedure 

1)  Trigger  oscilloscope  sueep  and  wait  until  sweep 
is  at  the  center  of  the  CRT. 

Z)  Close  the  Solenoid  pouer  supply  suitch  SI.  Uhen 
SI  is  shut,  the  oscilloscope  uill  shou  an 
additional  trace.  This  third  trace  is  used  to 
determine  the  closure  tine  relative  to  the 
presence  of  gas  jets  and  bubble  eruptions.  An 
example  of  a  typical  oscilloscope  trace  is  shoun 
in  Fig.  14. 

3)  Open  the  solenoid  trigger  suitch  SZ.  This 
reduces  the  current  to  the  solenoid  to  prevent 
overheating.  The  maximum  alloued  tine  to  let  SZ 
remain  closed  is  five  (5)  seconds. 

4)  Turn  on  the  air  supply  to  the  vacuum  system  and 
leave  on  one  (1)  minute.  The  exhaust  air  from 
the  venturi  on  the  vacuum  system  must  be  directed 
into  the  purge  line  in  an  oscillatory  manner. 

This  uill  agitate  the  particles  within  the  trap 


so  as  to  move  them  into  the  vacuum  line. 


5)  Turn  off  the  vacuum  air  supply. 

5)  Turn  off  the  solenoid  power  supply  switch  SI. 

7)  Remove  sample  from  sample  container  and  place  in 
specimen  bottle. 

8)  Photograph  trace  on  oscilloscope. 

3)  Repeat  from  A. 

Sample  Ueight  Determination 

Each  sample  obtained  was  weighed,  using  a  Torsion  Balance  Co. 
TORBflL  scale,  to  an  accuracy  of  0.01  grams.  The  average  value  and 
standard  deviation  was  then  determined  for  each  set  (specific 
height  and  velocity)  of  samples.  The  ueight  of  particles  in  a 
completely  filled  sample  trap  was  also  determined  for  void 
determination.  These  results  are  listed  in  Appendix  I. 


Determination  of  Particle  Size  Distribution 


To  determine  the  particle  size  distribution  which  occurs  at 
various  heights  for  a  specified  gas  velocity,  three  (3)  samples 
from  each  set  of  data  at  a  specified  height  and  velocity,  and  from 
the  bulk  bed  material  were  analyzed  using  a  Zeiss  Videoplan  Z 
Image  analyzer.  The  software  used  was  the  "Image  Analysis  System 
MOP-Videoplan”  distributed  by  KONTRQN  Electronics  Group. 

For  each  sample  selected,  a  microscope  slide  was  prepared 
using  double  sided  adhesive  tape  on  which  a  portion  of  the 
selected  sample  was  placed.  To  ensure  a  sharp  contrast  was 
achieved,  each  slide  uas  backed  with  white  paper.  The  sample  was 
then  placed  under  a  Zeiss  microscope  to  which  the  image  analyzer 
was  connected  via  a  vidiocon  tube.  The  magnification  used  was 
S0X  which  provided  an  average  view  of  about  8  particles  at  a 
time.  The  analyzer  uas  then  used  to  analyze  the  partices  present 
on  each  slide  of  interest.  The  mode  in  which  the  image  analyzer 
uas  used  uas  the  equivalent  diameter  mode.  This  mode  determined 
the  cross  sectional  area  of  each  particle  sampled  and  calculated 
the  diameter  of  a  circle  with  the  same  area.  Uhen  a  data  set  uas 
completed,  the  data  uas  analyzed  for  particle  size  distribution. 
The  resulting  output  (Appendix  J)  consists  of  a  particle  count  vs 
particle  size  histogram,  a  cumulative  frequency  plot  and  a 
classification  data  list.  A  gaussian  distribution  fit  for  the 
data  is  also  plotted  on  the  histogram  and  cumulative  frequency 


CHAPTER  IV 


COMPUTER  MODEL 


Introdiction 

This  chapter  descibes  the  theory,  logic  flou  and  testing  of  a 
particle  trajectory  conputer  nodel  to  predict  the  particle  loading 
in  the  freeboard  due  to  erupting  buubles.  The  progran,  listed  in 
Appendix  L,  is  written  in  HP  BASIC  2.0  and  uas  run  on  a  Heulett 
Packard  9816  nicroconputer .  The  output  fron  the  nodel  is  discussed 
in  chapter  VI  and  conpared  with  experinental  results.  An  in  depth 
analysis  of  the  progran  logic  and  structure  is  given  in  Appendix  L. 

Mqds.1  Thepcv 

The  nodel  developed  here,  is  based  on  calculating  the 
trajectory  of  a  single  particle  as  it  is  ejected  fron  the  bed 
surface  and  is  acted  upon  by  gravitational  and  drag  forces.  The 
drag  force  is  due  to  the  difference  in  absolute  particle  and  air 
velocities.  The  air  velocity  is  a  conbination  of  the  initial  jet 
velocity  produced  uhen  a  bubble  bursts,  and  the  superficial  bed 
velocity.  To  ensure  that  the  particle  drag  is  calculated 
accurately,  the  follouing  drag  coefficient  correlation  given  by 


Uhite  (171  uas  used. 


Cd  =  24  +  S  +  0.4 

Re  1+  /Re1  (3) 


where: 

Cd  =  Drag  coefficient  for  sphere 
Re  *  Reynolds  number 


Eqn  3  is  valid  over  the  range  0<Re<10^*.  To  calculate  the 
particles  position  and  velocity,  the  computer  uses  a  forward 
difference  method.  Using  Newton's  Law  (ZF=ma)  the  acceleration  of 
the  particle  due  to  gravity  and  drag  is  determined.  Inserting 
this  acceleration  into  Eqn  4,  the  particles  new  velocity  is 
determined. 


V  =  Vo  +  a  t 


(4) 


where: 

V  =  Particles  new  velocity 
Vo  =  Particles  present  velocity 
a  =  Acceleration  of  particle 
t  =  Time  increment  of  calculation 


To  determine  the  particles  new  position,  the  velocity 


calculated  in  Eqn  4  is  inserted  into  Eqn  5. 


(5) 


H  =  Ho  +  (Vo  t  V)  t 
Z 

I 

where: 

H  =  Particles  new  height 
Ho  *  Particles  present  height 
Vo  ~  Particles  present  velocity 
V  =  Particles  new  velocity 
t  =  Tine  increment  of  calculation 


Using  Eqns  3,4, and  5,  the  trajectory  of  the  particle  is  calculated 
fron  the  tine  it  initially  leaves  the  bed  until  the  tine  that  it 
returns  to  the  bed. 

These  calculations  are  repeated  over  a  range  of  particle 
dianeters  fron  80  to  570  nicrons.  By  determining  the  residence 
tine  of  each  particle  within  a  specified  height  increment  (A  H) 
above  the  surface  of  the  bed,  a  particle  density  distribution 
above  the  bed  is  determined.  The  height  increment  (A  H)  used  in 
the  program  is  Z  cm.  At  the  end  of  each  tine  step  when  the  height 
calculation  (Eqn  5)  is  completed,  the  counter  representing  the 
particular  Z  cm  height  increment  which  the  particle  is  in,  is 
incremented  by  one.  Each  particle  size  has  its  own  set  of 
counters  to  allow  individual  particle  analysis. 

The  calculated  density  distribution  is  then  weighted  with  the 


particle  size  distribution  of  the  bulk  bed  material  since  the 


probability  of  a  given  size  particle  being  present  at  a  specified 
height  is  dependent  upon  the  number  of  particles  uithin  the 
system.  This  is  accomplished  by  multiplying  each  height  counter 
of  a  given  particle  size  with  the  number  of  particles  for  that 
given  size  present  in  the  input  bed  distribution.  By  summing  the 
density  values  for  each  set  of  particle  diameters  at  a  given 
height  over  the  entire  freeboard  of  the  bed,  the  overall  particle 
density  above  the  bed  surface  is  determined. 

The  model  assumes  that  all  the  particles  are  ejected 
perpendicular  to  the  surface  of  the  bed  and  are  initially  at  a 
uniform  velocity.  Because  the  model  uses  single  particles  for  the 
analysis,  the  effects  of  multiparticle  interactions  are  not 
included  in  the  model. 

Testing  of  Proof an 

To  evaluate  the  validity  of  the  program,  tuo  tests  uere 
conducted.  The  first  test  compared  the  height  solution  produced 
by  the  computer  with  a  closed  form  solution.  The  second  test 
involved  running  the  program  uith  different  particle  diameter  and 
particle  distribution  height  intervals  to  ensure  that  a  valid 
sample  size  uas  being  used. 

Closed  Form  Solution 


To  determine  a  closed  form  solution  for  particle  height  as  a 


function  of  initial  particle  velocity  and  superficial  velocity  Uo, 
a  force  balance  uas  used.  The  forces  acting  on  a  particle  are 
gravitational  and  drag.  The  gravitational  force,  FI,  is  simply 
the  volume  of  the  particle  njltiplied  by  the  particles  density  and 
the  gravitational  acceleration,  and  can  be  written  as: 


3 

FI  =  -  PoTTO  g 

S  (6) 


where: 

FI  =  Force  due  to  gravity 
Pp  =  Oensity  of  particle 
0  =  Diameter  of  particle 
g  =  gravitational  acceleration 


In  order  to  get  a  closed  form  solution  that  did  not  involve 
non-linear  differential  equations,  Stokes  flou  uas  used  for  the 
closed  form  solution  only.  The  computer  model  used  the  Stokes 
equation  only  to  compare  results  with  the  closed  form  solution, 
afterward,  Eqn  3  uas  used.  Using  the  Stokes  drag  coefficient 
relation,  the  drag  force  on  a  particle  can  be  determined  as: 


FZ  =  3  (Uo  -  Up) ^  D  (7) 

uher  e : 

FZ  =  Force  due  to  drag 
Uo  =  Superficial  bed  velocity 
Up  =  Velocity  of  particle 
/+  -  Absolute  viscosity  of  air 
D  =  Diameter  of  particle 


By  inserting  Eqns  6  and  7  into  Neuton's  Law  (£  F=na)  and 
simplifying,  the  following  differential  equation  is  obtained: 

X  +  Cl  X  =  C2 

where: 

X  =  Acceleration  of  particle 

X  =  Velocity  of  particle 

Cl  =  18.jW  „ 

Pp  0X 

C2  =  Cl  Uo  -  a 

boundary  conditions: 

1 )  t=0  X=Uo 

Z )  t=0  X=0 


This  second  order  linear  differential  equation  can  be 
solved  using  the  given  boundary  conditions  with  the  resulting 
closed  form  solution  given  as: 


(3) 


X  =  CJ  +  ££  t  -  C3  exp<  -Cl  t  > 
Cl 


uhere: 


X  =  Height  of  particle 


Using  the  initial  conditions  listed  in  table  Z,  the  solution 
obtained  using  the  conputer  nodal  (naxinun  height3  1S.Z70  cnf  tine 
to  maximum  height3  0.Z75  sec)  uas  identical  to  three  decinal 
places  uith  the  solution  obtained  using  Eqn  3. 


TABLE  Z 


Superficial  Velocity  (Uo>: 

60.  SB 

cm/s  (Z  ft/s) 

Initial  particle  velocity  (Upo): 

304.8 

cm/s  (10  ft/s) 

Particle  diameter  •' 

Z00 

microns 

Particle  density: 

5000 

kg/n* 

Tine  increment  (computer): 

0.001 

sec 

Jet  velocity: 

0.0 

cm/s 

List  of  paraneters  used  to  check  conputer 
calculations  against  closed  forn  solution. 
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Sample  Size  Sensitivity  Test 

To  ensure  that  appropriate  sample  sizes  were  used  to  minimize 
errors  due  to  coarse  sampling  intervals,  tuo  sensitivity  tests 
uere  run.  One  test  involved  changing  the  particle  diameter 
interval  from  10  microns  to  S  microns.  The  second  test  changed 
the  height  sampling  interval  (A  H)  from  Z  cm  to  5cm  and  then  1  cm. 
For  each  of  the  tests,  the  same  Initial  conditions  uere  input  into 
the  program.  Table  3  shous  the  resulting  output  from  the  program 
listing  the  initial  conditions  and  the  resulting  diameter  versus 
maximum  height  data  calculated.  Fig.  15  shous  a  plot  of  the 
calculated  maximum  height  vs  particle  diameter  data  listed  in 
table  2.-  -  Fig.  16  shous  the  bed  particle  distribution  used  for 
each  of  the  tests. 

Fig.  17a  shous  the  entrainment  calculation  using  a  dimeter 
interval  of  10  microns  and  a  A  H  of  Z  cm.  Fig.  17b  shous  the 
same  calculation  using  a  diameter  interval  of  10  microns  and  a  A  H 
of  5  cm.  The  curve  is  not  as  smooth  but  still  retains  the  same 
general  shape.  The  peak  of  the  curve  shoun  in  Fig.  17b  occurs  at 
a  height  of  about  Z0  cm  uhereas  the  peak  in  Fig.  17a  occurs  at 
about  16  cm.  A  semi-log  plot  of  these  curves  uould  shou  that  the 
slope  of  the  line  to  the  right  of  the  peak  uould  be  larger  for  the 
data  represented  by  Fig.  17b.  The  effect  of  maintaining  A  H  at  5 
cm  but  decreasing  the  diameter  interval  to  5  microns  is  shoun  in 
Fig.  17c.  There  is  no  readily  detectable  difference  betueen  Fig. 


cm/s 

Mean  Bed  Velocity-  57.312 

Initial  Particle  Velocity-  97.2312 
Peak  Jet.  Velocity-  609. S 

Sas  Jet'  Ouration-  .02  s 


Diameter 

Maximum 

Heioht 

"Max  Ht 

T  ime 

Diameter 

Maximum 

Hexoht 

Max  Ht 
T  ime 

urn 

cm 

seconds 

urn 

cm 

seconds 

80 

55.1 

.282 

90 

57.0 

.290 

100 

57.2 

.297 

110 

56.8 

.302 

120 

55.8 

.305 

130 

54.4 

.307 

140 

52.7 

.307 

150 

50.8 

.305 

160 

48.8 

.303 

_  .  170 

46.8 

.299 

180 

44.7 

.295 

190 

42.7 

.290 

200 

40.7 

.285 

210 

38.8 

.280 

220 

37.0 

.275 

230 

35.3 

.270 

240 

33.7 

.264 

250 

32.2 

.259 

.  260 

30.8 

.254 

270 

29.5 

.249 

280 

28.2 

.245 

290 

27.1 

.240 

300 

26.1 

.236 

310 

25.1 

.232 

320 

24.2 

.228 

330 

23.3 

.224 

340 

22.5 

.220 

350 

21.8 

.217 

360 

21.1 

.213 

370 

20.4 

.210 

-  380 

19.8 

.207 

390 

19.2 

.204 

400 

18.7 

.201 

410 

18.2 

.199 

420 

17.7 

.196 

430 

17.3 

.194 

440 

16.9 

.191 

450 

16.5 

.189 

460 

16.1 

.187 

470 

15.7 

.185 

480 

15.4 

.183 

490 

15.1 

.  181 

S00 

14.8 

.179 

510 

14.5 

.178 

520 

14.2 

.176 

530 

14.0 

.  174 

540 

13.7 

.173 

550 

13.5 

.  171 

550 

13.2 

.170 

570 

13.0 

.168 

Table  3  Listing  of  Input  and  Resulting  Maximum 
Particle  Heights  with  Time  to  Maximum 
Height.  Thes  Values  Were  Used  During 
the  Increment  Sensitivity  Tests. 


Particle  Density  T  Relative  Particle  Density 


Fig.  17  Plots  of  Relative  Particle  Density  vs  Freeboard 
Height  shouin  the  effects  of  varying  Diameter 
Interval  and  Height  Interval. 


17c  and  Fig.  17b.  In  Fig.  17d,  the  effect  of  changing  A  H  to  1 
cm  and  diameter  spacing  to  S  microns  is  shown.  The  difference 
between  Fig.  17a  and  Fig.  17d  is  barely  not i cable  and  no 
detectable  change  in  the  slope  to  the  right  of  the  peaks  is 
present.  As  a  result  of  this  analysis,  the  program  was  operated 
with  a  A  H  of  2  cm  and  a  diameter  interval  of  10  microns.  This 
reduced  the  calculation  tine  to  half  of  that  required  uhen  using  a 
diameter  spacing  of  5  microns  when  the  same  total  diameter  span 
was  used  and,  as  was  seen  in  Figs.  17a  and  17d,  the  difference  in 
output  does  not  require  the  finer  increment. 


.  Fig.  ,_18  is  a  plot  of  naan  bed.. velocity  (Uo)  vs  pressure  drop 
through  the  bed  (Pb>.  The  velocity  values  uere  determined  using 
the  program  in  Appendix  6  and  the  pressure  data  listed  in  Appendix 
I.  From  this  plot,  the  minimum  fluidization  velocity  (Umf)  for 
the  bed  conditions  used  during  this  study  is  determined  to  be  0.15 
m/s  (0.5  Ft/sec). 

Entrainment  Analysis 

Table  4  lists  the  averaged  sample,  weights  and  their  standard 
deviations  for  the  samples  (Appendix  I)  collected  by  the  sampling 
apparatus.  Looking  at  the  standard  deviation  of  the  sample 
groups,  the  standard  deviation  is  fairly  large  compared  to  the 
average  values.  However,  visual  observations  of  the  fluidized  bed 
in  operation  would  suggest  that  a  larger  standard  deviation  would 
be  expected.  The  short  sample  cycle  <2  ms)  and  the  bubble  burst 
activity  in  the  bed  are  the  main  reasons  for  this  conclusion.  The 


average  standard  deviation  is  29  X  of  the  average  sample  weight 


and  iha  range  extends  from  68  Z  of  the  average  for  the  sanple  set 
numbers  186  thru  115  to  only  7  X  of  the  average  for  the  samples 
28S  thru  295.  This  variation  is  believed  to  be  due  to  the  nature 
of  the  bubble  activity  uithin  the  bed.  It  uould  appear  that  uhen 
a  bubble  erupts  at  a  time  just  prior  to  closure  of  the  trap,  a 
sample  weight  greater  than  the  average  would  be  obtained. 

However,  when  no  bubble  has  been  present,  the  material  caught  by 
the  trap  should  consist  mainly  of  particles  returning  to  the  bed 
and  therefore  be  below  the  average  weight  caught.  Table  4  also 
shows  the  density  of  each  sanple  collected.  These  values  were 
determined  by  dividing  the  averaged  sanple  weights  by  the  volume 
of  the  sample  trap  which  is  43.82  cubic  cm. 
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TABLE  4 


Superficial  Bed  Velocity  (Uo):  58.1  cn/s  (1.905  ft/s) 

Non- di mens ionalised  velocity  (Uo/Unf):  3.81 


Sample 

Height 

Average 

Standard 

Average 

Numbers 

Above  Bed 

Ueight 

Deviation 

Density 

(cm) 

( grams ) 

(grans) 

(grans/cn  ) 

286-29Si  -  -- 

;  3.8.3*! 

7.22K3  T-Bi' 

7  0.51 

1.68  (-01) 

24G-2S5 • 

7.6  - 

4.48 

0.72 

1.04  (-01) 

11S-1ZS 

12.4 

2.41 

0.76 

5.60  ( -02) 

166-175-5 

17.8 

1.08  — '• 

'  0.16 

2.51  (-02) 

7B-8S 

22.2 

1.32 

0.23 

3.07  (-02) 

17B-18S 

31.1 

0.31 

0.10 

7.21  (-03) 

Superficial 

Bed  Velocity 

(Uo): 

48.3  cn/s 

(1.585  ft/s) 

Non- di mens ionalised  velocity  (Uo/Umf): 

3.17 

Sample 

Height 

Average 

Standard 

Average 

Numbers 

Above  Bed 

Ueight 

Deviation 

Density 

(cm) 

(grams) 

( grans ) 

( grans/cm  ) 

Z76-28S 

5.1 

4.35 

1.05 

1.01  (-01) 

23G-Z45 

8.9 

1.88 

0.ZZ 

4.37  (-02) 

9B-10S 

13.6 

1.36 

0.32 

3.16  (-0Z) 

156-165 

19.1 

0.60 

0.12 

1.39  (-02) 

56-65 

Z4.1 

0.9Z 

0.36 

Z. 14  (-02) 

186-195 

31.8 

0.18 

0.04 

4.18  (-03) 

Superficial 

Bed  Velocity 

(Uo): 

39.5  cm/s 

(1.297  ft/s) 

Non-dimensional ised  velocity  (Uo/Umf): 

Z.S9 

Sample 

Height 

Average 

Standard 

Average 

lumbers 

Above  Bed 

Ueight 

Deviation 

Density 

( cm) 

( grams ) 

( grams ) 

(grams/cm  ) 

266-275 

5.7 

1.31 

0.37 

3.05  (-02) 

226-235 

9.5 

0.52 

0. 17 

1.44  (-02) 

126-135 

14.3 

0.28 

0.08 

6.51  (-03) 

136-145 

19.7 

0.08 

0.04 

1.86  (-03) 

86-95 

24.8 

0.09 

0.05 

2.09  (-03) 

196-205 

32.4 

0.06 

0.02 

1.39  (-03) 
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TABLE  4  (coni) 


Superficial  Bed  Velocity  (Uo):  35.4  cn/s  (1.161  ft/s) 

Non-dimensional ised  velocity  (Uo/Umf):  Z.3Z 


Sample 

Height 

Average 

Standard 

Average 

Numbers 

Above  Bed 

Ueight 

Deviation 

Density 

(cm) 

(grams) 

(grans) 

(grans/cm  ) 

Z56-ZE5 

6.4 

0.63 

0.13 

1.46  (-0Z) 

Z16-ZZS 

10.  Z 

0.3Z 

0.04 

7.44  (-03) 

106-115:7: 

'  14.  r  ^ 

0.2Z.n 

(-¥3 '0.15 

5.11  (-03) 

146-155  ' 

Z0.3 

0.02  ~ 

0.01 

4.65  (-04) 

66-75 

Z6.0 

0.1Z 

0.04 

2.79  (-03) 

Z06-Z1S 

33.0  ~  • 

0.01  - 

0.004 

2. 32  (-04) 

List  of  experimental  data  showing  sample  averages,  standard 
deviations,  heights,  and  velocity  conditions  measured. 
Oensity  values  are  calculated  by  dividing  the  average  sample 
weight  by  the  sample  trap  volume. 


Fig.  19  shows  the  relationship  between  the  average  particle 
density  caught  in  the  sample  trap  and  the  sample  trap  height  above 
the  bed  surface  as  a  function  of  Uo/Umf. 


For  the  relatively  low  fluidization  velocities  used  during 
the  data  measurements,  (maximum  Uo/Umf  s  3.81)  the  freeboard 
height  can  be  assumed  infinite.  Using  this  assumption  implies 
that  all  of  the  particles  return  to  the  bed  and  none  are 
elutriated,  ie  complete  reflux.  Under  these  conditions,  the 
following  equation  has  been  suggested  by  Lewis  et  al  1 31  and  Kumi 
and  Levenspiel  (11  to  model  the  particle  loading  within  the 
freeboard. 
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Particle  Oenaity  <Qn/cn*  ) 


Pr  =  Po  exp<-a  l> 


(10) 


where: 

Pr  =  Particle  density  at  height  1 
Po  =  Particle  density  at  (+)bed  surface 
a  =  Characteristic  particle  decay  length 
1  =  Height  above  bed  surface 


Table  5  shous  the  values  for  the  paraneters  Po  and  a  obtained 


by  linear  regression  for  the  curves  shown  in  Fig.  13. 


Uo/Umf 


TABLE  5 


Po  a  Correlation 

Coefficient 


3.81 

0.Z38 

0.1097 

0.379 

3.17 

0.137 

0.1 0Z3 

0.944 

Z.59 

0.040 

0.1181 

0.937 

Z.3Z 

..0.031 

0.1399 

0.857 

Pr  =  Po  expf-a  1> 

Results  of  linear  regression  analysis  for  particle 
loading  density  (grams/cm  )  vs  height  above  the  bed 
surface  <cm>. 


The  parameter  Po  physically  represents  the  particle  loading 
density  uhich  uould  be  obtained  if  the  sample  uere  taken  at  the 
surface  of  the  bed.  This  is  not  necessarily  the  case  as  is 
indicated  by  the  computer  model  uhich  is  discussed  in  chapter  VI, 
but  is  only  a  parameter  descibing  the  particle  loading 
distribution  in  the  region  of  the  data  obtained.  The  dashed  lines 
belou  freeboard  heights  of  4  cm  indicate  the  region  in  question. 
The  parameter  a  is  a  characteristic  length  of  decay  for  the 
particle  flux. 

A  correlation  between  the  values  for  Po  in  table  5  and  Uo/Umf 
is  shown  in  Fig.  Z0.  This  plot  shous  that  Po  is  closely  related 


with  Uo/Umf.  Po  varies  with  (Uo/Umf  -  1)  approximately  to  the  Z.9 


Po  (grama/cm^  ) 


power,  where  as  Ualsh  et  al  (10]  determined  the  coefficient  to  be 
about  2.1.  Uen  and  Chen  £51  reported  a  correlation  for  particle  flux 
which  is  proportional  to  bubble  diameter  and  ( Uo/Unf  -  1 )  to  the 
S/Z  power.  These  correlations  are  listed  in  table  S. 

TABLE  6 

H  -=j 3.- Presenter-work  1 ;  B  »  0 . 1 3  (Uo/Unf  -  1)  '  kg/  m  sec 

Uen  and  Chen  B  =  1.34  (04)  (Uo-Umf)  kg/  n  sec 
Ualsh  et  al  B  =  18  (Uo  -  Umf)  kg/  n  sec 

Comparison  of  least  square  fit  relations  for  Po 
as  functions  of  Uo  and  Umf. 

The  differences  in  these  correlations  are  due  to  the 
different  bed  configurations  in  which  the  data  was  taken  and  the 
measurement  technique  used.  The  present  uork  utilized  a  bed  with 
a  relatively  closely  spaced  tube  configuration  and  steel  grit 
(S.G.  8.1,  median  size  Z30  microns)  for  the  bed  mass.  Both 

ascending  and  descending  particle  fluxes  were  captured  in  the 
sample.  The  uork  of  Ualsh  et  al  (101,  used  a  bed  with  two  (Z) 
widely  spaced  horizontal  serpentine  tubes  and  Ottawa  sand  (S.G. 

Z.6,  median  size  755  microns).  Also,  only  descending  particle 
flux  was  used  in  determining  their  relations.  The  correlations  of 
Uen  and  Chen  153  are  a  result  of  studies  conducted  on  previous 
research  using  cylindrical  column  beds  and  low  mass  bed  materials 


(S.S.  0.8  -  Z.G).  The  data  for  these  analysis  is  based  mainly  on 

pressure  measurements. 

It  has  been  observed  by  Leuis  et  al  C31  and  Uen  and  Chen  (51 
that  a  is  not  a  strong  function  of  Uo.  Both  of  these  studies 
recommend  that  the  char act eristic  particle  decay  length,  1/a, 
could  be  approximated  by  an  expression  of  the  form: 


1/a  =  C  Uo 


(11) 


Table  7  is  a  list  of  correlations  obtained  by  other  studies 
and  in  the  present  uork.  The  study  by  Leuis  et  al  C 31  uas 
conducted  uith  75  micron  glass  spheres  in  a  cylindrical  bed.  Fig. 
21  shows  the  relationship  between  1/a  and  Uo  in  the  present  uork. 


TABLE  7 


Present 

uork 

1/a  = 

(0.13 

+ 

0.03 

s  )Uo 

n 

Leuis  et 

al 

i/a  = 

(1.4Z 

+ 

0.14 

s  )Uo 

n 

Uen  and 

Chen 

1/a  = 

(0.Z5 

+ 

0.09 

s  >Uo 

m 

Ualsh  et 

al 

1/a  = 

(0.3Z 

+ 

0.05 

s  )Uo 

n 

Comparison  of  least  square  fit  relations  for 
1/a  as  functions  of  Uo  ( m/ s > . 
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Two  methods  for  determining  particle  size  distribution  uere 
used.  The  sieve  method,  described  in  Appendix  E,  uas  used  for  bed 
material  analysis  only,  due  to  the  small  sample  sizes  obtained  from 
the  trap.  To  analyze  the  small  trap  samples,  an  image  analyzer 
uas  utilized.  The  procedure  used  uith  the  image  analyzer  is 
described  in  the  experimental  procedure  chapter. 

The  sieve  data  representing  bed  mass  as  a  function  of 
particle  size,  is  listed  in  Appendix  E.  Fig.  ZZa  shous  the  bed 
mass  data  converted  to  particle  number  as  a  function  of  particle 
diameter  uhere  the  interval  of  particle  diameter  betueen 
successive  measurements  is  10  microns.  These  plots  are  faired 
from  the  data  listed  in  Appendicies  E  and  J.  Fig.  ZZb  shous  the 
bed  mass  distribution.  The  conversion  from  a  mass  distribution  to 
a  particle  number  distribution  uas  calculated  assuming  that  each 
particle  uas  spherical  in  shape.  The  volume  corresponding  to  a 
given  particle  diameter  uas  multiplied  by  the  particle  density  to 
get  a  unit  particle  mass.  The  mass  fraction  of  the  sieve  analysis 
corresponding  to  the  specified  particle  diameter  uas  then  divided 
by  the  unit  particle  mass  to  obtain  the  representitive  particle 
number.  The  overall  resulting  particle  number  distribution  curve 
uas  then  normalized  uith  respect  to  a  maximum  value  of  130. 

The  particle  size  distribution  data  obtained  from  the  image 


analyzer  is  listed  in  Appendix  J.  The  inage  analyzer  uas  used  to 
analyze  the  bed  nass  and  sample  data  from  the  Uo/Umf  =  3.81  data 
set.  The  bar  graphs  and  the  cumulative  percentage  plots  shou  the 
number  of  particles  vieued  by  the  image  analyzer  plotted  as  a 
function  of  normalized  particle  diameter. 

Fig.  ZZa  shows  a  plot  of  the  particle  number  distribution  as 
a  function  of  particle  diameter  as  determined  by  the  image 
analyzer.  The  data  is  normalized  uith  respect  to  100  and  is 
compared  uith  the  data  as  determined  by  the  sieve  analysis.  The 
descrepancy  betueen  the  tuo  plots  can  be  explained  by  the  methods 
used  to  determine  the  respective  data.  The  analyzer  first 
determines  the  cross  sectional  area  of  the  vieued  particle.  A 
circle,  shaving,  the  same. cross  sectional  area,  is  then  calculated. 
This  results. in  averaging  the  smaller  minimum  diameter  uith  the 
larger  maximum  diameter  of  all  particles  uith  long  cylindrical  or 
ellipsoidal  shapes.  In  the  sieve  houever,  a  large  portion  of  the 
particles  uill  pass  through  the  sieve  screen  by  means  of  the  small 
cross  sectional  area  presented  by  their  longitudinal  direction. 

To  be  consistant  in  the  following  sections,  the  image  analyzer 
data  uill  be  used  to  correlate  all  particle  density  distributions. 

The  particle  size  distribution  as  a  function  of  height  above 
the  bed  surface  uas  evaluated  using  the  data  collected  at  Uo/Umf  ~ 
3.81.  The  image  analyser  uas  used  to  determine  the  distribution 
by  viewing  three  (3)  random  samples  from  each  of  the  six  (S) 
height  positions.  Table  8  lists  the  various  statistical  values 


for  the  particle  size  distribution  data  obtained  (Appendix  J). 


TABLE  8 


Freeboard 

Average 

Median 

Mode 

Diameter  at 

Height 

Diameter 

Diameter 

Diameter 

S0Z  of 
Cumulative 
Number 

Distribution 

(cm) 

( microns ) 

<  microns ) 

( microns ) 

( microns ) 

3.8 

199 

192 

142 

159 

7.S 

219 

214 

187 

203 

12.3 

204 

189 

162 

167 

17.8 

179 

164 

1S4 

179 

22.2 

197 

179 

152 

177 

~  31.1 

212 

197 

194 

190 

Staistical  values  for  particle  number  distribution 
as  a  function  of  freeboard  height.  A  complete 
listing  of  the  data  is  given  in  Appendix  J. 


Fig.  23  shous  these  values  plotted  against  the  bed  height  at 
uhich  they  uere  taken.  The  excursion  of  points  at  the  7.S  cm 
height  is  assumed  to  be  due  to  analysis  error.  The  linear 
regression  lines  for  the  average  and  median  values  shou  that  they 
are  ueak  functions  of  collection  height.  For  the  particle 
diameter  values  representing  the  mode  and  the  50Z  point  on  the 
cumulative  percentage  plot,  the  linear  regression  lines  shou  a 


stronger  dependence  on  collection  height. 


Freeboard  HelQht  (cm) 

Analysis  of  Variation  of  Particle 
Size  vo  Freeboard  HelQht. 


rigs.  24  thru  27  shou  examples  of  the  oscilloscope  traces 
given  in  Appendix  K.  Figs.  24  and  25  represent  the  traces 
obtained  during  lou  Uo  coditions.  Under  these  conditions,  the 
normal  bubble  probe  output  is  high  (+5  volts  dc).  When  a  bubble 
errupts,  particles  are  throun  from  the  surface  of  the  bed  and 
eclipse  the  light  path  at  the  tip  of  the  probe  (trace  goes  to 
zero).  Figs.  26  and  27  represent  the  traces  obtained  at  higher 
Uo  conditions  when  the  bubble  probe  is  normally  eclipsed  (bubble 

"  r 

probe  output  is  lou,  zero)  by  particles.  The  presence  of  a  bubble 

is  detected  by  a  high  outputfrom  the  bubble  probe  due  to  the 

bubble  creating  a  void  through  uhich  the  light  beam  can  pass.  The 

anemometer  output  just  above  the  bubble  probe  is  used  to  determine 

the  velocity  of  the  gas  jet  leaving  the  erupting  bubble.  As  seen 

in  the  figures,  there  is  not  always  a  bubble  associated  uith  a  gas 

jet  and  vice  versa.-  This  is  due  to  the  anemometer  probe  being 

slightly  offset  from  the  bubble  probe  as  described  in  the 

experimental  procedure  section. 

The  gas  jet  in  Fig.  24  is  delayed  5  ms  and  the  gas  jet  in 
Fig.  25  is  delayed  35  ms  from  the  point  uhere  the  bubble  traces 
begin  their  excursion  to  the  zero  ~(  eel ipsed)  condition.  Since  the 
seperation  between  the  bubble  probe  and  the  anemometer  is  1  cm, 
this  results  in  an  estimated  jet  velocity  of  200  cm/s  (Fig.  24) 


Fig.  Z6  Oscilloscope  Trace  of  Bubble  Probe, 
Anenoneter  Probe,  and  Solenoid 
Actuation  at  Higher  Uo.  _____ 


Fig.  Z7 


Oscilloscope  Trace  of  Bubble  Probe, 
Anenoneter  Probe,  and  Solenoid 
Actuation  at  Higher  Uo. 


and  28  cm/s  (Fig.  25).  Using  the  calibration  curve  in  Appendix 
M,  the  velocities  are  determined  to  be  518  cn/s  for  the  jet  in 
Fig.  24  and  782  cn/s  for  the  jet  in  Fig.  25.  The  large 
difference  in  these  values  could  be  due  to  a  bulge  in  the  bed 
surface  from  the  bubble  underneath  eclipsing  the  bubble  probe. 

If  this  is  the  case,  there  uauld  be  a  delay  betueen  the  tine  at 
uhich  the  bubble  probe  uas  eclipsed  and  the  eruption  of  the  bubble 
fron  the  bed  surface. 


An  analysis  of  Figs.  25  and  27  results  in  the  sane 
discrepancies  betueen  velocities  calculated  fron  delay  tines  and 
measured  by  the  anemometer  probe.  Fig.  27  shows  another  phenomena 
uhich  occurs  quite  frequently.  The  gas  jet  appears  before  the 
bubble  Is  detected. Therefore: .several  possibilities  uhich  can 
explain  these  occurances.  First,  the  bubble  occurs  off  center 
from  the  bubble  probe.  Under  this  condition,  the  bubble  nay  erupt 
and  initiate  a  jet  uhich  is  registered  by  the  anemometer.  The 
bubble  then  continues  to  rise  and  the  side  of  the  bubble  is 
registered  by  the  bubble  probe.  This  explanation  can  be  altered 
to  include  bubbles  coalescing  belou  the  surface.  Under  this 
condition,  a  bubble  may  be  pulled  into  the  vortex  of  an  already 
erupting  bubble  and  it  eclipses  the  bubble  probe. 


An  analysis  of  the  jet  velocities  and  durations  for  all  the 
oscilloscope  traces  in  Appendix  K,  resulted  in  the  data  listed  in 


table  9.  In  later  calculations,  the  values  to  be  used  for  the  jet 


velocity  and  duration  uill  be  609.6  cm/s  (20  ft/s)  and  20  ms 


respectively. 


TABLE  9 


Jet  Velocity 

Jet  duration 


Average 

G89  cm/ s 
Z2.S  ft/3 

21.25  ms 


Stnd  Qev 

533  cm/ s 
18.  ft/s 

13.22  ms 


Average  and  standard  deviation  of  jet  velocity 
determined  from  oscilloscope  traces  in 
Appendix  K. 


Sample  Ueiqht  Versus  Bed  Activity  Correlation 

The  main  purpose  of  inserting  the  bubble  and  anemometer  probe 
beneath- the- sampling' apparatus  uas  to  determine  uhether  a 
correlation  exists  betueen  the  sample  ueight  collected  and  the 
presence  of  bubble  eruptions  and  gas  jets.  To  evaluate  the 
photographs  taken  of  the  oscilloscope  traces  (Appendix  K),  the 
following  information  uas  required: 

A.  The  average  velocity  of  the  particles  as  they 
travel  from  the  bed  surface  to  the  trap. 

B.  The  distance  the  particle  must  travel  to  reach 
the  trap. 

C.  The  closure  time  of  the  trap  relative  to  the 
bubble  eruptions  and  gas  jets. 
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The  average  particle  velocity  uas  calculated  by  determining 
the  average  particle  size  and  using  the  height  output  from  the 
computer  trajectory  model.  The  average  particle  size  as 
determined  by  the  image  analyzer  uas  180  micons.  To  determine  the 
average  velocity  using  the  model,  the  following  initial  conditions 
uere  input  to  the  program: 

Superficial  velocity  (Uo)  =  57.9  cm/s  (1.9  ft/s) 

This  uas  the  actual  velocity 
measured  during  sampling. 

Jet  velocity  (Uj)  =  609.5  cm/s  (Z0  ft/s) 

This  value  uas  determined  from 
the  anemometer  data. 

Jet  duration  (t  jet)  *  Z0  ms 

This  value  uas  determined  from 
the  oscilloscope  traces. 

The  height  attained  by  a  180  micron  particle,  as  determined 
by  the  model  uith  the  above  conditions  is  44.7  cm  in  0.Z95 
seconds.  This  results  in  an  average  velocity  of  44.7/0. Z95  = 

151.5  cm/s  (4.97  ft/s). 

The  distance  a  particle  must  travel  to  reach  the  center  of 
the  trap  from  the  bed  surface  is  obtained  from  the  data  listed  in 
Appendix  I.  Since  trap  height  is  measured  from  the  bed  surface  to 
the  bottom  of  the  trap,  4  cm  (0.13  ft)  must  be  added  to  the  trap 


heights  to  obtain  the  distance  to  the  center  of  the  trap.  The 
total  tine  required  for  the  particles  to  leave  the  bed  surface  and 
arrive  at  the  center  of  the  trap  can  nou  be  determined.  Table  10 
shous  the  transit  tine  for  each  set  of  data  for  which  the 
oscilloscope  traces  were  photographed  (Appendix  K). 


TABLE  10 


Data 

Trap 

Trap 

Total 

Time 

set 

Height 

Height 

Transit 

Before 

Number 

Bottom 

Center 

Time 

Actuation 

( cn) 

( cn) 

( ns ) 

(ms) 

S6-GS 

24.13 

27.13 

179 

136 

SS-7S 

26.03 

30.03 

198 

155 

76- 8S 

22.23 

26.23 

173 

130 

86-9S 

24.76 

28.76 

190 

147 

96-105 

13.65 

17.65 

116 

73 

List  of 

transit  tines  for 

particles  traveling  from 

the  bed 

surface  to 

the  center  of  the  trap 

.  The  tine 

prior  to  actuation 

of  the 

sample  trap  is 

also  shown 

( Total 

tine  -  42.6 

ns) . 

The  total  closure  tine  of  the  sanpling  apparatus  (closure  of 
solenoid  power  supply  switch  SI  until  trap  is  shut)  is  deternined 
in  Appendix  D  to  be  42. B  ns.  This  value  is  subtracted  from  the 
particle  transit  tine  to  obtain  the  tine  before  actuation  in  which 
particles  leaving  the  bed  surface  will  be  caught  in  the  trap. 


To  deternine  whether  or  not  a  correlation  exists  between  the 


ueight  of  a  sample  and  the  amount  of  bed  activity  present  in  the 
bed,  a  weighted  analysis  uas  used.  Each  sample  was  evaluated 
three  ways.  First,  if  the  sample  ueight  uas  less  than  the  average 
sample  ueight,  a  ueight  factor  <U>  of  -1  uas  assigned  to  it.  If 
the  sample  ueight  uas  greater  than  the  average,  a  ueight  factor  of 
+1  uas  assigned.  Second,  the  oscilloscope  trace  uas  analyzed  at 
the  point  corresponding  to  the  “tine  before  actuation"  listed  in 
Table  10.  This  position  is  the  tine  before  the  solenoid  power 
supply  trace  is  present  (Fig.  Z8>.  If  the  trace  shoued  signs  of 
a  bubble  eruption  or  a  gas  jet  at  this  point,  the  position  factor 
(p)  uas  assigned  a  value  of  +1.  If  no  activity  uas  present,  the 
position  factor  uas  assigned  a  value  of  -1.  Third,  each  trace  uas 
analyzed  again  at  the  “before  actuation  tine”  but  bed  activity 
within  a  +/-  25  ns  region  uas  counted.  If  there  uas  bed  activity 
uithin  this  region,  the  region  factor  (r)  uas  assigned  the  value 
of  +1.  If  no  activity  uas  present,  the  region  factor  uas  assigned 
the  value  of  -1.  The  following  equation  uas  used  to  assign  a 
correlation  factor  to  each  sample : 

Q  =  U  p  +  U  r  (10) 

where: 

Q  -  Correlation  factor 

U  =  +1  if  sample  is  >  average  of  data  set 

-1  if  sample  is  (  average  of  data  set 

p  =  -1  if  no  activity  is  present  in  bed 

+1  if  bubble  or  gas  jet  activity  present 


r  =  -1  if  no  activity  is  present  +/-  25  ns 
+1  if  activity  is  present  +/-  25  ns 


Table  11  shows  the  correlation  factor  for  each  sanple.  The 
resulting  average  correlation  factor  for  all  samples  is  0.32. 

This  value  suggests  that  no  correlation  can  be  made  between  the 
sanple  weight  and  bed  activity.  There  are  several  explanations 
for  this.  Bubbles  not  detected  by  the  bubble  probe  or  debris  from 
previous  bubbles  returning  to  the  bed  surface  can  influence  the 
sanple  weight  by  increasing  the  amount  of  particles  collected. 
Also,  bubbles  of  smaller  size  have  a  smaller  velocity  and  would 
therefore  produce  a  particle  stream  uhich  is  either  slower  or  just 
delayed  in  leaving  the  bed  surface.  This  would  result  in  lower 
sanple  weights  than  expected. 


TABLE  11 


Sample 

Mjmber 

U 

P 

r 

56 

1 

1 

1 

57 

-1 

1 

1 

SB 

1 

1 

1 

53 

-1 

-1 

-1 

50 

1 

-1 

-1 

52 

-1 

1 

1 

S3 

1 

1 

1 

64 

-1 

-1 

-1 

65 

-1 

-1 

1 

67 

1 

1 

1 

69 

-1 

1 

1 

70 

-1 

1 

1 

77 

1 

1 

1 

79 

-1 

1 

1 

S3 

1 

1 

1 

84 

1 

1 

1 

86 

-1 

-1 

1 

88 

1 

1 

1 

89 

-1 

1 

1 

90 

-1 

-1 

1 

9Z 

-1 

-1 

1 

93 

1 

1 

1 

94 

-1 

-1 

1 

95 

1 

1 

1 

96 

-1 

-1 

1 

97 

-1 

-1 

-1 

98 

1 

-1 

1 

99 

1 

-1 

-1 

100 

1 

-1 

-1 

101 

1 

1 

1 

102 

-1 

-1 

1 

Average  - 


List  of  samples  and  their  correlation  parameters, 
resulting  average  value  for  Q  indicates  that  no 
correlation  can  be  made  from  the  data  obtained  to 
indicate  by  sample  weight  whether  or  not  any  bed 
activity  occured  below  the  sample  trap. 
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The  input  to  the  model  consists  of  the  following  five  (5) 
parameters: 

1)  Superficial  bed  velocity  (Uo) 

Z)  Initial  particle  velocity  (Upo) 

3)  Peak  gas  jet  velocity  (Uj) 

4)  6as  jet  duration  (tj) 

5)  Particle  distribution  of  the  bed  mass 

---  The  baseline  values  for  each  of  these  inputs  uas  determined 
to  be  as  close  to  the  actual  experimental  bed  conditions  as 
possible. 

Superficial  Bed  Velocity 

The  superficial  bed  velocity  uas  determined  directly  from  the 
experimental  data.  For  the  bed  conditions  discussed  in  this 
section,  the  superficial  velocity  used  is  the  same  as  the  highest 
velocity  condition  under  which  the  particle  sampler  uas  used. 

The  superficial  velocity  uithin  the  bed  for  these  samples  uas 


calculated  to  be  S7.9  cn/s  (1.9  ft/s>. 

Initial  Particle  Velocity 

To  determine  the  initial  velocity  of  the  particle,  it  was 
assumed  that  the  particle  uas  located  at  the  nose  of  a  bubble  and 
uould  therefore  have  the  bubbles  velocity.  To  determine  the 
bubble  velocity,  Kunii  and  Levenspiel  III  give  the  follouing 
equation: 


Ub  =  0.711 


'4 

Ig  DbJ  ♦  Uo  - 


Umf 


where: 


(11) 


Ub  =  Bubble  velocity 

Uo  =  Superficial  bed  velocity 

:  -----  -• - - --  Umf  =  Mini  mm  fluidization  velocity 

..  .  _  Ob  =  Bubble  diameter 

•  — rg~~  graviatational  acceleration 


Observations  of  the  bed  material  during  the  sampling 
operation  suggests  that  the  average  bubble  diameter  present  in  the 
bed  is  approximately  S  cm  (Z.4  in).  For  Uo  =  57.9  cm/s,  Umf  = 

1S.Z  cm/s  and  Ob  -  5cm,  the  bubble  velocity  is  calculated  to  be 
97. Z  cm/s  (3.19  ft/s).  This  value  uas  used  as  the  initial 


particle  velocity  for  the  base  line  data. 


Peak  Gas  Jet  Velocity 


To  deternine  an  average  peak  jet  velocity,  the  anenoneter 
output  on  the  oscilloscope  traces  (Appendix  K>  were  analyzed.  To 
deternine  the  relationship  betueen  the  voltage  output  fron  the 
anenoneter  (uhich  is  displayed  on  the  oscilloscope)  and  the  actual 
gas  jet  velocity,  the  anenoneter  uas  calibrated  in  a  uind  tunnel. 
The  resulting  calibration  curve  for  the  anenoneter  output,  as 
displayed  on  the  oscilloscope  trace,  is  given  in  Appendix  M.  The 
average  peak  velocity  of  the  anenoneter  traces  analyzed  in  chapter 
V  uas  deternined  to  be  approxinately  609.6  cn/s  (20  ft/s). 

Gas  Jet  OUration 

' Tode  tor-mine  an  average  gas -.jet  duration  tine  to  use  as  a 
baseline,  the  anenoneter  traces  in  Appendix  K  uere  used.  The  tine 
duration  of  each  gas  jet  analyzed  is  deternined  directly  fron  the 
oscilloscope  trace.  The  average  duration  of  a  gas  jet  uas 
deternined  to  be  approxinately  20  ns  in  chapter  V. 

Particle  Distribution  of  the  Bed  Mass 

Tuo  particle  distributions  for  the  bed  nass  uere  available 
for  use  as  the  baseline  data.  Both  particle  distributions  uere 
deternined  using  the  sane  sanple  but  analyzed  using  different 
analysis  techniques.  The  tuo  methods  used  are  the  sieve  method 
and  the  image  analyzer  method,  both  of  uhich  are  described  earlier 


in  chapter  V.  The  tuo  distributions  are  shown  in  Fig.  29. 

The  bed  distribution  deternined  by  the  inage  analyzer  was  used  to 
be  consistent  with  the  particle  distributions  neasured  at  various 
freeboard  heights. 


Fig.  30  shows  the  resulting  naxinun  particle  height  obtained 
as  a  function  of  particle  dimeter  for  the  baseline  conditions. 

The  general  shape  of  the  curve  is  deternined  by  the  doninating 
force  acting  on  the  particle.  For  large  particles,  the  doninating 
force  at  these  air  velocities  is  gravity.  Therefore,  nonentun 
(initial  particle  velocity)  is  the  controlling  factor  deterninlng 
the  naxinun  height  attained  by  the  particle.  As  the  particle  size 
decreases,  the  proportion  of  drag  force  to  gravitational  force 
becones  larger.  For  the  snaller  particles,  the  drag  force  becones 
the  doninate  force.  The  result  is  that  the  naxinun  height 
attained  by  a  particle  continues  to  increase  as  particle  size 
decreases.  This  continues  until  the  value  of  Uo  approaches  the 
terninal  velocity  of  the  snaller  particles.  At  this  point,  the 
particle  is  totally  doninated  by  the  air  flow  in  the  bed.  This  is 
seen  in  the  snail  decrease  in  height  attained  by  particles  less 
than  110  nicrons  as  a  result  of  the  particles  rapid  deacceleration 
to  Uo  after  the  gas  jet  has  stopped  uhile  the  larger  particles 
continue  a  little  higher  due  to  nonentun.  The  particle  height 
continues  to  decrease  for  a  short  tine  until  Uo  becones  equal  to 


or  greater  than  the  terminal  velocity  for  the  remaining  particles. 
These  particles  uould  be  elutriated  from  the  system. 

Figs.  31  thru  36  show  the  individual  particle  densities  at 
different  heights  above  the  bed.  All  of  the  heights  shou  an 
increase  in  particle  density  for  the  smaller  diameter  particles. 
This  is  the  result  of  the  smaller  particles  falling  at  their 
respective  terminal  velocities  uhich  is  slouer  than  for  the  larger 
particles.  The  result  is  that  the  smaller  the  particle  velocity, 
the  longer  the  particle  exists  uithin  a  given  height  region  and 
thus,  the  larger  the  particle  density.  It  is  also  due  to  the 
larger  number  of  smaller  particles  present  uithin  the  system. 

Fig.  37  shows  the  particle  density  distribution  above  the 
bed..  The  tuo  curves  represent  the  effect  of  varying  the  range  of 
particle. diameters  used  in  the  model.  The  two  ranges  are  80  -  570 
microns  and  80  -  1070  microns ^  The  peak  of  the  80  -  570  micron 
distribution  is  at  about  18  cm  of  bed  height  uhereas  the  other 
distribution  peaks  at  about  11  cm.  The  region  shown  on  the  graph 
below  these  heights  is  the  area  known  as  the  splash  zone.  The 
shape  of  this  part  of  the  curve  is  due  to  the  initial  acceleration 
of  the  particles  leaving  the  surface  of  the  bed  by  the  gas  jet. 

The  particles  then  begin  to  slow  down  due  to  drag  and  gravity 
resulting  in  an  increase  in  particle  density.  The  peak  on  this 
curve  coincides  with  the  maximum  height  attained  by  the  largest 
particles,  and  therefore,  it  can  be  assumed  that  the  change  in 
slope  is  due  to  the  loss  of  particles  as  they  return  to  the  bed. 
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Fig.  34  Relative  Particle  Number  vs  Particle 
Diameter  for  Baseline  Conditions. 
Freeboard  Height  of  18  cm. 
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Relative  Particle  Number  vs  Particle 
Oiameter  for  Baseline  Conditions. 
Freeboard  Height  of  31  cm. 
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Fig.  38  In  Particle  Density/Unit  Volume 
vs  Freeboard  height  for  Baseline 
Conditions. 


Fig.  38  is  a  seni-ln  plot  of  the  sane  data  as  plotted  in  Fig. 
37.  The  sharp  drop  at  the  right  end  of  the  plot  is  due  to 


analyzing  particles  uith  a  nininun  dianeter  of  88  microns.  Had 
smaller  particles  been  included,  the  slope  uould  have  approached 
zero  instead  of  infinity  uhich  uould  model  the  elutriation  of 
particles.  The  slope  of  the  line  to  the  right  of  the  peaks  shown 
in  Fig.  38,  uhich  are  located  at  11  and  18  cm  of  height  above  the 
bed  surface  are  the  sane  and  are  approximately  -0.117  grans/cn  . 

To  decrease  computation  tine,  the  80  -  570  micron  particle 
distribution  uas  used  in  the  following  sensitivity  analysis  as  the 


The  parameters  changed  and  the  values  used  are  as  follows- 


Superficial  velocity  (Uo>: 

Baseline  Uo:  57.3  cm/s  (1.3  ft/s) 

1.  30.5  cm/s  (1  ft/s) 

2.  31.4  cm/s  (3  ft/s) 

Initial  particle  velocity  (Upo>‘- 

Baseline  Upo^  11S.Z  cm/ s  (3.13  ft/s) 

1.  61.0  cm/s  (2  ft/s) 

2.  152.4  cm/s  (5  ft/s) 


3.  305.0  cn/s  ( 10  ft/s) 

Sas  Jet  velocity  (Uj): 

Baseline  Uj :  B03.S  cn/s  (20  ft/s) 

1.  305.0  cn/s  <10  ft/s) 

2.  457.2  cn/s  (15  ft/s) 

6as  jet  duration  (tj): 

Baseline  tj:  20  ns 

1 .  1 0  ns 

2.  30  ns 

Particle  distribution  of  bed  nass: 

.  Baseline  distribution:  Inage  Analyzer 

,  .  . .  -  .  1.  Sieve  analysis  data 

The  bed  characteristics  uhich  are  evaluated  against  the 
baseline  characteristics  are: 

1.  The  naximun  height  attained  by  each  particle  fron  80  to 
570  ni crons. 

Z.  The  density  of  each  particle  size  (80  -  570  nicrons)  at 
freeboard  heights  of: 

a.  4  cn  ( 1 .5  in) 


b.  8  cn  ( 3. 1  in ) 


c.  12  cn  (4.7  in) 


d.  18  cn  (7.1  in) 

e.  22  cn  (8.7  in) 

f.  31  cn  (12.2  in) 

3.  The  particle  density  distribution  in  the  freeboard. 


Variation  of  Superficial  Velocity  (Uo) 

Fig.  39  shows  the  influence  of  Uo  on  the  naxinun  height 
attained  by  the  particles.  As  uas  described  in  the  previous 
section,  the  larger  particles  are  doninated  by  nonentun  and  not 
drag.  -This  is~  shown  inFig.  ;  39  by  the  snail  change  in  naxinun 
height  attained  by  the  large  particles 'due  to  a  change  in  Uo.  As 
the  particle  size  decreases,  the  drag  force  becones  the  doninate 
force  and  the  effect  of  Uo  on  particle  height  increases.  As  seen 
in  Fig.  39,  a  change  in  Uo  produces  a  noderate  change  in  the 
naxinun  height  attained  by  the  snail  particles. 

Figs.  40  thru  45  shou  the  effect  of  Uo  on  individual 
particle  densities  at  different  heights  above  the  bed.  All  of  the 
heights  shou  an  increase  in  particle  density  for  the  snaller 
dianeter  particles  as  Uo  is  increased.  This  is  the  result  of  the 
snaller  particles  falling  at  their  respective  relative  terminal 
velocities  which  is  slower  for  snaller  particles.  A  decrease  in 
the  naxinun  particle  dianeter  present  at  a  given  freeboard  height 
as  Uo  is  decreased  can  be  observed  in  Figs.  40  thru  45.  This 
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decrease  is  due  to  the  decrease  in  lift  given  to  the  ascending 


particle  as  Uo  is  decreased.  As  a  result,  all  particles  achieve  a 
lower  maxima  height  when  Uo  is  decreased. 


Fig.  48  shows  the  effect  of  Uo  on  the  particle  density 
distribution  above  the  bed.  There  is  a  slight  increase  in  overall 
particle  density  with  an  increase  in  Uo.  As  expected,  there  is 
also  an  increase  in  the  naximn  height  attained  by  the  particles 
with  an  increase  in  Uo.  .  Fig.  47  is  a  seni-ln  plot  of  the  sane 
data  plotted  in  Fig.  48.  Table  12  lists  the  slopes  of  the  lines 


to  the  right  of  the  peak  shown  in  Fig.  47  which  is  located  at 
about  18  cn  of  height  above  the  bed  surface.  This  table  shows  an 
increase  in  the  particle  distribution  slope  as  Uo  is  decreased. 


TABLE  1Z 

Uo 

cn/s 

Slope 

gms/cn 

Intercept 

gns/cn 

30.5 

57.8 

91.4 

•0.136 

-0.117 

-0.098 

7.036 

6.758 

6.473 

Effect  of  Uo  on  the  slope  of  the 
particle  density  distribution  as  a 
function  of  height  for  the 
distributions  shown  in  Fig.  47. 
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Variation  of  Initial  Particle  Velocity  (Upo) 

Fig.  48  shous  the  relation  between  initial  particle  velocity 
(Upo)  and  the  nxiwn  height  attained  by  particles  of  different 
size.  This  figure  also  shous  that  the  larger  particles  are  nainly 
momentum  dependent  uhile  the  smaller  particles  are  drag  dependent. 
If  Upo  were  increased  further,  the  larger  particles  would  continue 
to  increase  their  maximum  height.  The  snaller  particles  would 
approach  a  height  linit  which  is  dependent  upon  the  superficial 
velocity  in  the  bed. 

Figa.~  49  thru  54  show  the  effects  of  increasing  Upo  on  the 
individual  particle  densities  for  increasing  heights.  In  general, 
the  individual  particle  distributions  undergo  the  sane  relative 
changes  from  the  bed  surface  to  the  maxinim  height  position.  The 
difference  being  the  height  above  the  bed  surface  at  which  the 
particular  distribution  is  present. 

Fig.  55  shous  the  effect  of  varying  Upo  on  the  particle 
density  distribution  in  the  freeboard,  fts  was  noted  previously, 
the  peak  density  occurs  at  the  point  where  the  largest  particles 
attain  their  naxiwn  height  above  the  bed.  fts  Upo  is  increased, 
the  curve  to  the  right  of  the  peak  becomes  shorter  and  steeper. 
This  trend  will  continue  until  the  the  effect  of  the  snaller 
particles  returning  to  the  bed  cause  the  slope  to  decrease.  Fig. 
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SS  is  a  semi-ln  plot  of  the  data  in  Fig.  55.  Table  13  confirms 


that  the  slope  is  becoming  steeper  as  Upo  is  increased  although 
the  influence  of  initial  particle  velocity  on  the  slope  of  the 
density  vs  height  curve  is  moderate. 


TABLE  13 


Upo 

Slope 

Intercept 

cm/s 

gms/cm 

gms/cm 

61.0 

-0.118 

6.058 

37.2 

-0.117 

6.758 

152.4 

-0.125 

8.500 

304.8 

-0.200 

18.20 

Effect  of  Upo  on  the  slope  of  the 
particle  density  distribution  as  a 
function  of  height  for  the 
distributions  shown  in  Fig.  SB. 


Variation  of  Jet  Velocity  (Uj) 

Fig.  57  shows  the  effect  on  maximum  particle  height  when  Uj 
is  varied.  As  in  the  two  previous  analysis',  the  influence  of 
momentum  and  drag  on  the  different  particle  sizes  is  apparent.  A 
decrease  in  jet  velocity  from  the  baseline  value  results  in  a 
large  drop  in  small  particle  height  and  a  small  change  in  large 
particle  height.  This  suggests  that  jet  velocity  has  a  strong 
effect  on  the  dispersion,  or  seperation  of  particles  of  different 
sizes  at  increasing  freeboard  heights.  Increasing  Uj  results  in  a 


stretching  effect  of  the  particle  size  distributions  above  the 


bed.  One  interesting  point  to  note  is  the  loss  of  the  rounded 
edge  at  the  heights  corresponding  to  the  smaller  particles.  This 
shous  the  extremely  large  effect  that  drag  plays  when  the  air 
velocity  is  very  close  to  the  terminal  velocity  of  the  particle. 

Figs.  S8  thru  S3  shou  the  effect  of  varying  Uj  on  the 
individual  partical  densities  at  varying  heights  above  the  bed. 
Again,  as  in  the  previous  parameter  analysis,  the  same 
distribution  shapes  can  be  seen  for  each  value  of  Uj ,  the  only 
difference  being  the  height  at  which  it  occurs. 

Fig.  S4  shous  the  change  in  particle  density  distribution  as 
Uj_is  varied.  Fig.  SS  is  a  semi-ln  plot  of  this  data  showing 
that  as  Uj  is  decreased,  the  slope  of  the  distribution  changes 
substantially.  Table  14  also  shous  this  effect. - 


TA6LE  14 


Uj 

Slope 

Intercept 

cm/s 

gms/cm 

gms/cm 

503.6 

-0.117 

5.758 

457.  Z 

-0.Z07 

7.380 

304.8 

-0.400 

7.300 

Effect  of  Uj  on  the  slope  of  the 
particle  density  distribution  as  a 
function  of  height  for  the 
distributions  shown  in  Fig.  S5. 
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Variation  of  Jat  Duration  (tj) 


Fig.  66  is  a  graph  showing  the  effect  of  varying  jet 
duration  on  the  naxirun  heights  attained  by  the  particles.  It  is 
easily  seen  that  a  snail  change  in  jet  duration  results  in  a 
substantial  change  in  particle  height.  This  effect  indicates  that 
the  distribution  present  in  an  actual  fluidized  bed  is  probably  a 
statistical  average  of  a  rapidly ^fluctuating  particle  distribution 
uhich  is  control ed  by  the  durations  of  the  jets  fron  the 
neighboring  bubble  eruptions  as  uell  as  local  eruptions. 


Figs.  57  thru  72  show  the  effects  on  the  individual  particle 
density  distributions,  at  different  freeboard  heights,  as  a 
function  of  jet  duration.  These  figures  also  shou  the  drastic 
change  in  density  distributions  caused  by  changes  in  jet  duration. 

Fig.  73  shous  the  particle  density  distributions  above  the 
bed  as  affected  by  jet  duration.  Fig.  74,  uhich  is  a  seni-ln 
plot  of  the  data,  shous  the  drastic  changes  in  the  decay  slopes. 
Table  15  also  shous  this  drastic  variation  in  slope. 


TABLE  IS 


tj 

Slope 

Intercept 

cm/s 

gms/cm 

gms/cm 

10.0 

-1.071 

1Z.50 

Z0.0 

-0.117 

6.758 

30.0 

-0.0Z0 

6.Z34 

Effect  of  tj  on  the  slope  of  the 
particle  density  distribution  as 
function  of  height  for  the 
distributions  shoun  in  Fig.  74. 


Variation  of  Particle-Distribution  in  Bed  Mass 

k  Fig., ,  75  shows . the _ tuo  bed  particle  distributions  used  in 
this  analysis.  The  determination  of  these  distributions  is 
descibed  earlier  in  this  uork. 


Figs.  76  thru  81  shou  the  effects  of  varying  the  bed 
particle  distribution  on  the  individual  particle  distributions  at 
different  heights  above  the  bed  surface.  Using  the  sieve  particle 
distribution  instead  of  the  image  analysis  distribution  produces  a 
shift  in  the  mean  particle  distribution  towards  the  smaller 
particles.  The  distributions  resulting  from  the  sieve  data  do  not 
shift  as  much  as  the  image  analyzer  data  uhen  the  freeboard  height 
is  increased.  The  small  shift  in  mean  diameter  exhibited  by  these 
plots  is  consistant  with  the  distributions  listed  in  Appendix  J. 


Fig.  8Z  shows  the  particle  density  distributions  in  the 


Appendix  E  and  ,J.  Interval  of  10  Microns. 
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Fig.  81  Relative  Particle  Number  vs  Particle 
Diameter  as  a  Function  of  Bed  Mass. 
Freeboard  Height  of  31  cm. 


freeboard  for  the  two  bed  rase  conditions.  The  effect  of  the 


sieve  distribution  (which  has  a  nuch  lower  average  particle 
dianeter)  is  to  decrease  the  initial  rate  at  which  particle  rass 
returns  to  the  bed.  This  is  dbe  to  the  snaller  particles  having  a 
nuch  snaller  terninal  velocity.  This  will  not  only  raise  the  bulk 
of  the  particles  to  a  higher  height,  but  will  also  increase  the 
tine  required  for  the  particles  to  return  to  the  bed. 


Conparison  of  Model  with  Exoerincntal  Results 

As  discussed  in  the  previous  section,  the  baseline  paraneters 
were  selected  as  the  experimental  conditions  present  in  the  bed 
when  sanpling  at  the  highest  Uo  setting.  Table  16  lists  these 
baseline  paraneters  again  for  review. 


TABLE  16 


Uo  =  57.9  cn/s  (1.9  ft/s) 

Upo  -  97. Z  cn/s  (3.19  ft/s) 

Uj  -  609.6  cn/s  (Z0  ft/s) 
tj  =  Z0  ns 

Baseline  paraneters  used  in  conputer  nodel . 


Fig.  83  shows  the  particle  density  distribution  predicted  by 
the  nodel  using  the  baseline  conditions.  Fig.  84  is  a  seni-log 
plot  of  this  sane  data.  A  conparison  of  the  slope  of  the 


distribution  lino  after  the  peak  value  in  Fig.  84  uith  the  slope 
obtained  fron  the  experinental  data  shous  very  good  agreenent 
between  the  two.  This  data  is  also  given  in  Table  17  along  uith 
the  data  obtained  for  the  other  experinental  data  sets  and  is 
shown  plotted  in  Fig.  85.  The  slopes  for  the  nodel  data  show  a 
steady  increase  in  slope  as  Uo/Unf  is  decreased.  This  trend  is 
also  followed  by  the  experinental  data  but  is  not  as  snooth.  The 
slope  of  the  experinental  data  is  also  changing  faster  than  the 
slope  of  the  nodel  data.  One  reason  for  this  is  that  all  the  nodel 
calculations  were  performed  using  the  sane  value  for  jet  velocity 
and  duration  of  689. S  cn/s  and  Z0  ns  respectively.  The  nagnitudes 
of-,  these:  values  should  decrease  with  decreaseing  Uo/Unf,  but  the 
values  to  be  used  for  these  different  conditions  could  not  be 
determined. 


TABLE  17 


Uo/Unf 

Experinental 

Slope 

gns/cm 

Model 

Slope 

gns/cn 

□iff 

Z 

3.81 

-0.1097 

-0.117 

8.6 

3.17 

-0.10Z9 

-0.1Z3 

19.5 

Z.59 

-0.1181 

-0.1Z9 

9.Z 

Z.3Z 

-0.1399 

-0.1 3Z 

5.6 

Comparison  of  slopes  for  the  particle  density 
distributions  above  the  bed  as  derived  from 
the  experimental  data  and  the  computer  model. 


The  slope  shoun  in  Fig.  84  for  heights  between  approximately 
4  and  18  cn  is  totally  different  fron  the  slopes  listed  in  table 
17.  This  is  also  the  lower  range  in  which  the  experinental  data 
was  collected.  As  was  shown  earlier,  when  the  particle 
distribution  analyzed  included  particles  fron  88  -  1870  nicrons, 
the  peak  was  located  at  11  cn.  This  is  the  height  where  the 
largest  particles  attain  their  naxinun  height  and  begin  to  fall 
back' to  the' beef.  It' does  not  seeHr  likely  that  including  the  feu 
particles  present  above  1070  nicrons  would  totally  account  for 
this  discrepancy.  Other  factors  which  the  nodel  has  not  accounted 
for  which  nay  explain  these  discrepancies  is  the  effect  of  varying 
jet^  velocities  and  durations ; which  uereki shown  in  chapter  V  to  vary 
greatly,  and  the  effect  of  particles  whose  velocity  vector  is  not 


perpiodicular  to  the  bed  surface 


CHAPTER  VII 


CONCLUSIONS  ftNO  RECOMMENDATIONS 


Conclusions 

Several  important  conclusions  can  be  nade  concerning  the 
particle  trajectory  computer  Model  and  the  experinental  data. 

-  -  ft  particle  traj ectory  nodal ,  even  as  simple  as  the  one 

discussed  in  this  work,  produces  results  which  closely  predict 
-  several  aspect par ticie'  acti vlrty  uittrin  a  fluidized  bed.  The 
nore  important  predictions  include  the  particle  density 
distribution  in  the  freeboard  and  the  height  distribution  of 
particles.  The  results  obtained  with  the  current  nodel  indicate 
that  further  work  on  improving  the  particle  trajectory  nodel  to 
include  particle-particle  interactions  and  paticle  velocities  not 
perpendicular  to  the  bed  surface  is  very  desirable. 

Z.  The  sensitivity  analysis  indicates  that  the  jet  duration 
and  jet  velocity  are  critical  parameters  in  determining  particle 
loading  conditions  in  the  freeboard.  Since  these  are  constantly 
changing  from  one  bubble  eruption  to  the  next,  a  statistical 
distribution  will  be  required  to  accuarately  model  freeboard 
particle  activity.  The  developement  of  a  basic  model  to  predict 


US 


the  jet  duration  and  velocity  as  a  function  of  bubble  size,  Uo,  Umf 
etc  is  therefore  needed. 

3.  Particle  distribution  analysis  of  experimental  samples  at 
increasing  heights  above  the  bed,  show  the  presence  of  large 
particles.  .  Many  of  these. are  above  the  maximum  calculated 
trajectory  heights  for  these  particle  sizes.  This  indicates  that 
either  the  jet  velocity  or  jet  duration  carried  these  particles  to 
these  heights,  or,  particle-particle  collisions  are  present  in 
numbers  great  enough  to  be  important  in  the  analysis. 

The -sampling -apparatus  designed  and  built  for  this  study 
appears  to  operate  satisfactory.  The  data  obtained  correlates 
with  uork  done  by  other  researchers  and  with  the  computer  model. 
Houever,  additional  uork  is  needed  to  correlate  sample  size  and 
particle  size  distribution  uith  freeboard  height. 

Recommendat ions 

1.  Further  experinental  data  is  needed  at  higher  velocities 
and  higher  heights  above  the  bed  using  the  sampling  apparatus 
designed  in  this  uork.  This  will  provide  additional  data  to 
evaluate  the  operation  of  the  sampling  apparatus  and  the  computer 
model.  Also,  a  scale  capable  of  measuring  quantities  of  samples 
less  than  0.01  grams  and  a  particle  removal  system  which  removes 
the  particles  uith  less  error  from  the  sample  trap  is  necessary 


for  accurate  uork. 


2.  A  sampling  device  which  hes  a  shorter  trap  height  should 


be  designed.  The  present  sample  trap  is  very  direction  oriented 
and  samples  particles  traveling  only  in  a  narrow  range  from  the 
vertical.  This  would  allow  more  accurate  work  to  be  done  in  the 
splash  zone  where  particles  are  more  likely  to  be  traveling  in 
directions  other  than  vertical. 

3.  The  computer  model  needs  to  have  encorporated  in  it,  a 
statistical  distribution  model  for  jet  velocity  and  jet  duration. 
This  will  allow  the  determination  of  the  effects  of  varying  jet 
velocity  and  duration  on  particle  distributions  in  the  freeboard. 

4.  A  correlation  for  jet  velocity  and  duration  as  functions 
of  bubble  diameter,  Uo,  and  Umf  should  be  determined  and  included 


in  the  model. 
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The  calculation  of  the  total  Moment  of  inertia  (I)  of  the 
paddles  can  be  broken  up  into  three  seperate  calculations.  First, 
the  Moment' of  inertia- (IT)  of  the  tub  aluminum  cylinders  used  to 
mount  the  paddle  arms  to  the  solenoid  shaft  is  calculated.  The 
second  calculation  (IZ),  accounts  for  the  moment  of  inertia  of  the 
paddles  themselves,  uhich  are  constructed  of  a  foam,  bassuood  and 
epox y- 1  ami netei:XFig^!~ i F i nal ly^Tthermonent  of  inertia  (13)  of 
the  harduood  mounting  ends  on  the  paddles  is  calculated. 

The  following  equations  uere  used  to  calculate  the  mass 
moments  of  inertia; 

ii  =  E._H  ri 

2  (Al> 

IZ  =  L  U  (L2  +  w2)  (PI  tl  +  PZ  t2?  . 

12  (AZ) 

13  =  2  L  U  t  P  <L2  +  U2) 

12  (A3) 


Table  A.l  is  a  listing  of  parameters  reauired  for  the  moment 


TABLE  A. 1 


Component  MKS  units  English  units 


ftliiainup  Cylinder 


Density  A1  £Pl'  " 

2780  " 

kg/m 

0.1004 

lbf/in 

Radius  Irl 

0.0171 

m 

0.675 

in 

Height  I  HI 

0.0381 

n 

1.5 

in 

Paddles 

Density  Foa«_IPl]  -  ■  —  r 

50 

kg/m 

0.00181 

lbf/in2 

Density  Basswood  lP2Tr~“ 

400 

kg/n* 

0.0144 

lbf/in3 

Thickness  Foam  I  tli.  5 

.0*0102., 

m  . . . . 

0.4 

in 

Thickness  Basswood * C  tZ T'^ 

0.001 S3 

n  ----- 

0.0625 

in 

Length  Foam  I  LI] 

0.5842 

n 

23.0 

in 

Length  Basswood  (LZ1 

0.5842 

n 

23.0 

in 

Uidth  Foam  (Ull 

0.0305 

n 

1.2 

in 

Uidth  Basswood  CUZ1  : 

0.0305 

m 

1.2 

in 

Hardwood 

Density  Maple  [PI 

650 

kg/n* 

0.0235 

lbf/in 

Thickness  Maple  ttl 

0.0102 

m 

0.4 

in 

Length  Maple  CL] 

0.0432 

n 

1.7 

in 

Uidth  Maple  IU] 

0.0348 

n 

l.Z 

in 

Listing  of  paddle  components  and  parameters. 


.  .s.  A 


Equation  (A.l)  is  used  to  calculate  the  nass  moment  of 
inertia  for  a  right  circular  cylinder  rotating  about  its  Z  axis 
(Fig.  A-Z).  The  height  H,  accounts  for  both  the  upper  and  the 
louer  solenoid  shaft  cylinders.  Equations  (A. 2)  and  (A.3> 
determine  the  nass  nonent  of  inertia  for  a  rectangular  prisn 
rotatingaboutitsiX-V  control  dal  axis~(Fig.  A-3).  The  length  L 
-in  equation  AZi accounts  f or;both_the  upper  and  louer  paddle 
lengths.  The  factor  of  Z  in  equation  (A. 3)  is  because  both  the 
upper  and  louer  harduood  sections  are  identical  and  can  be 
combined.  By  inserting  the  respective  values  from  Table  A.l  into 
equations  (A.  1 ) ,  ( A. Z ) ,  and_(A. 3) ,  the  value  for  each  inertia  ...  . 
component  can  be  calculated.  Adding  these  components,  the  total 
mass  moment  of  inertia  applied  to  the  solenoid  by  the  paddles  is 
determined. — These  results' are  listed  in  Table  A.Z. 
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TABLE  A.Z 


Moment  of  inertia  components: 

MKS  units  English  units 

II  =  1.422  (-05)  kg  m*  4.911  (-02)  Ibf  in 

r-'IZ  =  5.8Z3  (  -04)  kg-**-  1.982  (00)  '  lbf  in 

-13  =  2;5S&:<~0G»k*;*W  G.919  (-03)  lbf  in 

Total  moment:  I  =  II  +  12  +  13 

I  =  5.991  (-04)  kg  n*  2.037  (00)  lbf  in 


Calculated  values,  for^, mass  moment  of  inertia  of 
paddles,  and.  paddle'  components.  „  II  is  inertia  __ 
value  far:  th£ 'aluminum  cylinders.  12  i s  the 
inertia  value  for  the  paddles  themselves  and  13 
is  for  the  harduood  end  pieces. 


APPENDIX  B 


Error  Determination  of  Vacuum  Collection  System 

“—^"~'Tes'ts~ofr  the-vacuum  sample- removal~system  indicated  that  sorts 
'  particles- remained- in  the' sample- trap.  These  particles  were  ~ 
located  in  the  corners  of  the  trap  where  the  equalizing  air  stream 
could  not  agitate  then  enough  to  move  then  into  the  vacuum  stream. 
As  a  result.lt  became  necessary  to  determine  to  what  extent  these 
^remaining:  particles'  affected  theaccuracy  of  the  sample  attained.- 

The  procedure  used  involved  placing  samples  of  known  weight 
and  particle  size  distribution  (ecpiivalent  to  the  bed  material) 
inside  the  sample  trap.  Sample-  sizes  of-  10,  15,  Z0,  and  Z5  grams 
were  used.  The  sample  particles  were  then  vacuumed  out  and  their 
weight  determined.  The  difference  in  weights  of  the  samples  were 
then  calculated  along  with  the  percentage  differences.  These 
results  are  listed  in  Table  B.l. 

An  analysis  of  the  results  listed  in  Table  B.l  indicates  that 
the  average  difference  between  the  two  sample  weights  is  only  0.5Z 
Z  while. the  maximum  difference  observed  was  0.93  X.  It  can 
therefore  be  concluded  that  the  particles  remaining  in  the  sample 
trap  do  not  significantly  affect  the  accuracy  of  the  sample 


attained. 


il__ '.Que.to'.the  high  .velocity  of .the  .pact icles  uith  respect  to  the 
size  of  the  sample  trap,  it  uas  important  to  ensure  that  the 
sample  trap  uas  closed  in  a  very  short  period  of  time.  To  measure 
the  closure  time,  the  set-up  shoun  in  Fig.  C-l  uas  used. 

v=>;?v  ;An  electric  eye, "  consisting; of-  a  light  emitting  diode  and  a  ~ 
photo  transistor,  uas  used.  ~  The  diode  uas  pouered  by  a  G  volt 
battery  and  the  photo  transistor  uas  uired  directly  to  an 
oscilloscope.  The  electric  eye  circuit  uas  then  placed  at  a 
distance  of  0.ZG7  m  (10.5  in)  from  the  pivot  point  of  the  paddles. 
At  this  position,  the  closing  time  calculation  involves  only  a 
simple  proportion  relationship  betueen  the  paddle  and  the  sample 
trap  uidths.  The  eye  uas  also  positioned  as  close  to  the  trailing 
edge  of  the  fully  closed  paddle  as  possible  to  minimize  error. 

The  oscilloscope  sueep  uas  set  to  trigger  off  of  the 
initial  change  of  state  from  the  photo  transistor  when  the  paddle 
first  eclipsed  the  light  beam.  The  resultant  traces  were  then 
analysed  to  determine  closure  time.  Time  t=0  uas  set  equal  to  the 
initiation  of  the  trace  on  the  oscilloscope.  Time  t=tl  uas 
defined  as  the  point  at  uhich  the  trace  begins  its  excursion  back 


f* - ZS.7 - H 

(b)  Side  Vieu 


Fully  Activated 


Diagram  of  Closure  Tine  Determination  Set  Up. 
All  Dimensions  in  cm. 


Fig.  C-l 


to  the  base  voltage  as  shoun  in  Fig.  C-Z.  By  using  the  ratios 


ML  =  42 

tl  tZ  (C.l) 


Uhere-' 

U1  =  width  of  the  paddle 
uf  isrrpls  -b*i  UZ  =  width  of  sample  trap-—-  -  . 

-  -  :  ~ t l=-t i me  delay-  measured  on  oscilloscope 

_•  \;r;  of  :c  tZ  =  closure  tine  of  sample  trap 


the  closure  tine  of  the  sample  trap  can  be  determined.  Several 
trials  were  run  with  the  resulting  data  listed  in  Table  C-l.  All 
of- the  trials- ware-very .consistant^aslbe  resultant. determination  . 
for  the  average- closure  tine  was  1.44  ns. 


TABLE  C-l 


U1  =  0.0305  n  (l.Z  in) 

UZ  =  0.01Z7  n  (0.5  in) 

Trial  tl  tZ 

Mo.  ns _ ns 

1  3.4  1.4Z 

Z  3.4  1.4Z 

3  3.55  1.48 

4  3.5  1.46 

_§ _ 2LS _ 

Avg  3.47  1.44 


Sample  trap  closure  data-'  Ui  is  the  uidth  of 
the  paddle,  UZ  is  the  uidth  of  the  sample  trap 
and  tine  tl  is  eclipse  tine  of  paddle  through 
light  bean.  Tine  tZ  is  closure  time  of  Sample 
trap  as  calculated  using  equation  (C.l). 


I? 


The  goal  of  this  design  uas  to  shut  the  sample  trap  in  about 


1  millisecond  (ms),  ft-t  this  speed,  a  particle  traveling  at  10 
meters  per  second  (m/s>  would  travel  1  cm  or  approximately  11  X  of 
the  sample  trap  length.  A  velocity  of  10  m/s  is  at  the  upper 
limit  of:  the  particle  velacityidistribution  and  would  account  for 
only  a  small  percentage" of:.  the  particles- ejected  by  the  bed-  The 
majority  of"  the  particles,  for'  the  fluidization  conditions  used  in 
the  bed,  have  an  average  ejection  velocity  of  1-Z  m/s,  based  on 
data  from  George  and  Grace  (81. 


To  determine  the  minimum  required  torque  output  of  the 
solenoid,  two  analysis  were  performed.  The  first  analysis 
calculates  an  average  torque  required  to  produce  the  desired 
velocity.  The  second  analysis  uses  the  torque  data  from  the 
selected  solenoid  and  calculates  the  expected  closure  time  and 
suing  time  for  the  paddle  arms. 


».*%/ .  * 


•y  ' 


V  ’-O 


To  close  the  sample  trap  in  1  ms,  the  angular  velocity  of  the 


paddles  at  closure  can  be  calculated  as: 


(0.1) 


where: 


u  =  Angular  velocity  of  paddles 
i  UZ_ :  =  Uidth  of  sample  trap 
R  =  Paddle:  pivot- to; sample  trap  distance 
l  tZ_  =:  Closure  time  of  sample  trap  - 


Using  the  values  of  U2=  0.00197  m,  R=  0.0413  m,  and  tZ=  1  ms, 
equation  (0.1)  results  in  an  angular  velocity  of  47.70  rad/sec. 

The  radial  acceleration  required  to  attain  this  velocity  through  a 
deflection  of  **74  radians  (45  Deg),  is  calculated  by: 


Z 

a  -  u 
Z  0 


(0.Z) 


where:  a  =  Angular  acceleration  of  paddles 

u  =  Angular  velocity  of  paddles 
9  =  Angular  deflection  of  paddles 


Using  the  result  of  equation  D1  and  8=1*74  radians,  equation  DZ 
gives  the  required  angular  acceleration  of  the  paddles  as  1448.5 


rad/ s/s.  Combining  this  result  uith  the  result  for  the  total  mass 


moment  of  inertia  from  Appendix  A,  the  required  average  torque 
output  of  the  solenoid  can  be  calculated  as: 


T  =  1  a 


(03) 


-  where-vb  r.- a  =  Angular  acceleration  of  paddles 

-=  Mass  moment  of  inertia 

-  .  T  *  Required  solenoid  torque 


The  resultant  value  for  the  torque  (T)  is  0.87  N  m  (7.7  Lbf  in), 
ft  geometric:  .average  of  the  selected  solenoid's  torqueoutput.as 
shown  in  Fig.~-  D-la,~  indicates  an  average'  of  1.0  M  m  (8.8  Lbf  in). 
Therefore,  the  selected  solenoid  has  sufficient  torque  output  to 
achieve  the  desired  paddle  velocity. 


SECQNQ  RNQLY.5IS 

For  the  second  analysis,  it  will  be  assumed  that  the  torque 
output  of  the  selected  solenoid  can  be  modeled  as  a  linear  spring. 
Fig.  D-lb  shous  the  torque  output  of  the  solenoid  as  a  function 
of  angular  displacement.  From  Fig.  D-lb,  a  spring  constant  of  k: 
0.72  _N  m/rad  ( 0 .  U_  Lbf  _Ln/ deg ) - can  be  used  to  approximate  the 
torque  curve.  The  spring  constant  is  derived  from  the  slope  of 
the  torque-angular  displacement  curve. 


The  system  can  then  be  modeled  as  a  simple  rotational  mass  - 
spring  system.  The  differential  equation  of  notion,  general 
solution,  and  boundary  conditions  are: 


□iff  Eqn: 


I  9  +  k  0  =  9 


(0.4a) 


■  ■=  Gen  Soln~  '■  ~  8-.=:  Cl:  sin(n  t)  C2'  cos(n  t)  (0.4b) 

.  -  t  '  ~  Z  —  — ,  ;  .  ■  .  

-  — - uhere:— m-^ — -  K  - - — -  -  - 

I 

b.c. :  1)  t  =  0-  0=0 

Z)  t  =  9-  0  =  1 .65  rad 


Solving  Equation  (0.4b >=. using:, the: :gi.vapi  boundary  equations 


the' follouing  expressions  for  any  angular  displacement  (6) 
angular  velocity  (8): 


,  gives 
and 


0  =  1.65  cos(m  t)  (0.5a) 

0  =  -1.65nsin(nt)  (0.5b) 


Knouing  that  the  solenoid  stroke  is  1174  radians,  the  total  solenoid 
actuation  tine  can  be  obtained  from  equation  (0.5a)  by  setting  0  = 
<1.65  -TT/4).  The  total  actuation  time  uas  determined  to  be  Z9.Z 
ms.  Inserting  this  value  for  time  into  equation  (0.5b),  a 
rotational  velocity  of  48.6  rad/s  uas  found.  This  velocity  is 


very  close  to  the  desired  velocity  obtained  from  equation  (0.1). 


The  desire  to  determine  whether  the  debris  of  a  specific 
bubble  has  been  trapped  by  the  sampling  device  requires  the 
determination  of  the  total  actuation  time  for  the  apparatus. 
Knouing_thiai^yalue,^the:_delay  between.  Uhe  detection  of  a  bubble,, 
and  the  closure  of?  the  sample: .trap  can  be  determined,  and  hence, 
whether  the  par  t_i_cl  es  _f rom_the_detected  bubble  were  within  the 
vicinity  of  the  trap. 

4%.u  flppendixX.  datermlned-that  the  actu.al  closure  time  for  the 
sample  trap  was -1:44  ms  .:t:  Using  this  value  in  equation  ( D.  1 ) ,'  "a 
rotational  velocity- at -the  end  of  rthe  stroke  of  only  33. 1Z  rad/ sec 
Us  calculated.  Because  the. mass  moment  of  inertia  for  the 
solenoid  was  unknown  in  the  beginning,  the  analysis  ignored  it. 
Uith  the  determintion  of  the  actual  closure  tine  in  Appendix  C, 
the  total  moment  of  inertia,  including  the  solenoid  can  be 
approximated.  8y  iterating  equations  (0.5a)  and  (0.5b)  uith  0  = 
33.12  radians/sec,  a  new  inertia  value  of  1.27  (-03)  kg  n  is 
determined.  Using  this  inertia  value  in  equation  (0.5b),  the 
total  actuation  tine  is  42.6  ms.  This  increase  in  actuation  tine 
of  13.5  ms  is  relativily  close  to  the  manufactures  quoted 
actuation  time  for  the  solenoid  of  12  ms  and  suggests  that  the 
values  are  within  reason. 

Table  0.1  provides  a  summary  of  the  results  obtained  from  the 


dynamic  analysis. 


TABLE  0.1 


Mass  moment  of  inertia: 

Paddles  <!>  5.991  (-04)  kg  m] 

2.037  (00)  lb  in 

Xs'h  ■-*'  --n?  r  7  It  >  Paddled7+- Solenoid  fit)  1.27  (-03)  kg 

.  .  ‘  _4.32  (93)  i‘r  v.''  4.32  (00)  lb  in 

'll _ 1  Closure  tine: 

-  '  .  7  Sample  trap  (t2)  1.44  ms 

Total  suing  (t3)  42.8  ms 

vSummary„o€- dynamic  analysis  results. 


A  sample  of  particles  uas  obtained  from  the  fluidized  bed 
after  the  bed  had  been  operating  for  several  hours.  Sample  sizes, 
ranging  from  1_. Z  to  Z.8  kg,  were 'then  taken  from  the  central  area 
of  the  bed.  Each  sample  uas  then  sieved  through  a  series  of  14  US 
standard  uire  mesh  sieves  using  a  Tyler  Industrial  Products  Model 
RX-Z4  portable  sieve  shaker  for  30  minuites.  The  contents  of  each 
sieve  uas  then  ueighed,  using  a  Torsion  Balance  Co  TORBAL  scale, 
to  an  accuracy  of  0.01  gram.  The  resulting  average  particle  size 
distribution  is  listed  in  Table  E.l  and  shown  in  fig.  E-l. 


TABLE  E.l 


US 

Seive 

Weight 

Percent 

Cumulative 

Standard 

Size 

In  Seive 

of  Total 

Percent 

Seive  No. 

Ufl 

grans 

Sanple 

18 

1000 

8.49 

0.31 

0.31 

20 

850 

Z4.50 

0.S0 

1.Z1 

38 

600 

101.96 

3.73 

4.95 

35 

500 

110.85 

4.06 

9.01 

45 

355 

176.35 

6.46 

15.46 

sa 

ZS7 

430.77 

15.78 

31.24 

68 

Z50 

385.59 

14. 1Z 

45.36 

.70 

Z1Z 

395.32 

14.48 

59.84 

80 

180 

469. 20 

17.18 

77.03 

100 

149 

339.93 

1Z.45 

89.48 

1Z0 

1ZS 

1ZS.01 

4.58 

94.06 

140 

106 

109.45 

4.01 

98.06 

170 

90 

46.67 

1.71 

99.77 

Z00 

75 

5.91 

0.ZZ 

99.99 

<75 

0.30 

0.01 

100.00 

Particle  Specific  Gravity  =8.1 

Average  particle  size  distribution  of  bed  naterial  in  grans 
and  percentage  of  total  weight  using  sieve  analysis. 


APPENDIX  R 


Mean  Bed  Flou  Velocity  Determination 

To  determine  the  mean  velocity  of  air  flowing  through  the 
bed,  the  ASHE  report  on  fluid  meters  [171  uas  used.  Based  on  this 
papery  equation  (F.l)  uas  used  for  determining  the  mass  flou  rate 
of  air  through  the  bed. 


U  =  0.09970Z  ( 


K  Y  F  d  )  (  6  dP  ) 


(F.l) 


where: 

U  *  Mass  flou  rate  (lbm/sec) 
d  =  Orfice  diameter  (inches) 

K  =  Flou  coefficient 

F  =  Thermal  expansion  factor  (1  for  air) 
Y  =  Expansion  factor 

dP  =  Differential  pressure  (inches  water) 
C  -  Density  of  air  (ahead  of  orfice) 


The  flou  coefficient  (K)  is  a  function  of  Reynold’s  number,  the 
orfice  diameter  (d),  and  the  pipe  diameter  (D).  For  the  pressure 
tap  configuration  used  in  the  MIT  atmospheric  fluidized  bed  (1-D  / 
1/Z-O)  the  flou  coefficient  is  calculated  using  equation  (F.Z). 


I 


171 


flTH? 


<F.Z) 


where: 

Re  =  Reynolds  No. 

B  =  d/0  =  (orfice  die  /  pipe  dia) 

Ko  =  ( 0.S014  -  0.0135Z  0  ^  ) 

+  <0.3760  +  0.07ZS7  D~S 

(  9.99025  +  B*  +  1.5  B'‘  ) 

D  B  +0.000Z5  □ 

b  =  (0.000Z  +  0.0011  ) 

□ 

+  <0.0038  +  0.0004) 

D 

<  B*  +  <16.5  +  5  0)  s'4  ) 


The  expansion  factor  <Y)  is  detemined  from  equation  <F.3)  and  is 
a  function  of  the  diameter  ratio  (B),  the  ratio  of  specific  heats 
<s),  and  the  ratio  of  the  differential  pressure  to  inlet  pressure 


Y  =  1  -  <0.410  +  0.350  BH  )  dP 

PI  s  <F.3> 

where: 

PI  =  Inlet  pressure 

s  =  Ratio  of  specific  heats  (1.4  for  air) 


Mass  flou  rate  is  a  function  of  velocity,  and  therefore,  by 


definition  a  function  of  Reynold's  number.  Equation  (F.4)  is  used 
to  determine  mass  flou  rate  as  a  function  of  viscosity  and 
Reynold’s  number. 


U  =  u  02 

C  Re  (F.4) 


uhere  : 

u  =  Absolute  viscosity  (lbm/ft  sec) 
C  =  Constant  (15.28) 


This  requires  that  the  solutions  of  equation  (F.l)  and  equation 
(F.4)  be  iterated  until  the  Reynold’s  numbers  converge.  The 
determination  of  mean  air  velocity  is  obtained  from  the  mass  flou 
rate  (U)  using  equation  (F.5). 


V  =  U 

e  r 


(F.5) 


uhere: 

a  * 

v  =  Density  of  air  (lbm/ft  ) 
ft  =  Area  of  bed  (ft1) 
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To  perform  these  calculations,  a  computer  program  was  used.  The 


program  (APPENDIX  8)  is  written  in  HP  BASIC  2.0  for  running  on  an 
HP  9818  series  200  micro-computer.  Convergence  usually  required 
four  to  seven  iterations  and  should  be  correct  to  within  0.SX. 
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i  •*  .  lo*  •  1.  "  •». 


■ 


«  -  ■  -  mm 


. 


20 

30 

*0 

se 

50 

70 

80 

80 

100 

110 


. MAIN  Program . 

Metering  of  gases  oy  naans  of  the  ASrtE  souare-eageb  orfice  uith 
1-0  1/3-0  taos.  Reference  Fluid  Meters  Their  Theory  anc 
Application,  ASMS  Report  Sin  edition,  ASME  New  Yort.,  NY,  1371 
Progran  nusi  be  altered  for  Orfice  dianeter  other  than  7.G71  in 
and  pipe  diameter  other  than  10.03  in.  (see  line  130) 


INTESES  Answer 

REAL  01 ,03 ,U,K , Y1 ,Pi ,Pla,71 ,S,?_del .Sete,nrea_Bed, Vmu,fieynoids_no .Rhc ,Rhol 


,Patn, Velocity 


120 

Oin  Hone* C 2 1 , Clears! 2 1 

130 

0ATA  7.071,10.02.1.4,11.511,50000. ,75 

14.0 

REA0  02.01 , S, Area_bed.Reynolds_no ,Patn 

150 

Clear*-CHRS( 255 )4CHRs(  75 ) 

!  CLEAR  screen 

150 

Hone*-CHRS( 255 )4CHRS( 34 ) 

!  HOME  screen 

170 

l 

180 

!  Input  variables 

180 

i 

200 

INPUT  "Enter  static  pressure  PI  (c«  Hg)  =  ",P1 

210 

INPUT  "Enter  pressure  drop  dP  (inches 

water):  ",P_del 

220 

INPUT  "Enter  air  tenperature  (Degrees 

F>:  *  ,T1 

230 

i 

240 

!  Compute  data 

250 

! 

250 

OUTPUT  2: Hone*; 

!  Hone  display 

270 

OUTPUT  2; Clear*: 

!  Clear  display 

230 

Beta-02/01 

!  Calculate  beta 

290 

Ple-Pl+Patm 

f  Absolute  pressure 

300 

Pl-(Pla)/2.S4 

<  Change  to  inches 

310 

T2-T1+4S9.S7 

!  Convert  to  R 

320 

Rho-r NOens ity(Patn,T2) 

!  Air  density  bed 

330 

Rhol-FNOensity(Pla,T2  > 

!  Air  density  upstream 

340 

Vnu— rNViscosi ty( T2 ) 

!  Air  viscosity 

350 

Yl-FNExpan<8eta,Pl,S,P_del ) 

!  Expansion  factor  Y 

350 

K-FNFlow_coef f (Reynolds_no.8eta,Dl ) 

!  Flow  coefficient  K 

370 

U-FNflas  s_f  1  ou<  02 ,  K ,  Y 1 ,  Rho  1 ,  P_de  1  > 

!  Mass  flow  rate 

380 

Rd-FNReyn(U,Vnu,Q2 ) 

!  Reynolds  4  calculated 

390 

IF  AfiS (Rd-Rey no! ds_.no  )>'l .  THEN 

!  Check  accuracy 

400 

Reynolds_no-Rd 

410 

SOTO  350 

420 

ENO  IF 

430 

Uelocity-U/(Rho-Area_bed ) 

!  Velocity  calculation 

440 

PRINT  "Velocity-  ":Veloeity;"  ft/s" 

I  Output  results 

450 

INPUT  "Enter  (1)  to  continue,  (0)  to 

stop* .Answer 

450 

IF  Answer-0  THEN  STOP 

470 

IF  AnswerOl  THEN 

430 

PRINT  "Value  must  be  either  l  or  0: 

try  again’ 

430 

SOTO  450 

500 

ENO  IF 

510 

SOTO  150 

520 

ENO 

530  ! 


340  i  Ene  MAIN 
350  1 


Begin  FUNCTIONS 


see  FNE.NO 

sie  ! 

520  ! 

630  OEF  FNFlow_coef f (Reynoids_no,8eta,Dl )  !  Flow  coefficient  K 

S40  Ko».6014-.313S2*O1 -( -  .25 )+( . 376+. 37257*01 *(-.25) )*< .3002S/<Ol*ai*8eia*ee 

ta+ .  00025-01  >+8eta*4+l ,5*Seta* 16  ) 

650  3«.3002+.0011/Ol+< .0038+. 0004/01 >•( Beta*9eta+< 16.5+5*01 >»8eta"lS ) 

I  So0  K»!<o+1000*8/SOR(Seta*Reynoid3_no  > 

S70  RETURN  K 

S30  FNENO 
550  ! 

730  ! 

713  OEF  FNHass_f lou( C2,K ,Y1 ,Rhol ,P_dei )  !  Hass  flow  rate  U 

723  U«.099702*O2*O2*K*Yl»SQR(Rhoi»P_del > 

730  RETURN  U 

740  FNENO 
750  ! 

7S0  ! 

770  OEF  r NReyn< U , Umu , 02 )  !  Reynolds  4  calculated 

780  Rd«lS.28*U/<Umu*02  ) 

750  RETURN  Rd 

800  FNENO 
310  ■  > 

820  ! 

830  OEF  FN0ensity(Pain,T2 )  !  Calculate  air  density 

840  Rho“.522406*Pat«/T2 

850  RETURN  Rho 

j  360  FNENO 

870  ! 

880  ! 

890  OEF  FNUiscosity<  T2  ) 

900  Umu«7.303E-7*T2*l. 5/(72+198. S) 

910  RETURN  Umu 


!  Calculate  air  viscosity 


APPENDIX  H 


This  Appendix  contains  the  parts  list  for  the  sampling  system 
in  three  tables.  Table  H.l  is  a  list  of  components  for  the  sampling 
apparatus.  Table  H.Z  is  a  list  of  components  for  the  electrical 
circuit  uhile  Table  H.3  is  for  the  vacuum  system. 


Table  H.l 


lien  No.  Conponent 


Description 


1  Solenoid 


Z  Paddle 


3  Interface 


4  Solenoid 
Mounting 

5  Extension 
Bar 


LEOEX  Size  SS  45  Oegree  Right 
Hand  Stroke  Rotary  Solenoid 
Part  No.  S-8204-0Z9 
LEDEX  Inc. 

801  Scholz  Or,  P.0.  Box  427 
Vandal ia.  Oh  45377 
(513)  838-3521 

Foan,  Basswood  epoxy  laminate 
with  Maple  mounting  blocks 
0.5  X  1.2  X  11.5  inches 

Aluminum  Cylinders  uith  set  screu 
fastener.  R=  0.S7S  in  H=  1.5  in 

Aluminum  mounting  plate 
1/4  X  Z.8  X  2.8  inches 

Aluminum  bar 

1/4  X  3/4  X  12.5  inches 


6  Trap 

Mount 


Aluminum  Plate 

1/4  X  1.1  X  1.5  inches 


7 


8 


Sample 

Trap 


1/15  inch  Aluminum  plate 
Inside  Dimensions: 

0.5  X  1.5  X  3.5  inches 


Vertical 

Mounting 

Slide 


Aluminum  bar 
1/4  X  3/4  X  24  inches 
Adjustment  holes  drilled  every 
0.5  inch 


3 


Base  Tripod  Aluminum  Structure 

Structure 


List  of  components  for  sampling  apparatus. 
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TABLE  H.Z 


Item  No. 

Component 

Description 

1 

81 

Bridge  rectifier  assembly,  Silicon, 
LEDEX  Part  No.  1Z1011-001,  includes 
arc  suppressor  in  unit. 

Z 

01 

Arc  Suppressor,  not  needed  if  above 
rectifier  assembly  used. 

LEDEX  Part  No.  1ZZSSS-001 

3 

Ml 

Solenoid  (see  Table  H. 1) 

4 

R1 

Resistor,  Z  Mohm,  11,  Z  Uatt 

5 

RZ 

Resistor,  100  kohm,  IX,  1  Uatt 

e 

R3 

Resistor,  Z50  ohm,  10X,  50  Uatt 

7 

SI 

Suitch,  SPST,  Push  Button, 

10  A,  Z50  V 

8 

SZ 

Suitch,  SPST,  Toggle,  10  A,  ZS0  V 

List  of  components  for  solenoid  pouer  supply. 
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TABLE  H. 3 


Item  No. 

Component 

Description 

1 

VI 

3/8  inch  Ball  Valve 

2 

PI 

Eductor 

3 

FI 

Screen  Filter,  320  un  mesh 

4 

Cl 

Sample  Container,  Small  Plastic 
Bottle,  1  Pt 

5 

CZ 

Sample  Trap  (see  Table  H.l) 

‘  S 

LI 

1/4  inch  Polyflou  Tubing 

List  of  components  for  vacuum  system. 
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APPENDIX  I 


This  Appendix  contains  a  complete  listing  of  all  data 
obtained  during  this  study.  The  data  is  arranged 
acording  to  the  distance  of  the  trap  height  above  the 
distributor. 


Trap  Height  above  Distributor:  cm  (in) 

= 

30.96 

(12.19) 

Bed  Height  above  Distributor-'  cm  (in) 

= 

25.4 

(10) 

Trap  Height  above  bed  surface:  cm  (in) 

= 

5.56 

(2.19) 

Pi:  cm  Hg  (in  Hg) 

s 

5.8 

(2.28) 

dP :  cm  uater  (in  uater  > 

= 

1.234 

(0.486) 

Temp:  C  (F) 

= 

19.5 

(67.0) 

Ueight  of  Sample 
( grams ) 

3.16 
'  2.06 

1.88 

2.43 

2.48 

2.10 

1.84 

2.67 

2.57 

2.16 

2.16 
1.03 
2.46 

1.13  AVG  =2.09  grams 

1.21  STN  DV  =0.60  grams 


Trap  Height  above  Distributor:  cm  (in)  =  30.96  (12.19) 

Bed  Height  above  Distributor:  cm  (in)  =  26.67  (10.5) 

Trap  Height  above  bed  surface:  cm  (in)  =  4.29  (1.69) 

Pi:  cm  Hg  (in  Hg)  =  6.8  (2.68) 

dP:  cm  uater  (in  uater)  =  1.887  (0.743) 

Temp:  C  (F)  =18  (64.5) 


SAMP  t  Ueight  of  Sample 
( grams ) 

16  7.72 

17  9.12 

18  3.83 

19  4.04 

20  5. So 

21  3.88 

22  8.86 

23  4. 50 

24  6. 55  AVG  =  5.96  grams 

25  5.84  STN  OV  =  Z.1Z  grams 


SAMP  # 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 


Trap  Height  above  Distributor:  cm  (In) 

3 

37.3  L 

(14.59) 

Bed  Height  above  Distributor:  cm  (in) 

3 

25.4 

(10) 

Trap  Height  above  bed  surface:  cm  (in) 

3 

11.91 

(4.59) 

PI:  cm  Hg  (in  Hg) 

3 

5.8 

(2.28) 

dP:  cm  water  (in  water) 

= 

1.201 

(0.473) 

Temp:  C  (F) 

20.5 

(69.0) 

-  SAMP  * 

Ueight  of  Sample 

( grams ) 

-  --  • 

25 

0.62 

27 

0.55 

v'»*  28 

0.54 

29 

0.54 

-  30 

0.76 

31 

0.58 

32 

0.61 

•- .  .  33 

0.47 

34 

0.85 

AVG  =  0.6Z 

grams 

35 

0.61 

STN  DV  =  0.11 

grams 

Trap  Height  above  Distributor:  cm  (in) 

s 

37.31 

(14.69) 

Bed  Height  above  Distributor:  cm  (in) 

s 

25.57 

(10.5) 

Trap  Height  above  bed  surface:  cm  (in) 

= 

10.54 

(4.19) 

PI :  cm  Hg  ( in  Hg) 

= 

6.3 

(2.28) 

dP:  cm  water  (in  water) 

= 

1.897 

(0.747) 

Temp:  C  (F) 

= 

19 

(65.0) 

SAMP  # 

Ueight  of  Sample 

( grams ) 

36 

2.87 

37 

1.42 

38 

1.89 

39 

1.34 

40 

1.18 

41 

1 . 52 

42 

1.57 

43 

1.68 

44 

1.56 

AVG  =  1 . 56 

grams 

45 

1.57 

STN  DV  =  0.46 

gr  ams 

las 


Trap  Height  above  Distributor: 

cn  ( in) 

=  37.31 

(14.69) 

Bed  Height  above  Distributor: 

cn  ( in) 

=  Z7.94 

(11.0) 

Trap  Height  above  bed  surface: 

cn  (in) 

=  9.37 

(3.69) 

PI:  cn  Hg  (in  Hg) 

=  6.8 

(Z.68) 

dP:  cn  uater  (in  water) 

=1.887 

(0.743) 

Tenp:  C  (F) 

=  19 

(66.0) 

SAMP  %  Ueight  of  Sanple 
( grans ) 

4.3Z 
Z.48 
3.1Z 
Z.83 
4.46 
Z.37 
4.46 
4.95 

4.15  AV6  =  3.61  grans 
Z.96  STN  DV  =  0.95  grans 


46 

47 

48 

49 

-  50 

51 
SZ 

__  _  _  S3 

54 

55 


TRAP  HEI6HT  =  Z0.ZS  in 


Height  of 

Pressure 

Pressure 

Pressure 

Pressure 

Pressure 

for 

f  or 

for 

for 

Tap  above 

Sampl es 

Samples 

Samples 

Samples 

Distributor 

66-75 

86-95 

56-65 

76-85 

( i nches ) 

( cn  HZO) 

< cn  HZO) 

(cm  HZO) 

(cm  HZO) 

1.6 

5Z.9 

54.0 

56.1 

58.5 

Z.5 

48.8 

49.3 

51.7 

5Z.8 

3.6  —  - 

40.6 

-  -  -  41.3 

-  43.8 

44.7 

4.6 

34.9 

35.8 

38.9 

40.0 

5.6 

Z7.8 

Z9.7 

3Z.5 

33.6 

6.6 

Zl.Z 

ZZ.3 

Z5.5 

Z7.6 

7.6 

14.6 

15.1 

19.1 

Z0.7 

8.6 

8.3 

9.7 

12. S 

14.8 

9.6  5.S 

,  z.z 

3.5 

3.7  6.4 

8.5 

10.6 

0.0 

0.1 

1.5 

Z.4 

_  11.6 

0.0 

0.0 

0.0 

0.3 

1Z.6 

0.0 

0.0 

0.0 

0.0 

13.6  -  . 

0.0 

0.0 

0.0 

0.0 

14.6 

0.0 

0.0 

0.0 

0.0 

TRAP  HEIGHT  =  Z3.ZS  in 


Height  of 

Pressure 

Pressure 

Pressure 

Pressure 

Pressure 

for 

for 

for 

for 

Tap  above 

Samples 

Samples 

Samples 

Samples 

Distributor 

Z06-Z15 

196-Z0S 

186-195 

176-185 

( i nches ) 

(cm  HZO) 

(cm  HZO) 

(cm  HZO) 

(cm  HZO) 

1.6 

5Z.0 

53.4 

56.3 

57.7 

Z.6 

48.9 

49.6 

51.8 

5Z.Z 

3.6 

40.7 

41.5 

43.6 

44.7 

4.6 

34.7 

35.8 

38.8 

40.0 

5.6 

Z8.3 

Z9.3 

3Z.0 

33.3 

6.6 

Z1.5 

ZZ.9 

Z5.6 

Z6.6 

7.6 

14.9 

16.6 

18.8 

Z0.3 

8.6 

8.4 

9.8 

1Z.9 

13.8 

9.6 

Z.Z 

3.3 

6.  Z 

7.5 

10.5 

0.0 

0.1 

1.3 

Z.9 

11.6 

0.0 

0.0 

0.0 

0.3 

1Z.6 

0.0 

0.0 

0.0 

0.0 

13.6 

0.0 

0.0 

0.0 

0.0 

14.6 

0.0 

0.0 

0.0 

0.0 
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Trap  Height  above  Distributor:  cm  (in) 

2 

51.44 

(20.25) 

Bed  Height  above  Distributor :  cm  (in) 

2 

29.21 

(.11.5) 

Trap  Height  above  bed  surface:  cm  (in) 

= 

22.23 

(8.75) 

Static  Bed  Height:  cm  (in) 

= 

22.56 

(8.88) 

PI:  cm  Hg  (in  Hg) 

= 

8.3 

(3.21) 

dP:  cm  water  (in  water) 

= 

S.98S 

(2.75) 

Temp:  C  (F) 

2 

16.5 

(52. a) 

Weight  of  Sample 

-  (grams)  . . .  . 

1.47 
1 . 35 
1.24 
1.29 
«.  74 
1.38 
1.38 
1.47 

1.58  AVG  =  1 . 32  grams 
1.33  STN  DV  =  a. 23  grams 


SAMP  It 


75 

77 

78 

79 
88 
81 
8Z 

83 

84 

85 


Trap  Height  above  Distributor:  cm  (in) 

2 

51.44 

(20.25) 

Bed  Height  above  Distributor:  cm  (in) 

= 

27.31 

(10.75) 

Trap  Height  above  bed  surface-'  cm  (in) 

2 

24.13 

(9.5) 

Static  Bed  Height:  cm  (in) 

= 

22.56 

(8.88) 

PI :  cm  Hg  ( in  Hg) 

2 

7.3 

(2.87) 

dP:  cm  uater  (in  water) 

2 

4.910 

(1.933) 

Temp:  C  (F) 

= 

16.5 

(62.0) 

SAMP  It  Weight  of  Sample 


( grams ) 


56 

0.95 

57 

0.61 

58 

1.83 

59 

0.84 

50 

1.02 

51 

0.76 

SZ 

0.67 

S3 

1.06 

64 

0.60 

65 

0.88 

AVG  =  0.92  grams 
STN  OV  =  0.35  grams 


Trap  Height  above  Distributor:  cm  (in) 
Bed  Height  above  Distributor:  cm  (in) 
Trap  Height  above  bed  surface:  cm  (in) 
Static  Bed  Height:  cm  (in) 

Pi:  cm  Hg  ( in  Hg) 

dP :  cm  uater  (in  uater) 

Temp:  C  <F) 


s 

51.44 

(Z0.Z5) 

= 

ZS.S7 

(10.5) 

= 

Z4.76 

(9.75) 

s 

ZZ.56 

(8.88) 

s 

6.5 

( Z . 56 ) 

s 

3.Z61 

(1.Z84) 

= 

IS 

(51.0) 

SAMP  #  Ueight  of  Sample 
( grams ) 


0. 87 
0.88 
0.17 
0.04 
0.08 
0.04 
0.07 
0.13 
0.07 
0.17 


AVG 
STN  DV 


0.09  grams 
0.05  grams 


Trap  Height  above  Distributor:  cm  (in)  =  51.44  (Z0.Z5) 

Bed  Height  above  Distributor-’  cm  (in)  *  Z5.4  (10.0) 

Trap  Height  above  bed  surface:  cm  (in)  =  ZB. 03  (10.Z5) 

Static  Bed  Height:  cm  ( in)  ~  Z2.5B  (8.88) 

Pi:  cm  Hg  (in  Hg)  =  6.Z  (Z.44) 

dP:  cm  uater  (in  uater)  =  Z.581  (1.016) 

Temp:  C  (F)  =17  (B3.0) 


SAMP  #  Ueight  of  Sample 
( grams ) 


AVG 

STN  DV 


0.1Z  grams 
0.04  grams 


TRAP  HEIGHT  =  16.12  in 


Height  of 
Pressure 


Pressure 

for 


Pressure 

for 


Pressure 

for 


Pressure 

for 


"ap  above 

Sanples 

Sanples 

Sanples 

Sanpl es 

.stributor 

106-115 

126-135 

96-105 

116-125 

inches) 

Icn  H20) 

(cn  H20) 

(cn  H20) 

(cn  H20) 

1.6 

52.6 

54.4 

55.5 

58.5 

2.6 

48.7 

49.9 

51.3 

52.9 

3.6 

40.3 

.  41.8  -- 

43.3 

44.4 

4.6 

34.6 

35.9 

37.9 

"  39.6 

5.6 

27.9 

29.2 

31.3 

33.9 

6.6 

20.6 

22.8 

'  2S.2  - 

27.1 

7.6  ... 

13.8 

15.5 

18.9 

20.6 

8.6 

~~  7.8  ii.i 

9.4 

12.3 

14.6 

9.6  — '  ~ 

1.8  “ 

3.2 

5.3 

7.6 

18.6 

0.0 

0.3 

0.3 

2.5 

11.6 

0.0 

0.0 

0.0 

0.3 

12.6 

0.0 

0.0 

0.0 

0.0 

13.6 

0.0 

0.0 

0.0 

0.0 

14.6 

0.0 

0.0 

0.0 

0.0 

TRAP  HEIGHT  =  18.25  in 


Height  of 

Pressure 

Pressure 

Pressure 

Pressure 

Pressure 

for 

for 

for 

for 

Tap  above 

Sanples 

Sanples 

Sanples 

Sanples 

Distributor 

146-155 

136-145 

156-165 

166-175 

( inches ) 

(cn  H20) 

(cn  HZO) 

(cn  H20) 

(cn  H20) 

1.6 

52.7 

54.3 

55.8 

57.9 

2.6 

49.0 

49.9 

51.5 

52.6 

3.6 

41.0 

41.5 

43.3 

44.9 

4.6 

34.5 

36.1 

38.1 

39.9 

5.6 

28.0 

29.3 

31.6 

33.5 

6.6 

21.5 

Z2.5 

24.9 

26.6 

7.6 

14.9 

15.8 

18.8 

20.7 

8.5 

8.5 

9.3 

12.2 

14.5 

9.5 

Z.4 

3.4 

5.0 

7.8 

10.5 

0.0 

0. 1 

0.3 

2.5 

11.5 

0.0 

0.0 

0.0 

0.3 

12.5 

0.0 

0.0 

0.0 

0.0 

13.5 

0.0 

0.0 

0.0 

0.0 

14.5 

0.0 

0.0 

0.0 

0.0 

Trap  Height  above  Distributor:  cn  (in) 

2 

40.85 

(16.12) 

Bed  Height  above  Distributor:  cn  (in) 

2 

28.58 

(11.25) 

Trap  Height  above  bed  surface:  cn  (in) 

= 

12.37 

(4.88) 

Static  Bed  Height:  cn  (in) 

2 

22.56 

(8.88) 

PI :  cn  Hg  ( in  Hg) 

= 

8.0 

(3. IS) 

dP:  cn  uater  (in  uater) 

= 

5.858 

(2.7) 

Temp:  C  (F) 

2 

16 

(61.0) 

—  SAMP  It  Ueight  of  Sample 

(grans)  -  -  ~ 

4.18 
3.03 
2.32 
1.47 
2.34 
2.S4 
1.73 
2.12 

1.9S  AV6  =  2.41  grans 
2.38  STN  DV  =  0.76  grans 


116 

117 

118 
118 
120 
121 
122 

123 

124 

125 


Trap  Height  above  Distributor:  cn  (in) 
Bed  Height  above  Distributor:  cn  (in) 
Trap  Height  above  bed  surface:  cn  (in) 
Static  Bed  Height:  cn  (in) 

PI :  cn  Hg  (in  Hg) 

dP:  cn  uater  (in  uater) 

Tenp:  C  (F) 


= 

40.85 

(16.12) 

= 

27.30 

( 10.75) 

= 

13.65 

(5.38) 

= 

22.56 

(8.88) 

= 

7.4 

(2.91) 

= 

4.808 

(1.883) 

= 

16 

(61.0) 

SAMP  #  Weight  of  Sample 


( grans ) 


96 

1.33 

97 

0.81 

98 

1.37 

99 

1.56 

100 

1.42 

101 

1.48 

102 

1.27 

103 

2.05 

104 

1.21 

105 

1.11 

AVG  =  1.36  grams 
STN  DV  =  0.32  grans 


1  AD-A139  111  A  DETERMINATION  OF  PARTICLE  DENS I TV  DISTRIBUTIONS  ABOVE  3/% 
‘  FLUIDIZED  BEDS(U>  MASSACHUSETTS  INST  OF  TECH  CAHBRIDQE 

DEPT  OF  OCEAN  ENGINEERING  G  A  PIPER  MAR  83 
UNCLASSIFIED  N66314-70-A-8B73  F7G  1472  NL 


Trap  Height  above  Distributor:  cm  (in) 

2 

40.35 

(IB. 12) 

Bed  Height  above  Distributor:  cm  (in) 

S 

ZS.B7 

(10.S) 

Trap  Height  above  bed  surface:  cm  (in) 

= 

14. 28 

(S.62) 

Static  Bed  Height:  cm  (in) 

S 

22.  SB 

(8.88) 

Pi:  cm  Hg  (in  Hg) 

2 

6.S 

(2.60) 

dP :  cm  water  (in  uater) 

= 

3.2S3 

(1.283) 

Temp:  C  (F) 

2 

16 

(B1.0) 

,  .  .  .  SAMR  # 

.  Ueight.  of  Sample 

'  -•  ’  .  - 

(grams) 

-126 

0.36 

.....  127 

0.31 

3,21  128 

0.21 

129 

0.27 

130 

0.18 

131 

0.37 

132 

0.34 

133 

0.29 

134 

0.34  AV6  =  0.28 

grams 

135 

0.16  .  STN  DV  =  0.08 

grams 

Trap  Height  above  Distributor:  cm  (in) 
Bed  Height  above  Distributor:  cm  (in) 
Trap  Height  above  bed  surface:  cm  (in) 
Static  Bed  Height:  cm  (in) 

PI:  cm  Hg  (in  Hg) 
dP:  cm  uater  (in  water) 

Temp:  C  (F) 


2 

40  SS 

(16.12) 

2 

26.04 

(10.ZS) 

2 

14.91 

( S. 88) 

• 

ZZ.56 

(8.88) 

2 

6.3 

(2.48) 

= 

2.591 

(1.020) 

= 

16 

(61.0) 

SAMP  # 


Weight  of  Sample 


( grams ) 

106 

0.38 

107 

0.12 

108 

0.15 

103 

0.17 

110 

0.13 

111 

0.18 

112 

0.S6 

113 

0.05 

114 

0.18 

ns 

0.24 

AV6 
STN  OV 


0.ZZ 
0.  IS 


grams 

grams 


Trap  Haight  above  Distributor:  cm  (in) 

=  46.36 

(18.25) 

Bad  Height  above  Distributor:  cm  (in) 

=  28.58 

(11.25) 

Trap  Height  above  bed  surface:  cm  (in) 

=  17.78 

(7.08) 

Static  Bed  Height:  cm  (in) 

=  22.56 

(8.88) 

PI:  cm  Hg  (in  Hg) 

=  8.1 

(3.19) 

dP:  cm  water  (in  water) 

=  6.858 

(2.7) 

Temp:  C  (F) 

=  16 

(61.0) 

-  SAMP#- 

Weight  of 

Sample 

---  ..  .. 

-  (grams)  _  ___ 

—  -  - 

-- 

166 

1.13 

167 

1.18 

.  ..  _ 

168 

1.34 

163 

0.93 

170 

0.88 

171 

1.01 

172 

1.20 

173 

1.15 

174 

1.08 

AVG 

=  1.08 

grams 

175 

0.85 

STN  OV 

-  0.16 

grams 

Trap  Height  above  Distributor:  cm  (in) 

S 

46.36 

(18.25) 

Bed  Height  above  Distributor:  cm  (in) 

= 

27.30 

(10.75) 

Trap  Height  above  bed  surface:  cm  (in) 

2 

19.06 

(7.5) 

Static  Bed  Height:  cm  (in) 

= 

Z2.56 

(8.88) 

PI :  cm  Hg  ( in  Hg) 

= 

7.3 

(2.87) 

dP:  cm  water  (in  water) 

S 

4.808 

(1.893) 

Temp:  C  (F) 

= 

16 

(61.0) 

SAMP  8 

Weight  of  Sample 

(grams) 

156 

0.64 

157 

0.48 

158 

0.78 

159 

0.53 

150 

0.44 

151 

0.58 

162 

0.83 

163 

0.57 

164 

0.52 

AVG  =  0.50 

grams 

165 

0.6Z 

STN  DV  =  0.12 

grams 

Trap  Haight  above  Distributor:  cn  (in) 
Bed  Height  above  Distributor:  cn  (in) 
Trap  Height  above  bed  surface:  cn  (in) 
Static  Bed  Height:  cn  (in) 

PI:  cn  Hg  (in  Hg) 
dP:  cn  water  (in  water) 

Tenp:  C  (F) 


= 

46.36 

(18.Z5) 

s 

Z6.S7 

(10.5) 

= 

19.63 

(7.75) 

= 

ZZ.56 

(8.88) 

= 

6.6 

(Z.60) 

= 

3.Z44 

(1.Z77) 

3 

16.5 

(62.0) 

SAM3  ft 


Weight  of  Sanple 


^ grans) 

136 

-  0.03 

137 

0.13 

138 

0.06 

139 

0.05 

140 

0.07 

141 

0.06 

14Z 

0.05 

143 

0.05 

144 

0.11 

145 

0.14 

AV6 
STM  DV 


a.  08 
a. 04 


grans 

grans 


Trap  Height  above  Distributor:  cn  (in) 
Bed  Height  above  Distributor:  cn  (in) 
Trap  Height  above  bed  surface:  cn  (in) 
Static  Bed  Height:  cn  (in) 

PI:  cn  Hg  (in  Hg) 
dP:  cn  water  (in  water) 

Tenp:  C  (F) 


= 

46.36 

( 18.25) 

S 

Z6.04 

(10. Z5) 

s 

Z0.3Z 

(8.00) 

a 

ZZ.56 

(8.88) 

s 

6.3 

(Z.48) 

= 

Z.631 

(1.036) 

s 

16.5 

(62.0) 

SAMP  ft 


Weight  of  Sanple 


(grans) 

146 

0.01 

147 

0.01 

148 

0.0Z 

149 

0.0Z 

150 

0.03 

151 

0.0Z 

1SZ 

0.01 

153 

0.01 

154 

0.0Z 

155 

0.05 

AVS 
STN  DV 


).  0Z  grans 
).01  grans 


Trap  Height  above  Distributor:  c n  (in) 

= 

59.06 

(23. 2S) 

Bed  Height  above  Distributor:  cn  (in) 

3 

27.94 

(11.00) 

Trap  Height  above  bed  surface:  cn  (in) 

S 

31.12 

(12. 25) 

Static  Bed  Height:  cn  (in) 

3 

22.  SB 

(8.88) 

PI:  cn  Hg  (in  Hg) 

3 

8.2 

(3.23) 

dP :  cn  water  (in  water) 

3 

S.8S8 

(2.7) 

Tenp:  C  (F) 

= 

1S.S 

(62.0) 

— SAMP  8 

Weight  of  Sanple 

(grans) 

•-  .  -  ^ 

176 

0.32  - 

177 

0.22 

•  .  -  ^ 

178 

0.24 

179 

0.S2 

180 

0. 31 

181 

0.42 

182 

0.34 

183 

0.31 

184 

0.23  AV6 

=  0.31 

grans 

18S 

0.21  STN  DM 

=  0.10 

grans 

Trap  Height  above  Distributor:  cn  (in) 

3 

59.06 

(23.25) 

Bed  Height  above  Distributor:  cn  (in) 

3 

27.18 

(10.70) 

Trap  Height  above  bed  surface:  cn  (in) 

3 

31.76 

(12.50) 

Static  Bed  Height:  cn  (in) 

= 

22.56 

(8.88) 

PI:  cn  Hg  ( in  Hg) 

7.3 

(2.87) 

dP:  cn  water  (in  water) 

= 

4.813 

(1.895) 

Tenp:  C  (F) 

= 

16 

(61.0) 

SAMP  8 

Weight  of  Sanple 

( grans ) 

186 

0.25 

187 

0.17 

188 

0.14 

189 

0.17 

190 

0.12 

131 

0. 23 

132 

0.1S 

193 

0.25 

194 

0.18 

AVG  =  0.18 

grans 

195 

0.17 

STN  DV  =  0.04 

grans 

195 


Trap  Haight  above  Distributor:  cm  (in) 
Bed  Height  above  Distributor:  cm  (in) 
Trap  Height  above  bed  surface:  cm  (in) 
Static  Bed  Height:  cm  (in) 

PI :  cm  Hg  (in  Hg) 

dP :  cm  water  (in  water ) 

Temp:  C  (F) 


= 

S3. 06 

(Z3.Z5) 

= 

Z6.67 

(10.50) 

s 

32.39 

(12.75) 

= 

ZZ.S6 

(8.88) 

= 

6.S 

(Z.56) 

= 

3.Z36 

(1.Z74) 

= 

16 

(61.0) 

SAMP  * 


Weight  of  Sample 


__  .J  •  _  .  . 

(grams) 

196 

0.05 

197 

0.07 

138 

0.07 

139 

0.04 

200 

0.06 

201 

0.09 

202 

0.04 

Z03 

0.07 

204 

0.06 

205 

0.06 

AV6 
STM  DV 


0.06 

0.0Z 


grams 

grams 


Trap  Height  above  Distributor:  cm  (in) 
Bed  Height  above  Distributor:  cm  (in) 
Trap  Height  above  bed  surface:  cm  (in) 
Static  Bed  Height:  cm  (in) 

PI :  cm  Hg  ( in  Hg) 
dP:  cm  water  (in  water) 

Temp:  C  (F) 


S 

59.06 

( 23.25) 

S 

26.04 

(10.25) 

s 

33.02 

(13.00) 

s 

22.56 

(8.88) 

s 

6.3 

(2.48) 

s 

Z.619 

(1.031) 

2 

16 

(61.0) 

SAMP  * 

Weight  of 
( grams 

Z06 

0.01 

207 

0.0Z 

208 

0.01 

203 

0.01 

210 

0.01 

211 

0.01 

Z1Z 

0.01 

213 

0.02 

214 

0.01 

215 

0.01 

Sample 

) 


AVG  3  0.01  grams 
STN  DV  =  0.004  grams 


TRAP  HEIGHT  =  1Z.7S  in 


Height  of 

Pressure 

Pressure 

Pressure 

Pressure 

Pressure 

for 

for 

for 

for 

Tap  above 

Samples 

Samples 

Samples 

Samples 

Distributor 

ZS6-Z6S 

Z66-Z7S 

Z76-Z8S 

Z86-Z9S 

( inches) 

(cm  HZO) 

(cm  HZO) 

(cm  HZO) 

(cm  HZO) 

l.S 

S3. 7 

54. 8 

57. Z 

58.8 

Z.5 

48.  S 

49.9 

51.7 

S3.  Z 

3.  S  -  - 

40.8-  - 

4Z.4  -  " 

44.4 

45.8 

4.S 

_  34.9 

35. 7 

38. 7. 

40.4 

S.S  —  - 

Z8.1  -- 

Z9.7 

-3Z.S 

34.0 

6.6 

Z1.6 

Z3.8 

Z6. 1 

Z7.8 

-  7.6  - 

14. T 

17.  Z 

19.1 

Z1.6 

8.6 

8.1 

10.7 

13.5 

16.Z 

9.6 

Z.l 

4.5 

6.8 

9.5 

10.6 

0.0 

0.5 

1.6 

3.7 

11.6 

0.0 

0.0 

0.0 

0.7 

1Z.6 

0.0 

0.0 

0.0 

0.0 

13.6 

0.0 

0.0 

0.0 

0.0 

14.6  -*  ~ 

0.0 

0.0 

0.0 

0.0 

TRAP  HEIGHT  =  14. ZS  in 

Height  of  Pressure  Pressure  Pressure  Pressure 

Pressure  for  for  for  for 

Tap  above  Samples  Samples  Samples  Samples 

Distributor  Z1S-ZZS  ZZS-Z35  Z3G-Z4S  Z46-ZS5 

(inches)  (cm  HZO)  (cm  HZO)  (cm  HZO)  (cm  HZO) 


Trap  Height  above  Distributor :  cm  (in)  =  36. 20  (14.25) 

Bed  Height  above  Distributor:  cm  (in)  =  Z8.58  (11.25) 


Trap  Height  above 

bed  surface:  cm  (in) 

*  7.62 

(3.00) 

Static 

Bed  Height: 

cm  ( Ln ) 

=  22.56 

(8.88) 

Pi:  cm 

Hg  (in  Hg) 

=  8.2 

(3.23) 

dP :  CM 

uater  (in 

uater ) 

=  6.98 

(2.75) 

Temp: 

C  (F) 

=  16.5 

(62.0) 

SAMP  # 

Weight  of  Sanple 

( grams ) 

Z46 

3.34 

247 

5.36 

248 

4.64 

249 

5.71 

250 

4.17 

ZS1 

4.33 

ZSZ 

4.55 

253 

4.33 

254 

*8.44  AV6 

=  4.48  grans 

255 

3.90  STN  DV 

=  0.7Z  grams 

----  ■ 

*  value  not  used  in  coMputing  average. 

Trap  1 

height  above  Distributor:  cm  (in) 

*  36.20 

(14.25) 

Bed  Height  above 

Distributor:  cm  (in) 

=  Z7.30 

(10.75) 

Trap  Height  above 

bed  surface:  cm  (in) 

=  8.90 

(3.50) 

Static 

Bed  Height: 

:  cm  ( in) 

=  22 . 56 

(8.88) 

Pi:  cm 

Hg  ( in  Hg) 

=  7.3 

(2.87) 

dP;  CM 

uater  ( i n 

uater) 

=  4.768 

(1.877) 

Temp: 

C  (F) 

=  16.5 

(62.0) 

SAMP  tt 

Weight  of  Sample 

(grams) 

Z36 

1.55 

237 

1.88 

238 

1.92 

239 

1.80 

240 

1.69 

241 

2.05 

242 

1.84 

243 

1.80 

244 

1 . 87  AVG 

=  1.88 

grans 

245 

2 . 38  STN  DV 

=  0.22 

grans 

198 


Trap  Height  above  Distributor:  cn  (in) 

= 

36.  ZB 

(14.25) 

Bed  Height  above  Distributor:  cn  (in) 

S 

26. 67 

(18.58) 

Trap  Height  above  bed  surface:  cn  (in) 

= 

9.53 

(3.75) 

Static  Bed  Height:  cn  (in) 

s 

22.56 

(8.88) 

Pl=  cn  Hg  (in  Hg) 

= 

6.6 

(2.68) 

dP:  cn  water  (in  uater) 

s 

3.218 

(1.267) 

Tenp:  C  (F) 

= 

16.5 

(62.8) 

SAMP  « 

Weight  of  Sanple 

- - t-  - 

(grans) 

226 

8.93 

227 

8.68 

228 

8.61 

229 

8.56 

238 

8.81 

231 

8.52 

232 

8.38 

233 

8.57 

234 

8.75 

AV6  =0.62  grans 

235 

8.44 

STN  DV  =  8.17  grans 

Trap  Height  above  Distributor:  cn  (in) 

= 

36.28 

(14.25) 

Bed  Height  above  Distributor:  cn  (in) 

= 

26.84 

(18.25) 

Trap  Height  above  bed  surface:  cn  (in) 

S 

18.16 

(4.88) 

Static  Bed  Height:  cn  (in) 

= 

22.56 

(8.88) 

PI :  cn  Hg  ( in  Hg) 

= 

6.3 

(2.48) 

dP:  cn  uater  (in  uater) 

s 

2.624 

(1.833) 

Tenp:  C  (F) 

s 

16.5 

(62.8) 

SAMP  # 

Ueight  of  Sanple 

(grans) 

216 

8.32 

217 

8.35 

218 

8.36 

219 

8.38 

228 

8.32 

221 

8.35 

222 

8.24 

223 

8.38 

224 

8.34 

AVG  =  8.32 

grans 

225 

8.38 

STN  DV  =  8.84 

grans 

Trap  Height  above  Distributor:  cm  (in) 

= 

3Z.39 

(1Z.7S) 

Bed  Height  above  Distributor:  cm  (in) 

S 

Z8.58 

(11. ZS) 

Trap  Height  above  bed  surface:  cm  (in) 

S 

3.81 

(1.50) 

Static  Bed  Height:  cm  (in) 

= 

ZZ.56 

(8.88) 

PI:  cm  Hg  (in  Hg) 

s 

8.1 

(3.19) 

dP:  cm  water  (in  water) 

s 

6.985 

(Z.75) 

Temp:  C  (F) 

= 

17. 

(63.0) 

SAMP  *  Weight  of  Sample 
(grams) 


Z86 

6.55 

Z87 

7.51 

Z88 

7.46 

Z89 

8.05 

Z90 

6.9Z 

Z91 

7.38 

Z9Z 

7.Z3 

Z93 

7. 65 

Z94 

7.04 

Z9S 

6.38 

AV6  *  7.ZZ  grams 
STN  DV  =  0.51  grams 


Trap  Height  above  Distributor:  cm  (in) 
Bed  Height  above  Distributor:  cm  (in) 
Trap  Height  above  bed  surface:  cm  (in) 
Static  Bed  Height:  cm  (in) 
Pl:cmHg(inHg) 
dP:  cm  water  (in  water) 

Temp:  C  (F) 


= 

3Z.39 

( 1Z.75) 

S 

Z7.30 

(10.75) 

S 

5.09 

(Z.00) 

= 

ZZ.56 

(8.88) 

= 

7.3 

( Z . 87 ) 

4.836 

(1.904) 

s 

17 

(63.0) 

SAMP  0  Weight  of  Sample 
( grams ) 


AVG 

STN  DV 


4.35  grams 
1 . 05  grams 


APPENDIX  J 


This  Appendix  contains  the  output  fron  the  inage  analyzer. 
The  output  is  arranged  in  the  following  order. 

1.  Bed  distribution 

Z.  4  cn  freeboard  height 

3.  8  cn  freeboard  height 

4.  1Z  cn  freeboard  height 

5.  "18  cn  freeboard  height 

6.  ZZ  cn  freeboard  height 

7.  31  cn  freeboard  height 


Each  of  the  distributions  at  a  given  bed  height  applies  only 
to  the  Uo/Unf  =  3.81  condition.  For  each  of  the  seven  inage 
analyzer  outputs  listed,  the  following  three  fornats  are  used. 

1.  Histogran  of  absolute  particle  frequency 

vs  particle  dianeter. 

Z.  Cunulative  percentage  plot  of  particle 

size  distribution. 

3.  Table  listing  of  the  above  data. 


The  dashed  line  on  bath  the  bar  graphs  and  on  the  cutulative 
percentage  plots  represent  the  6eussian  distributions  uhich  fit 
the  given  set  of  data.  This  6aussian  distribution  should  be  used 
only  as  a  rough  estimate  of  the  data  because  the  analysis  included 
“particles'*  less  than  79  microns.  Sy  viewing  a  blank  slide  with 
only  the  tape  applied,  these  "particles”  were  confirmed  to  be 
bubbles  and  dirt  entrapped  in  the  adhesive  on  the  tape  used  to 
hold  the  sample  particles  in  place. 


Imq>  Anilyitr  Output 
Bod  Portlclo  Distribution 


CLASSIFICST 

:on  list  for 

©CIRCLE  IN 

CHANNELS 

1 

UNDERFLOW 

8  OVERFLOW  8 

- 

‘ 

CLASS' 

FROfl 

TO 

FREQUENCIES 

ASS 

REL  CUN , 

,  ass 

CUN.  REL 

1 

.  48 9 77 E- 01 

.14344 

13. 

4.89  2 

13. 

4.89  2 

2 

. 14344 

.24191 

24. 

7.82  Z 

39. 

12.79  2 

3 

.  2419.1 - 

.4^837 

12. 

3.91  Z 

51. 

16.61  Z. 

4  ■ 

■  m837 

.43484 

6. 

1 .95  5; 

57. 

"18.57  2 

* 

w 

.4-3484 

. 33 1 38  ~ 

18.  i  “ 

3.26  Si  _ 

67. 

21.82  2 

6 

.33138  * 

.62777 

17. 

5.34  2' 

84. 

27 . 36  2 

7 

.62777 

.72423 

28. 

6.31  Si 

184. 

33 . 88  2 

.  a  . 

:.7242S 

.82878  ... 

23.  _ 

7.49  Si 

127. 

41.37  2 

a  • 

.82878  . 

.91716 

23  ■ 

*7.49  Si-'— 

158. 

48 • 86  2 

18  . 

.91716 

1.8136 

28. 

6.51  Z 

179. 

55.37  2 

:  -11  - 

1.8136- 

1.1181  : 

14.-  ~ 

4.56  Si 

184. 

39.93  2 

12 

1.1181 

1.2966 

23. 

7.49  Z 

287. 

67.43  2 

is  " 

1.2866 

1.3938 

27. 

8.79  2 

234. 

76  •  22  2  - 

14 

1.3838 

1.3993 

12. 

3.91  Z 

246. 

88.13  2 

13 

1.3993 

1.4968 

17. 

5.34  Z 

263. 

85 . 67  2 . 

16 

1.4968 

1.3924 

7.. 

2.28  5 

278. 

87.95  2 

17 

1.3924 

1.6889 

7. 

2.28  St 

277. 

98.23  2 

.18--- 

1.6889 

1.7833 

11.  . 

3.58  Si 

288.  . 

93.81  2 

19 

1.7833 

1.8818 

4. 

1.38  2~~' 

292. 

93.11  2 

28 

1.8818 

1.9783 

1. 

.  .33  Z 

293. 

93.44  2 

21 

1.9783 

2.8747 

r. 

.33  Si 

294. 

95.77  2 

22 

2.8747 

2.1712 

5. 

1.63  2 

299. 

97.39  2 

23 

2. 1712 

2.2677 

i. 

.33  2 

388. 

97.72  2 

24 

2.2677 

2.3641 

9. 

8.88  2 

308. 

97.72  2 

23  ■ 

2.3641 

2.4686 

8. 

8.89  2 

308. 

97.72  2 

26 

2.4686 

2.3371 

1 . 

.44  2 

381  . 

98.85  2 

17 

2.5371 

2.6335 

1 . 

.33  2 

382. 

98. 4«  2 

28 

2.6335 

2.7388 

8. 

8.89  2 

302. 

98.37  2 

29 

2.7388 

2.8463 

9. 

9.88  2 

382. 

98  •  4*1  2 

38 

2.8465 

2.9429 

1 . 

.33  2 

383. 

98.79  2 

31 

2.9429  . 

3 . 9394 

9. 

8.88  2 

383. 

98.78  2 

32 

3 . 8394 

3. 1359 

1 . 

.44  2 

304. 

99.82  2 

3.1339 

3.2323 

9. 

8 . 88  2  . 

384. 

99.82  2 

3^- 

3.2323 

3.3288 

9. 

8.89  2 

304. 

99.02  2 

33 

3 . 3288 

3.4233  . 

8. 

8.88  2 

394. 

99.82  2 

36 

3 . 4253 

3.3217 

9. 

9.88  2 

304. 

99.02  2 

V  » 

3.3217 

4.6 182 

9. 

8.88  2 

384. 

99.82  2 

38 

3.6182 

3.7147 

9. 

8.99  2 

304. 

99.02  2 

39 

3.71 47 

3.8111 

1 . 

.  44  4 

383. 

99.33  2 

48 

3.811* 

3.9876 

2. 

.65  2 

307 . 

198.99  2 

v  r.  ?.-'i  .v.*. 


TT^TTTWTrrin-r^n 


Imq*  ftnelyzor  Output 
Height  Distribution 
4  en  Above  Bed 


curesif:-" 

u«ocan_otr 

CLASS 


1 

•y 


■s 
■  4 

"  r  :. 

7  ... 

8 

^  9  - 

10  ?. 

11 

12 
15 

14  . 

13 

-IB- 

17  : 

18 


f: 


I 


:-3H  LIST  FOR  OCIRCLS  IM  CHBHHELS 

8  ouerflow  8 


FRcn 


. 12788 
.21427 
• 08886 
.38743 
.47484 
.36864 
.8*723 
.73382 
.82841 
.88788 
. 93339 
1.8882 
1 . 1668 
1.2334 
1^3488 
1 . 4263 
1.3131 
1.3997 


TO 


.21427  .. 

. 06686 
.38743 
. 47484 
; 36864  . 

. 64723  - 
; 73382  — 
.82841  : 

-.98788 

".99339  *•. 

1 . 8882 
1 . 1668 
1.2334  .  - 
1 . 5488  ' 

1.4263 
1 .3131 
1.3997  - 
1.6863 


BBS 

9. 

.11 . 
3. 
11 . 
34.-  ■ 
54  r 

45. 

36 . 
27, 
25. 
18.* 
13. 
9. 
6. 
3. 
..  7. 
**  T. 

1 . 


FREOUEMCIES 
REL  _  cun.  BBS 


cun.  REL 


2.78 
3.48 
.93 
3.48 
18.49 
16.67 
13.89 
11.11 
8 . 33 
7.72 
*  5.56 
4.65 
2.78 
1.83 
1.54  X 
2.16  Z 
.  .31  Z 

.31  X 


X 

X 

X 

Z 

X 

z 

X 

X 

X 

be-! 

«• 

X 

Z' 

X 

z 


9. 
28. 
23. 
34. 
68. 
122. 
167.^- 
283. 
238. ... 
255. 
273. 
288. 
297. 
383 . 
088  • 

315. 

316. 
317^ 


2.78  X 

6.17  X 

7.18 
18.49 

28.99 

37.65 
_ ,51 .54 

62.65 

78.99 
78.78 
84.26 

;  88.89 
91.67  2 
93.32  2 
95.86  2 
97.22  2 
97.53  2 
97.84  2 


19 

1 . 6863 

- -1.7729  - 

2. 

.62  Z 

319. 

98.46  X 

28 

1.7729 

1.8395 

2. 

.62  X 

321 . 

99.87  X 

21 

1.3595 

1.9461 

1 . 

.31  X 

322. 

99.38  X 

1.9461 

2.8327 

8. 

8.88  2 

■04>4>  * . 

99.38  * 

23 

2.8327 

2.1193 

8. 

8.88  X 

^22 , 

99.38  X 

24 

2. 1193 

2.2859 

8. 

8.68  X 

799 

0^4.  1 

99.38  X 

23 

2.2839 

2.2925 

1. 

.31  X 

T9T 

99.69  X 

26 

2.2925 

2.3798 

8. 

8.88  X 

o2o . 

99.69  X 

27 

2.3798 

2.4656 

8. 

8.88  X 

•  323 . 

99.69  X 

28 

2.4656 

2.3322 

8. 

8 . 88  X 

323. 

99.69  X 

29 

2.3322 

2.6388 

8. 

8.88  X 

o2o . 

99.69  X 

36 

2.6388 

2.7234 

8. 

8.88  X 

323. 

99.69  X 

31 

2. 7254 

2.8128 

8. 

8.88  X 

o2o . 

99.69  X 

W4h» 

2.8128 

2.8986 

8. 

8.88  2 

323. 

99.69  X 

WW 

2.8986 

2.9832 

8. 

8.88  X  . 

o2o . 

99.69  X 

2.9832 

3.8718 

8. 

8.88  Z 

323. 

99.69  X 

T« 

-  3.8718 

3. 1584 

8. 

3.88  X 

323. 

99.69  X 

35 

3. 1384 

3.2449 

8.  . 

8 . 88  X 

323 . 

99.69  X 

37 

3.2449 

3. 331 5 

8. 

8.88  X 

323. 

99.69  X 

38 

3.3315 

3.4181 

8. 

9.99  X 

323  • 

99.69  X 

39 

7.4181 

3. 5847 

8. 

9.99  X 

323 . 

99.69  X 

*8 

•  T*  — 

«■  -• 

3.5913 

1 . 

.31  X 

324. 

188.88  X 

207 


,  *.»  •  .  v  _  • '  **  v"  •“  »*  ( 


,  >■.  «■. 


■  *  -s  v 


’  «  *»  *  >  *  .  1  •  *  .  * 


CUM*  f  9.E4UENCY 

190. 9- 

90.00- 

- V 

O’ 

• 

40.00- 

y 

29.00- 

/ 

0.000- 

/! 

asssftltj  s 

CUflSSSI 

^?TK'».3S 

. . . . 2. 91<S2l3 

st. 6ev. 

.  .  .,.8.522135 
MEDiftri 
...J.7«  1852 

•  •  •  •  *S?«  'J  -ww  v  ■•  1 


0C2RCLS 


Haight  Distribution 
t  cn  Abovo  Bob 


classificatic:*  .:r  for  dcircle  in  channels  i 


UNOESPwOH 

8 

OVERFLOW  8 

CLASS 

nton 

TO 

FREQUENCIES 

ASS 

REL 

CUN.  RBS 

cun.  REL 

1 

.14867 

. 28224 

2  ' 

.63  2 

2. 

.63  2 

-* 

4b 

.28224 

.25581 

2. 

.63  X 

4.  . 

1.26  2 

.  w 

.25581 

.38938 

3. 

.95  X 

7 . 

2.21  2 

.38938 

..  =  36294 

.  0. 

8.08  X 

7. 

2.21  2 

=■5!  3 

.36294 

=-57.41631 

-.  4. 

-  1.26  X 

11. 

3.47  2 

6 

.41651 

sr*-.  47888 

•  9- 

2.84  X 

20. 

6.31  2 

— :  r  ... 

-.47889 

'  -7  . 52365  .  7 

--28. 

r-~8.83  X  - 

48. 

15.14  2 

8  . 

.52363 

.57722 

21 . 

'  6.62  X 

69. 

21.77  2 

— .  .  8 

... 57722 

. 63878 

27. 

8.52  2 

96. 

39.28  2 

18 

.63878 

-7. 2768435 

’  32. 

"  10.89  X  _ 

128. 

44.38  2 

11 

. 68435 

—  73792 

21 . 

-  6.62  X 

149. 

47.08  2 

12 

.73792 

.79149 

22. 

6.94  X 

171. 

53.94  2 

13 

.79149 

.84586 

38. 

9.46  2 

281. 

63.41  2 

-  14 

.84586 

.89863 

17. 

3.36  2 

218. 

68.77.2. 

15 

.89863 

.95219 

19. 

3.99  2 

237. 

76.76  2 

16 

.95219 

1.8858 

11 . 

3.47  2 

248. 

78.23  2 

17 

1.0858 

7  ~ 1 .8593 

—  11. 

---  3.47  2  . 

259. 

81.70  2 

18 

1 . 0393 

1.1129 

9. 

2.84  2 

268. 

86.54  2 

15 

1.1129 

—  1 . 1665 

8. 

2.32  2 

276. 

87.87  2 

28 

1 . 1663 

1.2288 

18. 

3.13  2 

286. 

90.22  2 

21 

1.2288 

1.2736 

2. 

.63  2 

238. 

90.83  2 

1.2736 

1 .3272 

3. 

1.38  2 

293. 

92.43  2 

23 

1.3272 

1 . 3887 

7. 

2.21  2 

309. 

96.64  2 

24 

1 . 3887 

1.4343 

z . 

.95  2 

• 

95.38  2 

25 

1 . 434 3 

1.4879 

1 . 

.32  2 

304, 

95.90  2 

26 

1 . 4879 

1.3414 

4  « 

.63  2 

386. 

96.53  2 

27 

1.5414 

1 .3958 

1  . 

.32  2 

j0.'  . 

96.35  2 

28 

1 . 5938 

1.6486 

1  . 

.32  2 

308. 

97.16  2 

29 

1.6486 

1.7821 

1  . 

.32  2 

309. 

97.48  2 

38 

1 . 7821 

1 . 7357 

2. 

.  63  2 

311 . 

98.11  2 

31 

1.7557 

1 . 8893 

0. 

0.00  2 

311 . 

98.11  2 

32 

1 . 8893 

1.8628 

4< 

.63  2 

313. 

98.76  2 

wO 

1 . 8628 

1.9164 

4 

1  ( 

.32  2  \ 

314. 

99.05  2 

34 

1 .9164 

1 . 9708 

1  . 

.32  2 

313. 

99.37  2 

35 

1.9798 

2.0236 

9. 

9.08  2 

315. 

99.37  2 

36 

2.8236 

2.0771 

1 . 

.32  2 

316. 

99.68  2 

** 

2.0771 

2.1387 

8. 

9.08  2 

316 . 

99.68  2 

38 

2. 1307 

2. 1843 

8. 

8.09  2 

316 . 

99.68  2 

39 

2. 1843 

2.2378 

0. 

0.89  2 

316 . 

99.53  2 

48 

*  “t*: 

2.2914 

1  . 

.32  2 

sir. 

180.80  2 

ft3<5  .  F?.E4ljeKCY 
110.0-! 


1kooo  Anaiimr  Output 
Haight  Distribution 
12  cn  Above  Bod 


55.00 


8.000 


2.  888 


cuw.faeniSHCY  x 

108.0-1 


80.00-  /■ 

7/ 

00.00-  j; 

i: 

40.08-  j 

•* 

:\ 

28.88-  :l 

I f 

0. 800-k - 


2.088 


mihz  '  lit 

khhb  I 

CLSSSSl  48 

XN7SRV8L 

.  .  . .8. 12S249 

niriif-iuM 

.  .  .  .  0 . o 2 1? J0C-0 1 

MfiXittUM 

.....5.11212 

It. 6ev. 

•  *1  IT  » 


QCXRCLE 


S.  000 


.SS'JS5Scs  ii2 

8oi*ftSK  ,| 

CLhSScs  40 

INTE|Vl^''4« 

MiNinur? _ 

.  .  .  .8  .  o21  ✓  36c.-8  l 


•  I  •  1  V  •  V  *  4  f 

^7^11212 
"J?V  "19992 

lt.o|v.  ..... 

.  .  .  .8 . 

MEOi0M _ _ 

.  .  .  .0.6554.0 
mooe 

•  m  •  <  •  9  A  w  4  •  * 


9.000 


1.00 . 0T 


39. 9<H 


S0.001 


COUNTS  _ 
in  RBMCE 
IJHD.  FLOW 

CLBSSca 


ZZ9 

329 


IMTERyflU  „  a< 
.  .  .  .0  .  <i»l/  3«.c.-0  1 
MIHiMUM 
. . . . 9 . I 10231 

r**  ?r«?™ 

351330 


Inogc 

1  Anolyzor  Output 

Haight  Distribution 

18 

ch  Abovo  Bod 

CLASSIFICATION  LIST  FOR  DCIRCLE 

IN  CHANNELS 

-.1! 

-  t 

T 

UNDERFLOW 

0 

OVERFLOW 

0 

- 

■* 

CLASS 

FRQtl 

TO 

FREQUENCIES 

'  ' 

ASS 

REL 

CUN.  ABS 

CUN.  REL 

4 

■ 

.11028 

«  17445 

1 . 

.31  Z 

1  . 

.31  a 

2 

. 1744S 

. 23864 

8. 

2.50  Z 

9. 

2.31  a 

Z 

.23854 

. 38282 

1 . 

.31  Z 

18. 

3.12  a 

4 

.30282 

.35788 

5. 

1.56  Z 

15. 

4.59  a 

* 

.36780 

.43118 

4. 

1.23  a 

19. 

5.94  a 

5 

.43118 

•  49-J3S 

17. 

5.31  Z 

36 . 

11.25  a 

-  .  .  .7 

. 49536 

.55954 

16. 

5.00  Z 

52. 

16.25  a 

*  *  -  8 

.55954 

62372 

26 . 

8.12  Z 

73. 

24.37  a 

-  9 

.62372 

%  68789 

32. 

10.80  a 

118. 

34 . 37  a 

10 

.58789 

.75287 

31. 

9.69  a 

141 . 

44.06  a 

11  . 

. 35207 

.81625 

35. 

18.94  a  : 

176. 

55.38  a 

12 

.81625 

.88043 

21 .  _ 

5.56  a 

197. 

61.56  a 

13 

.88843 

.94461 

26. 

8.12  a 

223. 

69.-69  a 

14 

.94461 

1.0888 

20. 

6.25  a 

243. 

75.94  a 

15 

1 . 8888 

1.0730 

13. 

4.86  a 

256. 

30.08  a 

IS 

1 . 8  f  O0 

1.1371 

13. 

4.06  a 

269. 

84.06  a 

17 

1 . 1371 

1.2013 

7. 

2.19  a 

275. 

86.25  a 

19 

1.2013 

1.2655 

6. 

i.87  a 

232. 

88.12  a 

- 19 

1 . 2655- 

1.3297 

7. 

2.19  a 

239. 

98.31  a 

20 

1.3297 

173939 

2. 

"  .63  a 

291. 

90.94  a 

21 

1.3939 

1.4588 

4. 

1.23  a  . 

295. 

92.19  a 

22 

1 . 4588 

1.5222 

7. 

2.19  a 

382. 

94.37  a 

23 

1.5222 

1 . 5864 

4. 

1.25  a 

306. 

95.62  a 

24 

1 . 5864 

1.6506 

4. 

1.25  a 

310. 

96.87  a 

23 

1 .6586 

1.7148 

1  . 

.31  a 

311 . 

97. 19  a 

25 

1.7148 

1.7789 

£>, 

1.2s  a 

315. 

98.44  a 

27 

1 . 7789 

1.8431 

.63  a 

317. 

99.06  a 

29 

1.8431 

1 .90f  o 

0. 

0.88  a 

317. 

99.06  a 

29 

1 . 9073 

1.9715 

8. 

9.08  a 

317. 

99.06  a 

30 

1 .9715 

2.0357 

0. 

0.00  a 

317 . 

99.06  a 

?  4 

V  t 

2.8357 

2.0998 

0. 

8.88  a 

317. 

99.06  a 

VI  ^ 

2.0998 

2. 1640 

1  . 

.31  a 

318. 

99. Of  a 

vw 

2. 1648 

2.2282 

8. 

0.88  a 

313. 

99.37  a 

3^ 

2.2282 

2.2924 

8. 

8.00  a 

313. 

99.37  a 

^  C 

w 

2.2924 

2.3566 

8. 

8.08  a 

313. 

99.37  a 

38 

2 . 3565 

2.4207 

0. 

8.00  a 

313. 

99.37  a 

37 

2.4297 

2.4849 

8. 

0.00  a 

313. 

99.37  a 

38 

2.4045 

2.5451 

1  . 

.31  a 

319. 

99.89  a 

39 

2.5491 

-9  C  0  7? 

8. 

8.80  a 

315. 

99.69  a 

40 

2  5  133 

2 . 5774 

1  . 

.31  a 

“20. 

100.00  a 

A.3S*F?.EAIJ£NCT 
30.09 


40.00 


30.00 


20.88 


10.88-4 


8 . 880 


Inago  Analyzer  Output 
Haight  Distribution 
2Z  cm  Abovo  Bod 


*-*«*■** 


COUNTS  3t-l 

aN  RANGE  3i| 

UNO. FLOW  0 

OVERFLOW  8 

CLHSSts  -  .  40 

INTERVAL 

.  .  .  . 8.340214E-01 

MINIMUM 

. « .*0. 1 10^7  ^ 

mbximum 

.  • « . .2. 4309^ 

MEAN 

...  .0.3*9239 
St.O&V. 

....  0. 395442 

MEDIAN 

....  0 . 33^  s 8  + 

MODE 

_ 0 . *32213 


1 . 000 


l.  300 


3.  000 


4. 


□CIRCLE 


CUM.F  R.EG.UEN.CT  I; 
100.0- 


89.004 


50 . 00-1 


48 . 88-f 


>8.884 


8.888-4 


COUNTS 
.IN  RANGE 


/ 


/  * 

// 


/ 


•ill 

.0 

0 

40 


INTERVAL 

7. . . 0. 348214E-01 

MINIMUM 

....0.1200*5 

T?^T43093 

?T?%.  3*9239 
ST.OEV. 

....  8  .  -195444. 

MEDIAN  _ 

•  .  .  .0. 33/ f 34 
MOOE 

A  7  *  -4  •?  1  ? 

•  •  •  •  V  •  I  4  In  &  1  -4 


1 . 000 


000 


T.  800 


-  .  000 


W  *-r  _  7.  W  W  — 


Z14 


Imqo  Analyze  Output 
Haight  Distribution 
22  cn  Abovo  Bod 


CLftssirrcftTiofi  wisr  res.  dcircle  in  channels 


UNDERFLOW 


OVERFLOW 


kSS 
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9. 

9. 
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3. 

3. 

3. 

3. 

3. 

3. 
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73.83 

5.53 
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2  ‘ 
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3.47 

Z 

272.'  '  . 
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Z 

1.38  2 
3.13  Z 
.93  2 
.32  2 
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1.26  2 
.32  2 
.93  2 
.32  2 
0.88  2 

3.88  2 

8.88  2 
0.00  2 

.32  2 
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3.38  2 
0.33  2 
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2 
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314. 
313. 
316 . 
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APPENDIX  K 


This  Appendix  contains  the  photographs  of  the  oscilloscope 
traces  obtained  uhile  sampling.  Each  of  the  pictures  is  labeled 
uith  the  sample  number  for  which  it  represents. 
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APPENDI X  L_ 


This  Appendix  contains  a  complete  listino  of  the  particle 
trajectory  computer  model  used  in  this  research.  In  addition,  an 
in  depth  program  description  and  flou  diagram  is  included. 


This  section  describes  the  computer  progran  listed  in 
Appendix  L-  To  assist  in  understanding  the  progran  logic,  the 
progran  itself  contains  nunerous  connent  statements.  The  progran 
consists  of  a  nain  progran  and  seven  (7)  subroutines.  The  nain 
progran  controls  the  input  of  paraneters,  trajectory  calculations 
and  output  selections.  The  subroutines  control  the  actual  data 
output  in  either  table  forn  or  graphics.  The  follouing  is  an  in- 
depth  description  of  the  progran: 


Lines  Description 


50-150  This  section  is  used  to  explicitly  define 

the  najor  variables  used  within  the  progran. 
The  large  arrays  are  defined  in  connon  block 
forn  to  save  nenory.  Constants  used  in  the 
progran  are  also  defined. 

190-610  This  section  is  used  to  input  variable  data 

to  the  progran.  Three  options  exist  for  the 
bed  particle  distribution  input: 

1)  The  default  condition  sets  the  quantity 
of  each  particle  to  unity.  This  option  is 
generally  used  uhen  a  height  determination 
is  required  or  being  sought  after. 

Z)  The  particle  number  option  allows  the 
entry  of  bed  distribution  by  the  number  of 
particles  present  in  each  diameter  range. 
This  is  used  when  data  from  the  image 
analyzer  was  being  used. 

3)  The  last  option  is  used  uhen  the  particle 
size  distribution  is  determined  from  a  sieve 
analysis  and  the  data  is  measured  in  grams 
mass.  The  program  will  then  determine  the 


particle  number  density  based  on  the 
assumption  that  the  particle  are  spherical. 


650-690  These  lines  convert  the  velocity  input 

values  from  Ft/s  to  m/s . 

700-710  Determine  the  slope  of  the  triangular  jet 

using  the  input  amplitude  and  ciiration. 

720-800  Initialize  the  data  arrays  to  zero  (0>. 

810-1660  This  section  calculates  the  trajectory  for 

each  diameter  particle.  This  UHILE  condition 
contains  the  following  8  subsections. 

820-840  Initialize  the  height  of  the  particle  to 

pass  the  first  UHILE  statement.  Calculate 
the  particle  diameter  to  be  used. 

850-1340  Calculate  the  particles  trajectory 

parameters  while  the  particle  is  above  the 
bed  surface.  This  00  loop  contains  the 
following  6  subsections. 

860-940 . If  the  elapsed  time  since  the  particle  left 

i  the  bed  is  less'than  the  jet  duration  tine, 

the  add  the  jet  velocity  to  Uo. 

950-990  Calculates  the  relative  velocity  of  the  air 
with  respect  to  the  particle.  Determine  the 
sign  of  the  drag  force. 

1010-1060  Calculate  Reynold's  number. 

1100-1170  Calculate  the  drag  on  the  particle  and  then 
determine  the  particles  new  velocity  and 
position. 

1180-1220  If  the  particles  velocity  is  positive, 

record  the  new  max  height  and  elapsed  tine. 

1230-1320.  Add  1  count  to  the  probability  array  in  the 
storage  position  representing  the  particles 
height/2.  By  using  the  DIV  statement,  a 
height  window  (dH>  of  2  cm  is  created. 

1350-1370  Change  all  heights  from  Ft/s  to  m/s 

1380-1430  Multiply  probability  distribution  by  bed 
distribution  weighting  factor. 

1440-1500  Determine  the  maximum  value  and  normalize 
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the  data  such  that  the  maximum  is  100. 

1510-1560  Calculate  the  particle  entrainment  by 

sunning  the  total  volume  of  each  particle 
distribution  at  a  given  height. 

1570-1660  Normalize  the  entrainment  data  to  100  and 

calculate  the  In  value  if  not  equal  to  zero. 

1700-1820  01 splay  output  menu. 

1830-2210  —  Controls  the  selected  output. 

2220  End  of  MAIN  program 

2260-2520  Sub  0isplay_0ata  is  used  to  list  input 
parameters  used  and  the  maximum  height 
attained  by  each  diameter  particle. 

2560-3930  Sub  Oisplay_graphics  controls  the  processing 

and  output  of  graphic  information.  This 
subprogram  contains  the  following  4  sections 

2640-3120  Performs  scaling  for  determining  graphic 
dimensions  and  limits. 

3130-3240  ..  Creates  graphic  display  and  labels  the  X  and 
Y  axis. 

3250-3830  Controls  the  plotting  of  selected  graphic 
output. 

3840-3920  Prints  hard  copy. 

3980-4130  Function  routine  for  determining  scaling 

factor  for  graphic  display. 

4180-4280  Finds  the  limit  of  data  uithin  an  INTE6ER 
array  to  limit  the  X  axis  on  a  plot  to  the 
range  in  which  the  Y  values  are  non-zero. 

4330-4460  Determines  the  maximum  value  stored  in  an 
INTE6ER  array. 

4510-4640  Same  as  4330-4460  except  for  REAL  array. 

4690-4770  Same  as  4180-4280  except  for  REAL  array. 


ENO  of  PROGRAM 


RAIN  Progran 


OPTION  BASE  1 

INTEGER  Answer, Soecify, Point ,Uiew,8ed_nin, I 
REAL  Uo.Upo.Uj ,Dei_t,Part_height .Slope, Jet_t, Jt 

COR  INTEGER  Oistr ibut ion< SB, 80  ) ,  Qistrib_dens  i"ty< SB , SB ) .REAL  Height! 2, 58), 3 


*d_p«rt ic is( SB ) .Entrapment <2,88) 
88  OIH  HoneSC21,ClearSC2) 

188  C 1  e#rs«C}tf?S{  2S5  )&CHRS(  7S ) 

118  Hone*-CHR*< 2SS )8CHR*< 84 ) 

128  ViscosityJcin-l.*86E-5 
138  Oensity_air-i .281 
1*8  Oenait y_part-7 . 3SE+3 

1S8  6ravity*9.806SS 


I  CLEAR  SCR  key 
!  HOME  key 
!  air,  n*n/s 
I  kg/n»n*n 
I  kg/n»n»n 
!  n/s*s 


I  Input  varied las 


!  Bad  dist  wanted 


198  PRINT®  IS  1  I  Output  to  CRT 

280  INPUT  *£nt«r  naan  air  velocity  Uo  (Ft/*):  *,Uo 

218  INPUT  "Enter  anplituda  of  Jet  velocity  (Ft/s):  *,UJ 

220  INPUT  "Enter  initial  particle  velocity  Upo  (Ft/s):  *,Upol 

230  INPUT  "Enter  duration  of  jet  (a):  *,Jet_t 

Z40  INPUT  "Enter  tine  increnent  for  iteration  (a):  ",0el_t 

ZS8  INPUT  "Do  you  desire  to  input  the  bed  particle  distribution  (1  yes, 8  no):- 
.Answer 

ZS0  “  IF  Answer— 0  THEN  — =^.=. - *  ^ — ■=  |  gad  dist  not  wented 

270  FOR  I- l  TO  SB  I  Set  bed  dist  to  1 

280  Bed_particle< I )«l. 

290  NEXT  I 

308  ELSE 

310  IF  Answer*!  THEN  !  Bed  dist  wanted 

320  INPUT  "Input  will  be;  weight  in  grans  ( I )  or  nunder  (2):  ",  Specify 

330  IF  Specify*!  OR  Specify*2  THEN 

5A0  PRINT  "Enter  data  for  each  dianeter*  !  Enter  gran  or  * 

350  Min_b*d“l 

5S0  FOR  1*1  TO  S0  !  Each  dianeter 

370  PRINT  "0-  " : I • 10+70; *  un" 

380  INPUT  8ed.j>article<  I ) 

390  IF  Sed.jjerticl*<  I  ><Bed_pertiele(Min  bed)  THEN  Min  bed-I 

*08  NEXT  I 

*10  IF  Specify*l  THEN  !  If  weight  entry 

*20  FOR  1*1  TO  SB  !  Calculate  t 

‘30  8sd_pertiele< I >*5»8ed_particle< I  )/<PI*< 1*10+70 )"3*0ensity_part  > 

**0  IF  8ed_particie< I  )<Bed_particie(Min  bed)  THEN  Min  bed” I 

tsa  next  i 

*60  ENO  IF 


GOTO  22 3 
ENO  IF 
ELSE 

SOTO  2=3 
ENO  IF 
ENO  IF 
Ha*  Jaed«l 

FOR  I"!  TO  S3  i  Fmd  largest  t 

IF  Beo_particis< I '>8ed_par‘ ;c  l«<  Maxjied  )  THEN  MaxJjed-i 


!  determine  scaler 
!  Scale  data  max-108 


S70 

NEXT  I 

S80 

flax.value*  100/Bed_part  ic  le<  Max  bed  > 

!  determine  scaler 

530 

FOR  I-l  TO  50 

!  Scale  data  max-108 

S00 

Sedjarticle! I >-8ed.part icle! I >«Hax.value 

S10 

NEXT  I 

520 

! 

S30 

!  Compute  data 

540 

1 

550 

Ul-Uo  ‘ 

550 

Ugl-Upol 

570 

Uo-'Jo«  .3048 

!  Change  to  m/s 

580 

Upol-Upol* .3048 

530 

Ujl-Uj*.3048 

700 

IF  Jet_t-0  THEN  Jet.t-1 

710 

Slo0e-2.»UJ1/Jet_t 

!  Triangle  jet  pulse 

720 

FOR  I-l  TO  50 

730 

Height <1, I >-8. 

!  Zero  height  array 

740 

Height (2, I )-0. 

*  Zero  nunoer  array 

750 

FOR  J-l  TO  80 

750 

Distribution! I, J >-8 

!  Zero  dist  array 

770 

Entrainment! 1 ,J)-0. 

!  Zero  entrain  array 

780 

Entr ainment ! 2 , J >-8 . 

790 

NEXT  J 

800 

NEXT  I 

810 

FOR  I-l  TO  SO 

1  For  each  particle 

820 

Part.height - . 000001 

!  Set  for  WHILE  state 

830 

Uoo-Uool 

1  Initial  U  constant 

840 

Diameter  __part-<  1+7  >•! .0E-5 

•  Particle  Oie  un 

850 

WHILE  Part_helght>0 

860 

Jt-Jet_t/2. 

870 

Tine_nou»Height!2,I >»0el_t 

!  Time  of  part  flight 

880 

IF  Tine_nou<Jet.t  THEN 

!  part  in  jet? 

890 

IF  Height(2.I)*0el.t<Jt  THEN 

!  1st  half  of  jet 

900 

Uo-Uo+Slope*Haight< 2 , I >*Qel_t 

910 

ELSE 

'  * 

920 

Uo-Uo+Ujl-»Slope*Height<2,I  >+0el_t 

!  2nd  half  of  jet 

930 

ENO  IF 

940 

ENO  IF 

350 

Relative.vcl-Uo-Upo 

!  Rel  V  seen  by  part 

950 

Sign-1 

!  Set  sign  positive 

970 

IF  Relative_vel<0  THEN 

i  Is  U  negative? 

380 

Sign— 1 

!  Set  sign  negative 

930 

ENO  IF 

1000 

1 

Calculate  Re 

1010 

Reynolds.no-Relat ive_vel*Oiameter_part/Viscosity_kin 

1020 

Rayrtolds.no-ABS!  Reynolds.no ) 

1030 

IF  Reynolds .no-0.  THEN 

1040 

Acceleration— eravity 

1050 

SOTO  1150 

1060 

ENO  IF 

1070 

l 

Calculate  drag 

1080 

!  Uses  data  correlation  for  drag  coefficient 

good  for  Re<l£5 

1090 

! 

1100 

Oran  coa f  f -24/Reyno 1 ds_no+S / ( 1 +SOR ( Reyno 1 ds. 

no  >  >+.4 

1110 

Or sg-S ign«Or ag_eae f  f  *Oens i t y.e ir »Re 1 at i ve.ve 

i»Reiat i ve.ve 1*P I»0 tame 

_part 

•Oianeter_part/8 

1120 

i 

Calculate  acceierat 

1120 

Acceierat  ion-Orag/(  Oens  i  tyjart  »PI •Qiamet sr_part  *3/5  )-Sravi  ty 

1140 

1 

Calculate  oel  veloc 

1150 

Veloci ty*Uoo*Accal«rat ion*0«l_t 

1158 
1 178 
1180 
use 
1208 
1210 
1220 
1230 
1240 
1253 
1260 
1270 
1230 
1230 
1300 
1310 
1320 
1330 
1340 
1350 
1360 
1370 
1380 
1330 
1400 
1410 
1420 
1430 
1440 
14S0 
1460 
1470 
1480 
1430 
1500 
1510 
1520 
1530 
1S40 
1550 
1560 
1S70 
1530 
1530 
1600 
1610 
1620 
1620 
1640 
16=0 
1560 
1673 
1630 
1630 
1700 
1710 
1  f  23 
1730 
1740 
1750 


Calculate  new  position 

_t/2 

•  Pert  still  rising? 

!  Save  position 
!  Inc  t  for  nex  rise 

!  Set  new  U  for  next  inc 
!  2  en  wide  storage  bins 
!  Set  default  for  fatal 


!  Set  default  for  fatal 


!  Indicate  conp  working 
!  Save  t  tines  pert  in  height  bin 

Distribution! I .Point* l )-0i  str ibut ion( I , Point +1 >+l 
£N0  UHIL5 
NEXT  I 

FOR  1-1  TO  50 

Height! 1 , I )-Height! l, I >-100  !  Chang*  to  cn 

NEXT  I 

FOR  1-1  TO  50 

Ualue-9ed_particle! I )  !  Part  size  weight  factor  fron  bed  dist 

FOR  J-l  TO  80 

Oistrib_density< I , J )-Oistribution! I , J )>Oalue  !  Weight  dist  values 
NEXT  J 
NEXT  I 

N«x_dist»FNf1ax_intl0tstrib— density!  •), 50,80  >  !  Find  nex  value 

Faetor-I00.7f1ax_dist  (  Scale  for  100  nax 

FOR  I-l  TO  58 
FOR  J-l  TO  80 

Distrib_density!  I , J )-0tstrib  density! I , J  )»Factor  !  Scale  values 
NEXT  J 
NEXT  I 

FOR  I-l  TO  50  !  Ness  density/unit  arse  at  height  above  bed 

Uolune-PI»! ! 1+7 )/ 1000 )*3/S  !  Wolune  Cu  cn 

FOR  J-l  TO  80 

Entrainnent!  I ,  J  (-Entrainnent !  1 ,  J  )+Uolune-»Qistrlb_densi  ty!  I ,  J  ) 

NEXT  J 
NEXT  I 

rtax_antrain-l 

FOR  I-l  TO  30  i  Find  naxinun 

IF  Entrainnent!  1 , 1  )>Entrainnent!  1  .Hex  entrain)  THEN  flex  antretn-I 
NEXT  I 

Factor-100/Entrainnent!  I  .flax^antratn  >  !  Normalise  to  100 

FOR  I-l  TO  80 

Entrainnent! I , I )-Sntrainnent! 1 , I (‘Factor 
IF  Entrainnent! I, I  )<«0.  THEN  1S68 
Entrainnent! 2 , 1 )-LQ6\ Entrainnent! I , I > ) 

NEXT  I 

i 

1  Outsat  Central  Section 
i 

PRINTER  IS  1  1  Output  to  CRT 

PRINT  USING  *9.3/" 

PRINT  "1)  Oisolay  height  vs  oianeter  data' 

PRINT  *2)  Oisolay  height  vs  dianeter  grsoh' 

PRINT  "3;  Oisolay  density  vs  heiqnt  as  function  or  dia  graon" 

PRINT  *i  Oisolay  density  vs  dianeter  as  a  function  of  heignt  arson’ 


Pert  Jn»  i ght  -P  ar  t _ha  i gh t  +!  Va  1  oe  i  t  y +Upo  ) •( 
IF  Uoo>0  THEN 
Heignt ! 1 . I )-Pert_heignt 
Heignt! 2, 1  )-Heignt!2, 1 5+1 
ENO  IF 

Upo-Uelocity 

Point-! 100*Pert_heignt  )  QIU  2 
IF  Pomt>73  THEN 
Point-73 
ENO  IF 

IF  ?oint<-0  THEN 
Point-0 
ENO  IF 

PRINT  I, Point 


1769  PRINT  "S )  Oisoley  sens  as  3  but  with  bed  density* 

1779  PRINT  *3)  Oisoley  sene  es  i  but  with  bed  density* 

1789  PRINT  *7)  Oisoley  density  vs  dieneter  of  bed  ness* 

1739  PRINT  *3 )  Oisoley  sntrainnent  density  ebove  bed* 

1899  PRINT  *3)  Oisoley  Ln  sntreinnent  density  above  bed* 

1819  PRINT  *19)  EXIT  PR06SAH" 

1829  INPUT  ‘Enter  minder  of  desired  display:  *, Answer 
1839  SELECT  Answer 

18*9  CASE  >1 

13S9  CALL  Oisolay_data<UI ,Uoi ,Uj , Jet  t.Oel  t, Hones, Clear*) 

1869  CASE  *2 

1879  CALL  Oisnlav  or aoh< Clear S.HonaS- Answer. View) 

1389  CASE  *3 

1839  INPUT  'Enter  particle  sice  to  be  viewed  <88-378  un)(9  for  all):  *,Uie 

1399  IF  U lew-9  THEN  1329 

1319  IF  Uiew<39  OR  Utew>S79  THEN  1839 

1329  Uiew-(Uisw  OIU  18 )-7 

1338  CALL  Oisoley  ar aoh< Clears. Honei. Answer .View) 

13*9  CASE  -4 

13S9  INPUT  'Enter  desired  height  above  bed  surface  (8-158  cn>:  '.View 

1369  IF  Ulew<8  OR  Uiew>153  THEN  1359 

1379  View-< View  OIV  2»l 

1389  CALL  Oisolay_graph<Cleerd,HoneS,Answer , View  ) 

1338  CASE  -5 

2909  INPUT  "Enter  particle  size  to  be  viewed  <80-379  un>(9  for  all):  *,Vle 

2810  IF  V lew-0  THEN  2039 

2020  IF  Vlew<80  OR  Vlew>570  THEN  2900 

2030  Vlew-< View  OIV  10 )-7 

2040  CALL  Oisoley  araoh<ClaarS.Hone*. Answer .View > 

2050  CASE  -6  '  ~ 

2969  INPUT  ‘Enter  desired  height  above  bed  surface  (0-158  cn>:  ".View 

2070  IF  Uiew<0  OR  ViewMSS  THEN  2060 

2080  View-< View  OIV  2)4-1 

2030  CALL  Oispiay_gra0h<ClearS,HoneS,Answer, View  ) 

2100  CASE  -7 

2110  CALL  Oisolay_gragh<ClearS, HoneS, Answer, View) 

2120  CASE  -8 

2130  CALL  Qisolav  oraoh( Clears .Hone* .Answer .View ) 

2140  CASE  -3 

2159  CALL  Oisoiay_graoh(CIearS,HoneS,Answer,Viaw ) 

2160  CASE  -19 

2170  STOP 

2189  CASE  ELSE 

2139  SOTO  1829 

2200  ENO  SELECT 
2210  SOTO  1700 
2229  ENO 

2239  ! 

2240  ! 

2250  « 

22S3  !  Sue  used  to  list  incut  parameters  and  max  heignt  at  tine  t  car  dianeter 

2270  ILS  ji sc lay_data<  U1  .Uol  ,'Jj  ,  Jet_t  .Cei.t , HoneS, Clears  ) 

22S0  CCA  IN7SSE?.  Qistr ibut ion(  •  )  ,0isiri3_cena i ty<  •  1  ,?.E.-,L  -e ignt  <  *  >  ,3eo_part  ic  1 
< •  )  .Entrainment ( •  ) 

2290  OUTPUT  2:HomeS;  i  Home  and  Clear  screen 

2390  OUTPUT  CiClearSi 
2310  • 

2320  PRINT  ' 

2330  PRINT  "Mean  9ed  Velocity- 


cm/ s " 

. Ul *30 . 49 


Print  outgut 


2358  PRINT  'Peak  J«x  Velocity-  *;UJ*30.48 

2360  PRINT  *Sas  Jet  Ouration-  *;J«t  t 

2378  PRINT 

2388  PRINT  *  Oianetar*,*  Haight", *  Tina*,*  Oianatar*;*  Haight*:* 

Tina* 

2398  PRINT  *  «"*,*  cn  *,*  seconds*,*  un * , *  cn",*  seconds 

2488  PRINT 

2418  FOR  >1  TO  58  STEP  2 

2420  PRINT  USIN6  *2(5X ,20, SX, 30.0, SX, 0.000, SX  )*: < 1+7 )• 10, Heights  1, 1 >, Height (2 

, I  )*0al_t , ( 1+8  >*10,Heigni< 1 , 1+1 > .Height (2 . 1+1 >*0el_t 

2430  NEXT  I 

1448  PRINTER  IS  1 

2458  INPUT  'Print  hard  cooy?  < 1 )■  yes,  <0>»  no  =  *, Answer 

2460  IF  Answer- 1  THEN 

2478  PRINTER  IS  701 

2488  60T0  2320 

2490  ELSE 

2508  IF  Answer<>0  THEN  2458 
2510  ENO  IF 
2528  SUBENO 
2538  ! 

2548  ! 

2558  I 

2560  !  Sub  used  to  control  graohics  outout  of  data 

2570  SU8  Oisolav  or aoh< Clears. Honed. INTESER  Oata.set .View ) 

2588  COM  .INTESER  Oistr ibuiion< ») ,Olstrlb_density< •>, REAL  Heightt • >,8ed_particie 
( • )  ,Entrainnent( • ) 

2530  ' REAL'Xnax , Ynax ,  X  t  ick ,  Yt ick , Xnln, Ynin 
2608  OUTPUT  2: Honed; 

2618  OUTPUT  2:Clear*; 

2620  6INIT  !  Initialize  graohics 

2530  SRAPHICS  ON 

2640  SELECT  Oata.set  !  Scale  plot  routines 

2658  CASE  -2 
2668  Xnax-608. 

2670  Ynax-FNMax_raal<Height<  *>,1,S0> 

2688  X tick-58. 

2698  Yt ick-FNScal*< Ynax  ) 

2700  CASE  -3 

2710  Ynax-FNNax.int (Oistr ibut ion< •  ), 50,30 ) 

2720  Xnax-FNGaia_linit<Oistribution( * ),50,30 ) 

2730  Xtick-FNScaie<Xnex ) 

2740  Yt ick-FNScale< Ynax  ) 

2750  CASE  -4 
2768  Xnax-500 . 

2770  Ynax«FNHax_int(Oistr ibut ion< • >,50,30  > 

2780  X tick-50. 

2790  YUcfc-FNScale<Ynax  > 

2800  CASE  *5 

2810  Xnax-rNOata_linit<Oistrib_density< ♦ ) ,50,30 ) 

2920  Ymax-100. 

2920  X t ic*--NScaie< Xnax ) 

2340  Yt icK-rNScsl s: Ynax ) 

2958  CASE  -5 
2960  Xnax-500 . 

2870  Ynax-FNhax_int ( Oistr ib_aensi ty< • ) ,50,30 ) 

2880  Xtick-50. 

2390  Yt ick-FNScale< Ynax ) 

2900  CASE  -7 


2310 

Xnax-600. 

2320 

Ynax-100. 

2330 

Xtick-50. 

2340 

Ytick-10 

23S0 

CASE  -8 

2350 

Xnax-FNOata_l  ini tr<  Entrainnent ( * >,1,30) 

2370 

Ynax-100. 

2380 

Xt  ick-FNSca ia( Xnex  ) 

2330 

Ytick-10. 

3000 

CASE  -3 

3010 

Xnax-FNOata_linitr(Sntrainnent( ♦  ), 1,30  > 

3020 

Ynax-5. 

3030 

Ynxn— 3. 

3040 

Xtick-FNScale(Xnax  ) 

30S0 

Ytick-.S 

3060 

Y_axis— 2.5 

3070 

60T0  3110 

3080 

ENO  SELECT 

3030 

Y_axis-0. 

3100 

Ynin— Ytick 

3110 

Xmn— 2.  »Xtick 

3120 

X_axis-0. 

3130 

UINOOU  Xnin,l.l*Xnax,Ynin, 1. l»Ynax 

!  Set  graohics  scale 

3140 

AXES  Xtick, Ytick, X_axis,Y_axis 

<  Set  graoh  axis 

3  ISO 

L0R6  S 

!  Label  X-axis 

3160 

FOR  I-X_axis  TO  Xnax  STEP  Xtiefc-2 

3170 

HOWE  I,Y_axis 

3180 

LABEL  I 

3130 

NEXT  I 

3200 

L0R6  a 

!  Label  Y-axia 

3210 

FOR  I-Y_axla  TO  Ynax  STEP  Ytick 

3220 

HOWE  X_axis,I 

3230 

LABEL  I 

3240 

NEXT  I 

32S0 

SELECT  Oata.set 

!  Plot  data 

3260 

CASE  -2 

3270 

HOWE  80, Height ( 1,1) 

3230 

FOR  1-2  TO  50 

3230 

ORAU  ( 1+7 )• 10, Height (1,1 > 

3300 

NEXT  I 

3310 

CASE  -3 

3320 

Begin-Wieu 

3330 

Fimsh-Wlew 

3340 

IF  Wlew<0  THEN 

!  Plot  all  if<0 

3350 

Segm-1 

3360 

Fjniah-50 

3370 

ENO  IF 

3380 

FOR  I-8egin  TO  Finish 

3330 

HOWE  0, Distribution* I, 1 ) 

3400 

FOR  J-l  TO  Xnax/2 

3413 

ORAU  2-J, Distribution! I.J  ) 

3423 

NEXT  ; 

3428 

NEXT  : 

3440 

CASE  «4 

34S0 

HOWE  30  .Oistnbut  ion(  1  ,Wlau  ) 

34S0 

FOP.  1-2  TO  50 

3470 

ORAU  ( I+7)»l0,Olstri bution( I, View  > 

3480 

NEXT  I 

3490 

CASE  -5 

3530 

Segm-Wieu 

Ztia 

!  Plot  all  if  <0 


3510  Fmish-Ute« 

3520  IF  Uiew<0  THEN 

3550  3egin-l 

3540  Finish-50 

3550  EnQ  IF 

5550  FOR  I-8eqin  TO  Finish 

3570  nave  0,Oistrib_density<I,l> 

3580  FOR  J-l  TO  Xmex/2 

3550  J)RAU  J*J^Qiatrib_d»n*ity( I ,J  ) 

3600  NEXT  J 

3610  NEXT  I 

3620  CASE  -5 

3630  HOME  90,0iatrib_density( 1 ,Uleu > 

3640  FOR  1-2  TO  50 

36S0  ORAU  <I+7>»10,Oistrib_density<I,Uiew) 

3680  NEXT  I 

3670  CASE  -7 

3630  MOVE  30,3ed_particle< 1 ) 

3690  FOR  1-2  TO  50 

3700  ORAU  ( I+7)*l0,3ad_jJartici*< I  ) 

3710  NEXT  I 

3720  CASE  -8 

3730  MOVE  0, Entrainment! 1,1 1 

3740  FOR  1-2  TO  80 

3750  ORAU  I »2, Entrainment! 1,1 > 

3760  NEXT  I 

3770  CASE  -9 

3780  MOVE  0, Entrainment! 2, 1 ) 

3750  FOR  1-2  TO  30 

3800  IF  Entrain*»nt<2.n-0.  THEN  3840 

3810  ORAU  1-2, Entrainment!  2, 1 ) 

3820  NEXT  I 

3830  ENO  SELECT 

3840  INPUT  “Print  graoh  (1)  yea,  (0)  no:  *, Answer 
3850  IF  Answer- I  THEN 
3860  QUMP  QEUICE  IS  701  '  Ouno  to  printer 

3870  DUMP  6RAPHICS 

3880  5CLEAR 

3890  ELSE 

3900  IF  Anawer<>0  THEN  3840 

3910  ENO  IF 

3920  SCLEAR 

3930  SU8ENO 

3940  ! 

3950  ! 

3550  < 

3970  !  Function  determines  graorucs  axis  acaleing 

3980  QEF  rNScale! Ymax ) 

3990  IF  Ymax> 100  THEN 
4000  Tick-20. 

4010  else 

4020  IF  Yn«x<-13C  ANO  Ynax>23  THEN 

4030  Tick-10. 

4040  ELSE 

4050  IF  Ymax<-20  ANQ  Ynax>4  THEN 

4060  Tick-2. 

4070  ELSE 

4080  Tick«.S 

4050  ENO  IF 

4100  END  IF 


245 


4110  END  IF 
4120  RETURN  Tick 
4130  FNENQ 
4140  ! 

41S0  ! 

4160  ( 

4170  !  Function  deterninea  nax  extent  of  data  for  X-axis  Unit 

4180  OEF  FNOata_iinit( INTE6ER  Oata_l<« ),Row_nax,Col_nax ) 

4130  Col-1 
4200  REPEAT 
4210  Sun-0. 

4220  FOR  1*1  TO  Rou_nax 

4230  Sum-5um'rOeta_l(  I, Cal  > 

4240  NEXT  I 

4250  Col-Col+t 

4260  UNTIL  Sun-0  OR  Coi-Col_nax-l 
4270  RETURN  CoI*2 
4230  FNENO 
4250  ! 

4300  I 
4310  I 

4320  !  Function  deterninea  nax  value  in  integer  array 

4330  OEF  FNflax_int(  INTESER  Oata_l<  *  ),Rou_nax,Col_nax  ) 

4340  R_nax-1 

4350  C_nax-1 

4360  FOR  1*1  TO  Row_nax 

4370  FOR  J-l  TO  Col_nax 

4380  -  —  IF  Oata_i< I.J )>0ate_l( R_nax , C_nax >  THEN 

4380  R_nax*I 

4400  C_nax-J 

4410  ENO  IF 

4420  NEXT  J 

4430  NEXT  I 

4440  Nax*0ata_l <  R_nax , C_nax ) 

4450  RETURN  Nax 
4460  FNENO 
4470  j 
4480  I 
4430  ! 

4500  !  Function  deterninea  naxinun  value  in  real  array 

4510  OEF  FNMax_rsal(Oata_l< ♦), INTESER  Rou_nax , Col _max ) 

4520  R_nax-l 

4530  C_nax-1 

4S40  FOR  I-l  TO  Rou_nax 

4530  FOR  J-l  TO  Col .max 

4550  IF  Oata_l( I. J )>Oata_l(R_nax,C_nax )  THEN 

4570  R_nax”I 

4530  C_nax-J 

4530  ENO  IF 

4630  NEXT  J 

4610  NEXT  I 

4623  flax  -Oat  a_l  ( R_nax,  C_nax  ) 

4630  RETURN  flax 
4640  FNENO 
4653  i 
4650  ! 

4670  ! 

4630  1  Function  deterninea  nax  extent  of  data  for  X-axis  Unit 

4630  OEF  FNOaia_l ini tr< Data_l< • ) ,Rou, Col_nax ) 

4700  Col -l 


Z4B 


ftPPENDIX  M 


This  Appendix  contains  the  calibration  data  for  the 
anemometer  probe.  The  calibration  uas  conducted  in  a  small  uind 
tunnel  using  a  pitot  tube  connected  to  a  micromanometer  capable  of 
measuring  pressures  to  uithin  0.001  ins.  of  water. 


END 

FILMED 

11-85 

DTIC 


